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ABSTRACT 



The subjects of this thesis are the invisible axion and the more general family of axion- 
like particles. 

The invisible axion is a hypothetical elementary particle and a cold dark matter can- 
didate. I present an improved computation of the constraints on the parameter space 
of the cold dark matter axion in the standard cosmology, that includes the contributions 
from anharmonicities in the axion potential and from the decay of axionic strings. In this 
scenario, I update the value of the mass of the cold dark matter axion, finding the value 
(67 lb 17)^eV, approximately one order of magnitude larger than previous computations. 

The effect of nonstandard cosmological scenarios on the parameter space of axion cold 
dark matter is studied for the first time. In particular, I consider the cases of low- 
temperature reheating and kination cosmologies, and I show that the mass of the cold 
dark matter axion can differ from the value in the standard cosmological scenario by orders 
of magnitude. 

Finally, I consider the family of axion-like particles, assuming that these particles serve 
as the infiaton in the context of warm inflation. I find that the axion energy scale /, which in 
the standard inflation scenario is of the order of the Planck mass, can be lowered to the much 
safer Grand Uniflcation Theory scale / ~ lO^^GeV. I also constrain the parameter space 
and the amount of gravitational waves from this model, using results from the Wilkinson 
Microwave Anisotropy Probe 7-year data. 



A mia moglie Erika 
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CHAPTER 1 



INTRODUCTION 

My interest in the invisible axion, a hypothetical pseudo-scalar particle arising in the 
QCD sector of the Standard Model, addresses the fact that this particle might play the role 
of the Cold Dark Matter (CDM). This possibility has been intensively studied for more than 
thirty years; original work on the subject by myself and my advisor prof. Paolo Gondolo is 
exposed in Chapters [5] and [H drawn from Refs. |185| and |186j . respectively. 

In Chapter [SJ I considered the invisible axion as the CDM particle in the standard 
cosmological scenario, updating the bounds on the axion parameter space in light of the 
most recent cosmological data, and improving the treatment of anharmonicities in the axion 
potential and axionic strings. In Chapter [6l it is shown how the invisible axion can be used 
as a probe to study the properties of the universe before Big Bang nucleosynthesis (BBN) 
took place. In fact, we only know that the universe has to be radiation dominated (standard 
cosmology) at the time of BBN, but a plethora of different possibilities could have taken 
place before this period. We studied the properties of the axion CDM in two different 
nonstandard cosmological scenarios, the low temperature reheating (LTR) cosmology and 
the kination cosmology. 

The last chapter of the thesis deals with the role of axion-like particles, a generalization 
of the invisible axion theory, as the inflaton in the specific scenario of warm inflation, and 
refers to the paper in Ref. |184j . Axion-like particles have long been successfully considered 
as the inflaton particle in the standard inflaton, with the model going under the name of 
Natural Inflation (NI). One problem of the NI model resides in that the energy scale at 
which the axion symmetry spontaneously breaks is of the order of the Planck mass, which 
makes it difficult to embed NI into larger Grand Unification Theories (GUT). In Chapter [71 
I show that this problem can be overcome if we consider the warm inflation scenario in 
place of standard inflation, thus developing the natural warm inflation model (NWI). 



CHAPTER 2 



CONCORDANCE COSMOLOGY 

A quantitative approach to cosmology began in the early 20th century, after Einstein's 
General Relativity (GR) was formulated. In fact, cosmology as we know today is firmly 
established on the basis of the GR theory, as we will explore in detail in Sec. 12. 2[ Before 
relativity, space and time were considered fixed, an idea that directly followed from the New- 
tonian concept of an absolute framework in which physical processes would take place. The 
birth of GR provided scientists with a new concept, the fact that space-time is a dynamic 
entity whose evolution is described by Einstein equation. Within the GR theory, gravity 
can be understood as a geometric property of the space-time. Soon after the introduction 
of GR in its final version, the scientific community began speculating on the origin of the 
universe, merging observations with the new framework provided by the GR theory. The 
cosmological solutions of general relativity were first found, independently, by A. Friedmann, 
G. Lemaitre, H. P. Robertson, and A. G. Walker and constitute what is known as the 
Friedmann- Robertson- Walker (FRW) model (see Refs. [1081 1192j for exhaustive reviews on 
the subject). 

2.1 General relativity 

The insight behind GR consists of equating the effects of the gravitational force on a 
test particle with the acceleration that such a test particle would experience on a nontrivial 
manifold. This idea is encoded in the weak equivalence principle, which states that local 
effects of motion in a curved space-time are indistinguishable from those of an accelerated 
observer in a flat space-time, and it is depicted by the famous elevator thought-experiment. 
With this key concept in mind, Einstein derived the famous formula 

n^,^- g^^n = S^GT^^+g^^K. (2.1) 

The term A appearing in Eq. (|2.ip is a constant that might play an important role in the 
description of the universe at present times. We will discuss the effects of this term in depth 
in Sees. 12.2.2.21 and WIM 



We first discuss the left-hand side of Eq. (j2.ip . Here, g^y is the metric of the manifold 
considered, from which we can derive the Ricci tensor IZ^i, via the Riemann tensor T^JJ^p as 

n^u = n%,- (2.2) 

The Riemann tensor 71"^ is defined through the Christoffel symbol (Levi-Civita connection) 



as 



Here, the Levi-Civita connection is given by 

^% = \ f [du gp^. + g,p - dp gp,) , (2.4) 

and dp means differentiation with respect to the manifold coordinate x^. 

In Eq. (12. ip , the curvature is 7^ = g^^ 'R-pu , with the inverse of the metric g^'^ satisfying 

gpag""" = 6';,. (2.5) 



Sometimes the combination 



Gpu = TZpu - gpu Tl, (2.6) 



is referred to as the Einstein tensor. 

On the right-hand-side of Eq. (|2.ip . the term Tp^j describes the energy- momentum 
content at each point of space-time and is known as the stress-energy tensor. Since Tp^, 
depends implicitly on the metric tensor, in GR a feedback mechanism appears in which the 
metric structure of space-time and the matter content of the universe mutually influence 
each other, a statement that can be summed up with the paraphrase by Wheeler: space-time 
tells matter how to move; matter tells space-time how to curve. In the following, we 
will always use the case for an isotropic perfect fluid of density field p, pressure p and 
four- velocity Up, whose stress tensor reads 

Tpy = pUpUy +p{upUy - gpy). (2.7) 



We remark that the density p and the pressure p are time-dependent quantities, although 
in the text, we do not indicate such dependence explicitly. In Sec. \2.2.2\ we will discuss how 
the density and the pressure for a specific fiuid are related by a specific equation of state. 



2.2 General relativity and cosmology 
2.2.1 The Priedmann- Robertson- Walker metric 

In the FRW metric, the Copernican principle that demotes Earth's position in the Solar 
system is extended in a cosmological framework, by stating that no point in the universe is 
special with respect to the other ones. This concept, called the cosmological principle, can 
be restated by posing that the universe is isotropic with respect to any point, or equivalently, 
that it is homogeneous and isotropic with respect to one point. Isotropy and homogeneity 
can be understood in terms of cosmic inflation, see Sec. 12.31 and lead to the FRW universe, 
with line element 



+ d9'^ + sin^ i 



1 - kr^ 



(2. 



Here a{t) is the scale factor of the universe depending on the cosmic time t, and the term 
in square brackets indicates the spatial metric, which is specialized to represent a three- 
manifold of constant curvature k. The only possibilities for a spatial manifold of constant 
curvature are zero curvature with k = (flat universe), positive curvature with k = +1 
(closed universe), and negative curvature with k = —1 (open universe). Sometimes the 
conformal time is defined, 

dr = (2.9) 
a{t) 

so that the FRW metric is written in the conformal form, 

ds^ = aHt{T)) (dr' - -r^de^ -r' sin'edA . (2.10) 

\ 1 — /e J 

2.2.2 The Priedmann equations 

The Friedmann equations are a set of two equations that relate the expansion of the 
universe to its energy density and pressure content. In this section, we expose their 
derivation from the Einstein equation and some of their uses. 

In the formulas used for computations, we trade Newton's gravitational constant G, 
appearing in Eq. (|2.1|) . with the Planck mass Afpi, using the relation 

G = ^. (2.11) 



Mi, 



For numerical computations, we will use Mpi = 1.221 x 10^^ GeV 



2.2.2.1 Derivation of the Friedmann equations 

It can be shown that the 00 and 11 components of the Einstein Eq. (j2.ip . expressed on 
the FRW metric in Eq. ()2.8p and with the source in Eq. ()2.7p . give the relations 

and 



a(t) 47r , , A , , 

A dot over a quantity will always indicate a total derivative with respect to the cosmic time 
t, so a{t) = da{t)/dt. Because of isotropy, the 22 and 33 components of Eq. (j2.ip give the 
same expression as Eq. (I2.13p . while the off-diagonal terms are all zero due to homogeneity. 
The role of the term A, also known as the vacuum energy term, will be discussed in detail 
in Sec. 12.2.2.21 

In a realistic cosmological model, more types of fluids are present. Examples include the 
nonrelativistic matter fluid, of density pm and pressure pm = 0, and the radiation fluid, of 
density pr and pressure pr = Prl'^- When more than one fluid are present in the theory, we 
define the total energy densities ptot and the total pressure ptot of the system as 

Ptot = X] Ptot = X] '^^■^^ 

i i 

where pi and pi are the density and the pressure of the i-th fluid. To completely specify the 
system, we also need an equation of state linking the density and the pressure of each fluid, 

Pi = WiPi, (2.15) 

with the parameter Wi specifying the type of fluid. For example, Wm = is for a nonrela- 
tivistic fluid and Wj. = 1/3 is for a relativistic fluid. 

In the multiple fields case, we substitute p — ptot and p ptot in the Friedmann 
Eqs. (I2.12P and (j2.13p . This is the case we will discuss from now on. Defining the Hubble 
rate 

Hit) = ^ (2.16) 
a{t) 

we obtain the Friedmann Eqs. ()2.12p and (j2.13p in the form 
and 

ptot = -3H{t) [ptot + Ptot). (2.18) 
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Unless otherwise specified, we will always derive our results in the next chapters in the 
case of flat geometry = and negligible cosmological constant A = 0, which are suitable 
approximations in the cases we consider. For these reasons, we will refer to the Friedmann 
equation in the form 

We will keep the parameters k and A for the rest of the discussion in the current chapter. 

2.2.2.2 The term A 

In cosmology, the constant term A is usually identified with the vacuum energy, respon- 
sible for the accelerated expansion we observe at the present time. Using Eq. (I2.12p . we can 
associate to the cosmological constant A the energy density and the pressure 

PA = —^^ PA = -PA. 2.20 
Stt 

From now on, we include, in the definition of the terms ptot and ptot, the contributions from 
p\ and PA) respectively. With these definitions, the Friedmann Eq. (j2.17p is rewritten as 

From Eq. (j2.15p . we see that the equation of state for the cosmological constant term has 
w\ = — 1. Models in which only a cosmological constant is present have been studied since 
the early days of relativistic cosmology, and fall under the name of Lemaitre models. 

To sum up, in order to include the contribution from the vacuum energy in the Friedmann 
equations, one has to add the energy density pA to the sum ptot, and similarly the pressure 
Pa to the sum ptot, keeping in mind that the equation of state for the vacuum energy is 
PA = -PA- 

2.2.2.3 Equivalent forms of the Friedmann equations 

The Friedmann Eq. (j2.2ip is particularly suitable for solving specific problems in which 
the content of the universe is specified and for explaining the thermal history of the universe. 
Defining the critical density. 

Pent = (2.22) 

and the cosmological density ratios. 

Pi o u\ /'tot 



n,{t) = -^, ntot{t) = ^ = J2ni{t), (2.23) 

Pcrit Pcrit 



we cast the Friedmann Eq. (|2.17p in the alternative form 

J^^=H\t){QMt)-l)- (2.24) 
Eq. ()2.24p shows that the universe is 

• open (fe = — 1) when r2tot(i) < Ij 

• flat {k = 0) when Qtot{t) = !> 

• closed (k = +1) when ^Itotit) > 1- 

The second, alternative formulation of the Priedmann equation we show now is partic- 
ularly suitable for computing the scale factor a{t) given the energy content of the universe 
at a specific time, like for example at the present time to- We indicate the values of the 
density ratios at the present time with the index 0, and the present value of the critical 
density is 

Pc = Pcrito = ^^kMli = 1.878 X 10-29 g/cm3, (2.25) 
svr 

where the Hubble expansion rate at present time and the value of h are |110] 

Hq = lOO/ikms-^Mpc^^ = (70.2 ± 1.4) kms"^ Mpc'^ (2.26) 
Introducing the density of curvature 

Eq. (|2.17p is cast in the form 

^ = ^r(^y + n„j^)\nj^)\n^ + nj^) ''^'''^^ . (2.28) 



Hi \a{t)J "'\a{t)J \a{t)J Wt) J 

Here and in the following, we have used the notation fij = r2j(to)- We have included the 
possibility to include in the theory a fluid of density ratio 0^ with generic equation of state 
p = w p and w unspecified. Notice that setting a{t) = a(to) = clq in Eq. (j2.28p gives the 
constraint 

Qr + ^ra + ^k + ^A + ^w = 1- (2.29) 



2.2.3 Redshifts 

The redshift z describes the lowering in the frequency of cosmological photons that have 
been emitted in the past with some frequency v and are received at the present time with 
frequency fo, and is defined as 

l + z = — . (2.30) 

Since it can be shown that the momentum of free particles moving on the FRW metric 
decreases with l/a(t), as well as lengths and wavelength stretch proportionally to a(t), the 
redshift at some time t is related to the scale factor at that time by 

l + z = ^. (2.31) 

ait) 

A precise relation between redshift z and time t is given below in Eq. (I2.39p . Since the scale 
factor a(t) — )• for t — )• (Big Bang) and a(t) — )■ oq at present time, we also have z — )• +oo 
when t — 7- and z = at the present time. 

Using Eq. (j2.3ip . one can trade the time or scale factor dependence appearing in some 
equation with a z-dependence. In this view, a useful form of Eq. (|2.28p that provides the 
Hubble expansion rate in terms of the redshift z is 

R{z) = Ho klr (1 + zY + (1 + zf + Ofc (1 + zf + + 0^ (1 + zf^^+'^A ^'^ . (2.32) 

2.2.4 Particular solutions of the Priedmann equation 

We show the explicit solution of Eq. (I2.19P in various cosmologically interesting scenarios. 
In general, we expect the energy density pi that dominates the universe at some particular 
time to depend on the scale factor as 



(^)^^^"^^ (2.33) 
The Friedmann Eq. (j2.19p with p = pi is written in the form 



Pi = PiO 

' a 



with solution 



Exp [//oi\/^] 1 with Wi = —1. 
In particular, the solutions when the universe is dominated by matter (wm = 0), radiation 
{wr = 1/3), and vacuum {wa = —1) are, respectively, 

a{t) DC t^^^ for matter domination, (2.36) 
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a(t) oc t^/^ for radiation domination, (2.37) 

a(t)oce*^^ for vacuum domination. (2.38) 

Notice that when Wi = —1, the scale factor expands exponentiahy. A family of cosmological 
models that use such accelerated expansion in the very early stages of the universe has been 
used in order to solve serious problems linked with the standard Big Bang model and fall 
under the name of inflation. We will discuss these problems and the solution posed by the 
inflationary scenarios in Sec. 12.31 

When more than one fluid is present, a general expression for the cosmological time t as 
a function of the redshift z can be found by integrating Eq. (|2.32p , 

we derived Eq. (j2.39p setting our initial time t = at redshift z = +oo. The age of the 
universe to is given by Eq. (|2.39|) by setting 2: = 0, so 

For a universe where only matter is present we have 



Ho to 



Ar- 2(i-">/^ ArcCosh(^-l), forO„,<l, 
2 

3' 

+ ArcCos (t^-1], for f]^ > 1. 



for Qm = 1, (2.41) 



This equation, for = 1 and with i^o = 70 km/s/Mpc yields to = 9.3 Gyr, which 
underestimates the effective age of the universe by almost 35%. To obtain the correct result 
for to) we need to include the dark energy contribution in Eq. (j2.39p : in fact, considering 
a flat (k = 0) universe where matter ilrn and vacuum energy 1 — ilrn are present (with 
negligible radiation), we obtain 



Ho to = — ■j^=^= In 



2(1 + yi - Qm) ^ 



(2.42) 



Eq. (12.421) with Ho = 70 km/s/Mpc and Qm = 0.25 gives the correct result for the age of the 
universe, to = 13.7 Gyr. In Sec. 12. 6. 3^ we will discuss how current observations motivate a 
nonzero value of ^a- 
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2.3 Cosmic inflation 

The standard Big Bang theory yields striking successes in explaining a large number of 
cosmological observations. However, in order for this model to actually be consistent, the 
universe has to emerge from the Big Bang with very specific initial conditions, in order to 
match measurements of the quantities we observe at the present time. The most challenging 
problems that the standard Big Bang theory faces are the horizon, flatness, and unwanted 
relics problems. The inflationary scenario was introduced in order to solve these specific 
problems [102^ 1165^ [78| 1151^ EJ I118j . and it also proved to be a valid model in which the 
seed of local inhomogeneities form [1371 [79l [571 11661 118j . For reviews of the inflationary 
mechanism, see Refs. [1201 1108] . 

2.3.1 The flatness problem 

If the initial value of the total energy density of the universe slightly differs from the 
critical density, or |r2tot(i) — 1| 7^ in Eq. ()2.24p . any deviation from unity will eventually 
be amplified by the present time; alternatively, it is required that the universe be extremely 
flat at early times, in order to explain the close-to-flat geometry we observe today. 

To illustrate the flatness problem in cosmology, we assume that the universe had an 
initial deviation from unity r^tot.i — 1- Indicating with Otot = ^^tot(io) the value of Vttot{t) at 
the present time Iq, we compute the deviation of Otot from unity by considering Eqs. (|2.24p 
in the form k = a? (Otot,i — 1): assuming radiation domination at all times a{t) ~ t^/^, we 
obtain 

|Otot,i - 1| = l^ot - 1| , (2.43) 

where tpi = 1/Mpi = 5.3 x 10"'^^ s is the Planck time and to = 13.7 Gyr is the age of the 
universe. Numerically it is 

|J^tot,i-l| « 10^^°|^^tot-l|. (2.44) 

Measurements of the anisotropies in the Cosmic Microwave Background Radiation (CMBR) 
by the WMAP7+BA0+SN data constraint |Otot - 1| at 95% C.L. as [IIO] 

- 0.0178 < S^tot - 1 < 0.0063; (2.45) 

we state that |iltot — 1| ^ 0.01, implying that iltot,i initially differed from one by one part over 
10^^, any initial deviation of ritot.i from unity being amplified at the present time according 
to Eq. ()2.43p above. Summing up, we have that the primordial universe has to exit from 
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the Big Bang being extremely flat in order to explain the tiny anisotropy in the density 
observed. Inflation solves this problem because, during the inflationary period, the value 
on the right hand side of Eq. (I2.24P drops exponentially due to the exponential increase of 
the scale factor. Inflation has to last sufficiently long in order to solve the flatness problem: 
a parameter used to describe such a requirement is the number of e- folds A'e, defined as the 
logarithm of the ratio of the scale factors at the time when inflation ends tend and at the 
beginning of inflation ti , 



The horizon problem deals with the fact that a large degree of homogeneity is observed in 
the sky, although most of the patches have never been in causal contact before. For example, 
inhomogeneities in the CMBR are observed only at a 6T/T ~ 10^^ scale, even though the 
sky contains about a hundred patches that never interact causally in the standard picture. 
A very simplified view of this problem can be sketched by computing the age of the universe 
using Eq.([232]) with i^o = 70 km/s/Mpc and = 0.25, which yields to = 20 Gyr. This 
computation overestimates the effective age of the universe to = 14.3 Gyr and implies that 
in the past, there was not enough time for the universe to be in causal contact. An early 
inflationary period would solve such a problem because in this scenario, an initial region 
that was in causal contact before inflation is then expanded and appears at the present time 
as distinct patches not causally connected. 



Generic Grand Uniflcation Theories (GUT) predict heavy particles and topological 
defects to be copiously produced in the early universe (see Ref. |182| for a review), to 
the extent that these relics and the particles emitted by them would eventually dominate 
the expansion of the universe [22]. Since such relics are not observed in the present universe, 
it is believed that some other mechanisms have diluted their number to the present day. One 
viable explanation for this fact is that defects were produced before or during inflation, with 
monopoles being separated from each other (or "washed out") as the accelerated expansion 
progressed and their number density being consistently reduced to a safe cosmological value. 




(2.46) 



Sufficient inflation requires iVg > 60. 



2.3.2 The horizon problem 



2.3.3 Unwanted rehcs 
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2.3.4 Small-scale structures 



Since the standard FRW model describes a homogeneous and isotropic universe, it does 
not account for the structures we observe today at various scales, like stars, galaxies, and 
clusters of galaxies. For this reason, perturbations in the FRW metric have been long 
studied, see Refs. pTlll06| . and following the evolution of these perturbations, it is possible 
to explain the spectrum of inhomogeneities we observe in the CMBR. The initial power 
spectrum of perturbations that describes such variations in the densities will be reviewed 
in Sec. El 



A period of inflation is defined as a period of accelerated expansion: using Eq. (j2.13p 
with p = w p, it follows that it must be 



If the condition in Eq. ()2.47p is satisfied, a solution to the Friedmann Eq. (j2.12p when k = 
and p is constant is (see Eq. (I2.38P above) 



An exponential growth like in Eq. ()2.48p can be obtained from the Friedmann equation 
whenever the density, and thus the Hubble, are constant. This situation can be achieved in a 
A-dominated universe, the so-called Lemaitre models, or when we have some form of energy 
that is dominating the expansion rate of the universe whose energy density p is constant. 
It is generally believed that a constant term A might account for the accelerated expansion 
at the present time, see Sec. 12.6.31 whereas in order to give a microscopic explanation of 
primordial inflation, the early domination of an exotic form of energy has been invoked. It 
is usually postulated that the expansion rate of the universe during the inflation epoch is 
dominated by a hypothetical scalar field called the inflaton (p. In most inflation models, 
the inflaton is initially stuck in a high energy state and it is the slow release of energy that 
governs the inflationary period. The dynamics of the inflation is governed by some flat 
potential U{(j)), so that the energy density and pressure associated with this fleld are 



2.3.5 Inflation building 



1 + 3u; < 0. 



(2.47) 






(2.49) 



P4> = l4>^-U{4>)- 



(2.50) 
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If the field does not develop enough kinetic energy, = —p^ and = —1, meeting the 
condition in Eq. (12.470 for inflation to occur. If the value of the potential U{(j)) remains 
approximately unchanged during the field evolution, then the density and the pressure do 
not change as well during the inflationary period and Eq. (I2.48|) applies with = U{cj)). 
A precise mechanism for inflation is yet to be found, the major problems facing inflaton 
models being new fine-tuning troubles concerning the self-interaction of the inflaton field 
itself. Moreover, the inflaton field is, in the simplest models, a scalar field, and motivations 
from string theories justify the use of such fields, although fundamental scalar fields are yet 
to be discover experimentally. 



As mentioned in Sec. 12.3.41 one attractive feature of inflation is that scalar and tensor 
perturbations emerge during this epoch: these features later evolve into fluctuations in 
the primordial density and gravitational waves that might lead an imprint in the CMBR 
anisotropy and on the large scale structures [137^ [79l Wl\ 1166^ [18]. Each fluctuation is 
characterized by a power spectrum and a spectral index, respectively A^(A;), Ug for density 
perturbations and /^^-{k), nx for tensor perturbations. 

2.4.1 The scalar power spectrum 

The spectrum of the adiabatic density perturbations generated by inflation is specified by 
the power spectrum A|j(A;) that depends mildly on the comoving wavenumber k accordingly 
to a spectral index Ug and its tilt dus/dlnk as [111] 



The function A'^{k) describes the contribution to the total variance of the primordial density 
perturbation due to perturbations at a given scale per logarithmic interval in k |109] . 
The WMAP collaboration reports the combined measurement from WMAP7-|-BA0-|-SN 
of A^(/co) at the reference wavenumber k = ko = 0.002 Mpc~^ |110j . 



2.4 Fluctuations during inflation 




(2.51) 



(ko) = (2.430 ± 0.091) x 10 



(2.52) 



where the uncertainty refers to a 68% likelihood interval. The RHS of Eq. (I2.5ip is evaluated 
when a given comoving wavelength crosses outside the Hubble radius during inflation, and 
the LHS when the same wavelength re-enters the horizon. In Eq. ()2.5ip . we have used the 
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notation in Refs. \109\ 1110] for the density perturbations. Other authors use the symbol 
Pll{k) for our d3j^{k) and 1Z\ for our P'ji{k), and might differ by factors of 27r^. 

From a theoretical computation that describes the fluctuations in the inflaton field (p, 
the scalar power spectrum has the form 

Al{k)=(^y m'), (2.53) 



where <p is the time derivative of the inflaton field and describes the variance of 

the fluctuations in the inflaton scalar field, related to the spectrum of fiuctuations. The 
LHS of Eq. (|2.53p is computed at the time at which the largest density perturbations on 
observable scales are produced, corresponding Ne e-foldings before the end of inflation. For 
any massless and nearly-massless fields, the theory of quantum fiuctuations predicts [79] 

(|<5.^|'>quantum = ( ^) • (2-54) 

Since the Hubble rate H is approximately constant during inflation, see Sec. 12.3.51 we can 
substitute H in Eq. (j2.54p with its value when inflation ends Hj, with the index / standing 
for "Inflation": 
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(|<5</.|2)quantum = (^^J • (2.55) 

The Hubble expansion rate at the end of inflation Hj is bound by the WMAP measurements, 
as we will discuss in Sec. I2.4.4| and parametrizes the effectiveness of inflation. Using 
Eqs. ()2.53p and (j2.55p . the scalar power spectrum when quantum fluctuations in the inflaton 
fleld dominate is 



^niko)={^) . (2.56) 



2.4.2 The scalar spectral index 

The scalar spectral index Ug describes the mild dependence of the scalar power spectrum 
on the wavenumber k, as in Eq. (I2.5ip . We expand the spectral index around the reference 
scale ko as 

1 k 

ns{k)=ns + -Tln—, (2.57) 

2 ko 



where Ug = ns{ko), and the spectral tilt r is 

dns{k) 



dink/ko 

Using Eq. ()2.5ip . the scalar spectral index is [111] 



(2.58) 

k=ko 



9 , Al(k) 
olnk/ko A|j(fco) 



15 

2.4.3 The tensor power spectrum 

Fluctuations in the gravitational wave field are statistically described by a power spec- 
trum for tensor perturbations A|(A;). Writing this spectrum in a similar fashion to A^(A;o), 
we have 

^rW-^-^ = ^r(M (2.60) 

where the tensor spectral index ut is assumed to be independent of k, because current 
measurement cannot constraint its scale dependence. 

WMAP does not constraint A^(/co) directly, but rather the tensor-to-scalar ratio 



A^(fco) 
A2,(fco)' 



''-=f7?T- (2-61) 



which qualitatively measures the amplitude of gravitational waves per density fluctuations. 
The WMAP5+BA0+SN measurement constrains the tensor-to-scalar ratio as |110j 

r < 0.20 at 95% C.L. (2.62) 

2.4.4 Bounds on the Hubble rate at the end of inflation 

We combine the results from the WMAP-7 plus BAO and SN in Eqs. ([232]) and ([2:62|) 
to obtain an upper bound on the spectrum of primordial gravitational waves, 

A|(A;o) < 4.86 x lO^^". (2.63) 

Expressing A|(fco) in terms of if/, 

Off 2 

A?.(^-„) = (2^64) 

leads to an upper bound on i?/. 

Hi < 6.0 X 10^^ GeV. (2.65) 

A lower limit on Hi comes from requiring the Universe to be radiation-dominated at T ~ 
4 MeV, so that primordial nucleosynthesis can take place [98] [99[ [85] . Equating the highest 
temperature of the radiation 

Tmax ~ {TIhHiMpiYI\ (2.66) 
to the smallest allowed reheating temperature Tun = 4 MeV gives 

Hi > H{Trh) = 7.2 x lO'^^ GeV. (2.67) 
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2.4.5 Exiting inflation 



After a certain period in which the universe experienced an exponential growth, a 
transition towards the standard cosmology occurred. In the literature, this transition period 
goes under the general name of "reheating", the designation referring to the idea that it 
was at this time that most of the particles and radiation that formed the primordial soup 
were created. Since the detailed mechanism behind inflation is still obscure, the theory of 
the reheating process is yet to be specified in its entirety. Generally speaking, inflation 
ends when the potential of the inflaton field is no longer flat enough for the exponential 
solution to occur: the inflaton field starts rolling down and oscillates around a minimum 
of its potential U{(f)), with damped oscillations. The energy stored in the inflaton field is 
transferred through these damped oscillations into Standard Model particles and possibly 
other exotic particles, which make up the primordial soup [561 [T]. An alternative mechanism 
is the decay of the inflaton through a broad parametric resonance into intermediary particles, 
which then decay into Standard Model particles [55l 11761 1107] . In the standard cosmological 
scenario, the universe quickly becomes radiation-dominated, with Eq. (|2.37p describing the 
growth of the scale factor with time, a{t) ~ t^/^. 



Right after the end of inflation and the subsequent re-ionization, the early universe was 
filled with a hot plasma of Standard Model particles and possibly dark matter and other 
exotic particles and forms of energy. At such high temperature, most interaction rates were 
capable of keeping these constituents in thermal equilibrium; a specific particle i would go 
out of the thermal equilibrium when its annihilation rate into other particles Fj falls below 
the Hubble expansion rate H{T) at some temperature Tf^ec,i defined via 



From that moment on, the number density of such a relic is fixed to its value at temperature 



2.5.1 Statistical mechanics at thermal equilibrium 

Because of the thermal equilibrium existing among particles participating in this primor- 
dial soup, we can use statistical tools for describing the properties of each species. Here we 
review these methods, following the treatments in Refs. |108t I29j. In statistical mechanics. 



2.5 Thermal history of the universe 




(2.68) 
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the properties of a species i of mass rrii in a thermal bath at temperature T are described 
by a distribution function over the momentum p, 



Exp (^) ill" 



(2.69) 



where Ei = y + "T-f i is the total energy, is the chemical potential, and the minus sign 
(plus sign) specifies bosons (fermions). An additional number that characterizes the species 
is the number gi of degrees of freedom, describing the possible number of polarization of 
the species. In terms of these quantities, the number density rii and energy density pi are, 
respectively, 

ni = ^ / ^'^'^^(P)' (2.70) 

and 

Pi = ^ I d'pEihiv). (2.71) 
The pressure of the i fluid is obtained from the equation of state 

Pi = WiPi, (2.72) 

where Wi is specified by the type of fluid itself, see Sec. 12.2.2.11 

These quantities can be readily computed in the nonrelativistic T <^ mi and ultrarela- 
tivistic T ^ rui cases. In the former case we obtain, for both boson and fermion particles 

rii = 9i e-™'/^, for T « m^, (2.73) 



and 

Pi = mi Ui = gi mi 



r\ 3 



^ e"""'/^, forT^mj. (2.74) 



27r 

In the ultra-relativistic case, different results are obtained for the two statistics. 



■C(3) 



gi , for T ^ mi and Bose-Einstein statistics. 



Ui = { \ (2-75) 

I i 9i ^ ' ^ ^ ''^i Fermi-Dirac statistics, 

where C{z) is the Riemann zeta function of argument z, and 



?(T gi T^, for T mi and Bose-Einstein statistics. 



30 
7 

8 30 



fn 5* T ^ mi and Fermi-Dirac statistics. 



(2.76) 



Notice that the number density in the Fermi-Dirac statistics differs from the one in the 
Bose-Einstein statistics for a factor 3/4, and similarly the result for the energy density in 
the two cases differs by a factor 7/8. 
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2.5.2 Conservation of the number and entropy densities 

We now quote some important results concerning the evolution of a system of particles 
in thermal equilibrium, remanding to Ref. [29] for further details. The second Friedmann 
Eq. ()2.18p can be cast in the form 



d 
It 



0, (2.77) 



where the expression between square brackets is identified with the total entropy of the 
universe within a volume and p = ^.iPi- Introducing the entropy density 

s{T) = P±^, (2.78) 

we have that the entropy density scales with the Hubble volume a^. Using the fact that 
for a relativistic bath of particles p = p/3 and using the expressions for the various pi in 
Eq. ([2776]) . we obtain 

27r2 

s{T) = —g^s{T)T^ (2.79) 

with the entropy degrees of freedom 5*5 (T) at temperature T being defined in a similar way 
as g^{T) in Eq. ^SB, 

iSbosons j'Sfermions 

For T >1 MeV, we essentially have g*s{T) ~ g*{T). From the conservation of the entropy 
in a comoving volume a^{T) and using Eq. (|2.79p . we derive the relation between the scale 
factor and temperature as 

5* (T) (T) = constant , (2.81) 

valid only when no release of entropy occurs. 

If there is no release of entropy and g*s{T) is constant, the total number density is also 
conserved, 

dn 

- = 0. (2.82) 

The two conditions for the conservation of the number and the entropy densities can be 
cast in a single expression |108| . 

5 = 0. (2.83) 
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2.5.3 Application to the radiation-dominated universe 

We now turn our attention to a thermal bath of particles in cosmology. Considering an 
ensemble of relativistic particle species in thermal equilibrium, we evaluate the Hubble rate 
Hrad(T) from Eq. (I2.19|) using the expression for the energy density in Eq. (I2.76P as 

i 

where we defined the number of relativistic degrees of freedom at temperature T for an 
ensemble of species i that decouple at temperature T^eci in Eq- (|2.68p as 

E + l E *(%^)'- (2-85) 

iSbosons jGfcrmions 

Numerically, the function that gives the number relativistic degrees of freedom as a function 
of temperature for values of T around the QCD scale Aqcd can be approximated by a step 
function, 

'61.75, for r> Aqcd, 
5*(T) = <^ 10.75, for Aqcd > T > 4MeV, (2.86) 
^3.36, forT<4MeV. 

When all standard model particles can be treated as relativistic, we have (?*(T) = 106.75, 
while the number of relativistic degrees of freedom in the Minimal Supersymmetric Standard 
Model (MSSM) is g^{T) = 228.75. 
We simplify Eq. ()2.84p by writing 



i^,^ = \/*m|^|i = 1.06v^^. (2.87) 

This expression describes the Hubble rate when the universe is dominated by radiation. 

2.5.4 Matter-dominated universe 

As the universe cools down, relativistic particles lose momentum due to the redshift effect 
and eventually become nonrelativistic. For a stable nonrelativistic particle, the number 
density scales with n ~ ~ , with the latter relation coming from Eq. (|2.8ip . Knowing 
the matter energy density at the present time pM, the energy density at a time t is 

Pmatter = PM \~ J ' (2.88) 

from which the Hubble rate for a matter-dominated universe scales as H{T) ~ T^/^. 
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2.6 Content of the universe at the present time 

Even though we are still lacking a fundamental theory for inflation, the inflationary 
scenario has been embedded in the history of the cosmos due to its capability in solving all 
of the problems posed in Sec. 12.31 and make predictions on the amplitude of the seeds for 
inhomogeneities. This paradigm on the history of the universe explains the exceptionally 
flat universe we live in, Otot,o ~ 1, as evidenced from distinct measurements on the CMBR, 
the baryon acoustic oscillations and the redshift of supernovae. These same measurements 
also point out that the expansion of the universe is accelerating at the present time, a fact 
that is in sharp contrast with the naive expectation that the universe be matter-dominated 
at the present time. Clearly, the specific content of the universe is yet to be determined, 
although in the last decade we have been able to determine some general features of the 
different fluids that make it up. Here, we review the most important components of the 
present total energy density. 

2.6.1 Baryons 

In cosmology, the term "baryon" indicates the totality of the Standard Model species, 
and not only the color-neutral bound system made of three quarks in the acceptation of 
particle physics. The WMAP-7 data constrain the density of baryons today as 

l^b = 0.0458 ± 0.0016. (2.89) 

2.6.2 Dark matter 

It has long been known that baryons only account for a very small fraction of the 
present energy density. This conclusion was first obtained from considering the rotation 
curves of outer objects in galaxies, which reveal that the average galactic mass not only 
consists of dust and gases, but also and for its most part of a nonluminous halo of unknown 
composition, hence dubbed dark matter. 

The majority of the dark matter observed has to be in the form of CDM, which means 
that this exotic component has to be nonrelativistic at the time of galaxy formation. The 
fact that dark matter is a nonrelativistic field has been established with the first results 
on the CMBR anisotropy from COBE. In the CDM theory, small structures clump and 
grow hierarchically "from the bottom up", forming larger structures. This scenario is 
opposite to the Hot Dark Matter (HDM) paradigm, in which larger structure form earlier 
and subsequently fragment, following a "top - down" evolution. The predictions of the 
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CDM model are in general agreement with the observations, whereas the HDM paradigm 
disagrees with large-scale structure observations. The WMAP 7-years data, once combined 
with the BAO and SNe data, yield the value 

QcDM = 0.229 ± 0.015. (2.90) 

for the totality of CDM observed in the present universe. 

Although the CDM paradigm explains current data and the evolution of large-scale 
structures, there are some major discrepancies with observation of other features within 
galaxies and clusters of galaxies. In particular, CDM models predict that the density 
distribution of dark matter halos be much more peaked than what is inferred by the rotation 
curves of galaxies, the so-called cuspy halo problem. Moreover, CDM models predict a large 
amount of low angular momentum dust, in contrast with observations. 

2.6.3 Dark energy 

The dark matter and baryon components are not sufficient to explain the flatness 
observed, since the total abundance of nonrelativistic matter only accounts for about 26% 
of the total content of the Universe. Instead, what comes out of the measurements on the 
content of the universe reveals that a large fraction of the present energy density is due 
to the so-called dark energy, responsible for the current period of accelerated expansion. 
In some models, dark energy is identified with the constant A appearing in the Friedmann 
Eq. ()2.17p and which would lead to an accelerated expansion as discussed in Sec. 12.2.41 For 
anthropic reasons for this choice, see Ref. |191] . Another popular explanation for the dark 
energy introduces a new light scalar field whose equation of state, similarly to the inflaton 
field (/), resembles that of a cosmological constant. Mechanisms of this latter type include 
quintessence |148t 1193] . or the landscape of string theory (see Ref. |170j ). 

Evidence that the universe is experiencing a period of accelerated expansion comes 
from measuring the distances of type la supernovae (SNe). In fact, the lifespan of a SN 
is directly correlated with its luminosity, so that SNe can be used as standard candles to 
measure distances of neighboring stars to us via the "cosmic distance ladder" technique. 
Experiments reporting distances of SNe have shown that a nonzero cosmological term Vl\ 
better fits data than a vanishing Vt\. For a flat universe, the favored region of the parameter 
space has 

= 0.725 ±0.016. (2.91) 



CHAPTER 3 



ELEMENTS OF QUANTUM FIELD 
THEORY AND QCD 

3.1 Elements of group theory 

Group theory is a branch of mathematics devoted to the study of groups. A group G 
is a mathematical structure consisting of a set V together with an operation • that acts on 
a pair of elements in V . For this reason, the particular group structure is often explicitly 
indicated by writing G = G{V^ ■). In order for G to qualify as a group, the operation • must 
satisfy the following conditions: 

1. Closure: y x,y £ V, X • y £ V; 

2. Associativity: y x,y, z £ V, (x ■ y) ■ z = x • {y • z); 

3. Identity: 31 e = x ■ e = e ■ x = e\/x £ V ; 

4. Invertibility: \/x £ V3\ x^^ £ V = x ■ x~^ = x~^ ■ x = e. 

If further the commutativity relation holds, the group G{V, •) is called Abelian. 

Of central importance in the context of Quantum Field Theory (QFT) is the concept of 
the compact Lie group. The adjective compact refers to the compactness of the topology of 
the group, namely the topological space associated to the group G is compact; the notion 
of the Lie group refers to the smoothness of the differentiable manifold associated to G. 
Due to the importance of the subject, the theory of compact Lie groups is particularly 
well-developed. 

The properties of a compact Lie group assure that the group G contains elements x{a) 
that are arbitrarily close to the identity e for small values of the group parameter a. For 
continuous groups such as Lie groups, we can write 



x{a) = e + ia''T'' + 0{a^), 



(3.1) 
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where the index a runs over the dimension of the group and the T^'s are cahed the generators 
of the group. These generators satisfy the commutation relation 

[T'',T^]=if"-^''T^, (3.2) 

and the coefficients f"'^'^ are cahed the structure constants of the group. If the set of 
generators cannot be further divided into two subsets of mutuahy commuting generators, 
the generated group is called simple. For compact Lie groups that are simple, a complete 
classification of these groups is known [5], and it is indicated as the Cartan classification 
system. Given two A^-dimensional complex vectors u and v, a general transformation is, 
respectively, a 

1. unitary transformations if it preserves the inner product u*^Va- A unitary transfor- 
mation belongs to the so-called special unitary group in N dimensions, SU{N). The 
generators t"" of this group are N x N traceless Hermitian matrices, and there are 
iV^ — 1 of these generators. 

2. orthogonal transformations if it preserves the inner product Ua dab Vb- An orthogonal 
transformation belongs to the so-called orthogonal group in N dimensions, SO{N). 
The elements of this group are N x N orthogonal matrices of determinant one, and 
there are N{N — l)/2 generators. 

3. symplectic transformations if it preserves the inner product Ua EabVb, with Eab being 
the symplectic matrix. A symplectic transformation belongs to the so-called symplec- 
tic group in N dimensions {N even), Sp{N). There exist N{N + l)/2 generators. 

There also exist five exceptional Lie groups that are simple but do not belong to the above 
classification. These exceptional groups are of interest in unified theories, but they will not 
be treated in this thesis. 

Summing up, the dimension of the simple Lie groups, not including the exceptional 
groups, can be written as 



d{G) 



- 1 for SU{N), 
N{N-l)/2 for SO (N), (3.3) 
^N{N + l)/2 ioT Sp{N). 



In the following, we will specify the relations in Eq. (j3.2p for some groups of interest in 
particle physics, namely SU(2), S0(3), SU(3), and the homogeneous Lorentz group S0(l,3). 
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3.2 Lagrangian functions for scalar and spinor fields 

3.2.1 The Lorentz group 

As discussed in Sec. 13. H a group can be defined through the algebra of its generators. 
For the case of the Lorentz group S0(l,3), we have that the six generators for the angular 
momentum Jj and for the boost Ki, with i G {1, 2, 3}, satisfy the commutation relations: 

[Ji,Jj] = ieijkJk, (3.4) 

[Ki,Kj] = -ietjk Jk, (3.5) 

[Ji,Kj] = ieijkKk. (3.6) 

We see from Eq. (13. 4p that the three Ji form a closed subalgebra of the Lorentz group, more 
precisely the algebra that defines the SU(2) or S0(3) Lie groups. Representations for both 
these groups are three hermitian matrices of dimension 2n + 1, with n assuming integer 
values when referring to a representation of the S0(3) group, or n taking half-integer values 
when representing the SU(2) group. Nevertheless, each specific value of n in physics refers 
to a different particle that is associated with this group. For each different choice of n, the 
wave function for such a particle will transform differently when a Lorentz transformation 
is applied. We will focus only on the cases n = (boson) and n = 1/2 (Weyl spinor), with 
other values of n not being treated here. 

3.2.2 The Lorentz transformation 

In the following, we denote with A^^y the matrix describing the homogeneous linear 
transformation of a four- vector x^^, 

x^" ^ x'f" = A^'^x''. (3.7) 
Imposing the invariance of the pseudo-length of the four-vector, 

2_/j_ 1/ _ 2 2 2 2 Q\ 

X — X X ^ — ^fiu XX — Xq X-^ X2 X^ , (^O.o J 

we find the relation 

Aa^ Ap"" ri^^ = ri^p. (3.9) 

A Lorentz transformation, represented by the matrix A, belongs to the Lorentz group 
S0(l,3), and thus the product of two Lorentz matrices is a Lorentz matrix; the identity 6u 
belongs to the group and a specific transformation A has an inverse A~^. 
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Taking the determinant of Eq. (j3.9p , we see that 

DetA = ±l, (3.10) 

so that Lorentz transformations are called proper if Det A = +1 and improper if Det A = — 1. 
A further, independent subdivision can be made by noticing that Eq. (j3.9p implies that 

(A%)' = l+(A^o)', (3.11) 

so that a transformation is orthochronous when A'^q ^ !> and nonorthochronous when 
A% < -1. 

Transformations that are both proper and orthochronous are connected to the identity; 
in particular, we can write an infinitesimal transformation as 

A^ = 6i^ + e^, (3.12) 

where the matrix of infinitesimal quantities e^i, is forced by Eq. ()3.9p to be antisymmetric: 

(^^lu = -(-v^. (3.13) 

These are six parameters, corresponding to the three spatial rotations and the three boosts. 
Explicitly, a Lorentz transformation that rotates the system by an infinitesimal angle 59 
about the direction n and gives an infinitesimal boost along n with rapidity 6rj to the system 
has parameters 

(-ij = -€ijk fij 69, eio = hi 6r]. (3.14) 

3.2.3 Equation of motion for a scalar field 

We consider a generic complex scalar field ^ = <I>(x). The term "scalar" indicates that, 
when a Lorentz transformation A is applied to the coordinate four- vector x as in Eq. (j3.7p . 
the scalar field transforms as 

$(x) ^ $'(x) = $(A-^ x'). (3.15) 
The Lagrangian density for a classical scalar field $ moving in a potential U{^) is 

£ = a^$* -[/($), (3.16) 

and the equation of motion is 

^'^5.^-^ = 0. (3.17) 



In particular, when the potential is quadratic in the field, 
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U{^) = m^\^\'^, (3.18) 
Eq. (|3.17p is cast as the Klein-Gordon equation 

{df'd^-m'^) ^{x) = 0. (3.19) 
A similar equation holds for the complex conjugate field ^*{x). 

3.2.4 Equation of motion for a spinor field 

3.2.4.1 Representation of dimension two for the Lorentz 
transformation 



We focus again on Eqs. (j3.4p - (|3.6p that give the algebra of the Lorentz group in terms 
of the generators. We define two sets of linear combinations of these generators, 

N, = ^iJi-iKi), (3.20) 



and 



Nj = ^iJi + iK,), (3.21) 



so that the Lorentz algebra is rewritten as 

[N„Nj]=ieijkNk, (3.22) 

[iV/,iVj] = -^e,,,iV,^ (3.23) 

[Ni,N^^]=0. (3.24) 

This means that the algebra of each set Ni and A'^^ is separately closed, and the Lorentz 
algebra consists of two sets of SU(2) algebras, whose generators are hermitian conjugate. 
For this reason, a Lorentz transformation is specified by two indices ni, n2, and the 
representation has dimension (2ni + 1) x (2n2 + 1). The convention for Weyl spinors is that 
the choice (ni = 1/2, n2 = 0) leads to a left-handed transformation, while (ni = 1/2, 712 = 0) 
yields a right-handed transformation. These are two distinct representations of the Lorentz 
group with dimension two. 
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3.2.4.2 Weyl spinors 

Because of the two different representations of the Lorentz group with the same dimen- 
sion of two, two distinct type of spinors are possible, cahed the left-handed and right-handed 
Weyl spinors. These two representations are distinguished by the action of the boost 
transformations Ki, which in this representation are 

Kf = ±'-ai, (3.25) 

with the plus (minus) sign acting on the left (right) handed Weyl spinor. For both spinors, 
the generators of the angular momentum are given by 

Ji = ^ai. (3.26) 

The nonzero entries of the left-handed generator T^'^ are defined as 

n^ = e'^'jfc = ^e*^'Vfc, and T^' = K'^+ = '-a^. (3.27) 

We now discuss the action of a generic Lorentz transformation A on these spinors, 
considering first a left-handed Weyl spinor. In the following, we write the left-handed 
spinor as Xa, with the spinor index a raised and lowered by the spinor ni6tric tensor Cq,i^. 
Under the transformation in Eq. (j3.7p . the left-handed Weyl spinor transforms as 

Xa{x) ^ x'aix) = Lj{A) X/3(A-^x). (3.28) 

Here, L{A) is a matrix in the (1/2, 0) representation of the Lorentz group: for the infinites- 
imal Lorentz transformations in Eq. (I3.12p . it is 

L/(l + e) = 5/ + i(rr)/6,., (3.29) 

where we have indicated explicitly the spinor structure of the generator Tl, that is, the 
entries of the specific matrix with given fi and i^. 
The right-handed Weyl spinor is defined as 

xi = [Xa{x)]l (3.30) 

A dotted index a is used for the right-handed spinor, in order to distinguish it from 
the undotted index for the left-handed spinor. The metric tensor used for raising and 



28 



lowering dotted indices is e"^. Under a Lorentz transformation, the right-handed Weyl 
spinor transforms as 

Xa{x) ^ = rJ{A) X^(A^'^), (3.31) 

where R{A) is a matrix in the (0,1/2) representation of the Lorentz group. For an infinites- 
imal Lorentz transformation, we have 

Rj{l + e) = 5j+'-{T^'').^e,., (3.32) 

where T^'^ is a set of matrices, linked to the set of T^'^ matrices by 

m/ = -[(rr)/]*- (3-33) 

3.2.4.3 Equation of motion for the Weyl spinor 

Under specific conditions, it is possible to study the equation of motion for just one 
spinor, say the left-handed spinor Xa{x). The Lagrangian is |164j 

Cl = ixi {cy'-f^ d,xp + ^ (x"Xa + X^^X^ , (3.34) 
with the set of matrices 

a^ = (/2,fT), a^ = (/2,-<T). (3.35) 
From Eq. (j3.34p . we derive two equations of motion, one for Xa and one for xjj as 

i{(7napd^X^^ -^Xa = 0, (3.36) 

and 

iian"^ d^X^-mx^'' = 0. (3.37) 
These last two equations can be combined in the matrix form 



0. (3.38) 



3.2.5 Equation of motion for the Majorana spinor 

We define the four-spinor 

= ( ) ' (3-39) 

which is known in the literature as the Majorana spinor. The key element characterizing 
the Majorana spinor is that it is invariant under the charge conjugation operation, that is. 
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when we invert the electric charge of the particle, we obtain the same spinor. In order for 
this to be consistent, the electric charge of a Majorana field must be zero. To see this in 
more detail, we introduce two operations that act on generic four-spinors ^. 
A bar over a spinor ^ indicates the operation 

^ = -^^ (3.40) 



/3=(,". ?|. (3.41) 



*M= U-",,\t . (3^42) 



with the matrix 



For the Majorana field, 



We also introduce the charge conjugation operator C, whose matrix representation in four 
dimensions reads 

Given a generic four-spinor the charge-conjugated spinor is 

For a Majorana field, we obtain indeed = ^m- Examples of Majorana particles are the 
chargeless super symmetric partners of the Higgs, and photon. A question mark is still 
posed on the nature of the neutrino, which may also be described by a Majorana field. 
Defining the Dirac matrices in the chiral representation 

;"). (-) 

we have the Dirac equation for a Majorana field 

- m) = 0. (3.46) 
This equation of motion can be derived from the Lagrangian 

Ui = ^ C 7'^9^*M -^m^ljC ^M, (3.47) 
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3.2.6 Equation of motion for the Dirac spinor 

In order to describe spin one-half particles that possess charge, such as electrons, we 
consider a four-spinor in which two different and unrelated Weyl spinors Xa, C^" appear. 
We construct the Dirac four-spinor as 

= ( p ) . (3.48) 

Under charge conjugation, the Dirac field does not transform into itself; in fact we have 

*g = ( ) / ^D, (3.49) 

so that nonzero values of the electric charge are allowed. The Dirac field satisfies the 
equation 

which is put in the compact form 

{i-i^d^, -m)^D = 0. (3.51) 

The fact that the Dirac field allows us to describe particles with a nonzero charge can be 
seen from the Dirac Lagrangian, 

CD = i^Dl^d^^D-m^D^D. (3.52) 
which is invariant under a U(l) transformation of parameter w, 

^D^e^'^'^D, ^D^e-'^^^D- (3.53) 

3.3 Application of group theory to QFT 

We now turn our attention to gauge theories in Quantum Field Theory (QFT), focusing 
on Dirac fields. For this reason, we drop the index D and simply denote with ^ a collection 
of one or more Dirac fields. The Dirac Lagrangian in Eq. (j3.52p . describing a free field 
theory for a massless complex- valued collection of Dirac fields, reads 

C = i^{x)Y d^.-^ix). (3.54) 

We assume that the collection of fields ^(x) is invariant under a global symmetry that 
generalizes Eq. (|3.53p . 

*(x) ^ C/(a'')^'(x), (3.55) 




0, 



(3.50) 



31 

where the transformation of the field can be written in terms of a set of constant parameter 

a"- as 

U{a) = exp^ia^t") . (3.56) 

For example, if ^{x) represents a single Dirac field, then the Lagrangian in Eq. (|3.54p is 
invariant under the U{1) transformation in Eq. (|3.53p . while if ^(x) represents a doublet of 
Dirac fields then the Lagrangian in Eq. (j3.54p is invariant under a SU{2) transformation, 

^'(x) exp (^ia ^'(x), (3.57) 

where a = {a^,a'^,a^} are the usual 2x2 cr-matrices. 

In QFT, a general recipe is adopted in order to include the interaction of the Dirac field 
^(x) with a gauge field A^. We first promote the symmetry from global to local by making 
the parameters a depending on the space-time coordinate x, 

a^a{x). (3.58) 

The transformation of the Dirac field ^{x) under this local symmetry reads 

^(x) ^ Uia^ix)) ^{x) = exp (ia^ix) t") ^-(3;). (3.59) 

This local transformation resembles that in Eq. (13.550 . the main difference being that the 
parameters a"'{x) now depend on the specific space-time point x at which the transfor- 
mation takes place. This justifies the use of the adjective "local" to describe such kinds 
of transformation. In order for this transformation to be imposed on the Lagrangian in 
Eq. (j3.54p . we define a covariant derivative 

D^ = d^-igAie, (3.60) 

where g is the coupling of the Dirac field to the gauge field. The Lagrangian in Eq. (j3.54p 
then reads 

£ = ^!{x)iY ^'(x) = ^{x)i-/>' {df, - ig A"^ ^'(x). (3.61) 

We see that the introduction of the covariant derivative introduces in the Lagrangian above 
the interaction term 

C,nt = g^{x)i^Aie^>{x). (3.62) 

The requirement that the Lagrangian in Eq. (j3.6ip be invariant with respect to the local 
transformation in Eq. (|3.59p imposes the gauge field A'^ to transform as 

Al^ Al + -d^a''e + A^^a^, (3.63) 
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and f"'^'^ are the structure constants defined in Eq. ()3.2p for the group G with generators 
t". 

Finally, we can write the complete Lagrangian for a Dirac multiplet ^{x) of mass matrix 
M and belonging to an irreducible representation of a gauge group G, known as the Yang- 
Mills Lagrangian, as 

Cym = ^{x) {i r - M) ^{x) - \ rr F;,. (3.64) 

In the Yang-Mills Lagrangian, the index a is intended summed over the generators of the 
group G, while we have introduced the field strength 

f;, = d,Al-d,Al + g 4 Al. (3.65) 

In Eq. ()3.64p . we have included the mass term M which is manifestly invariant under 
the symmetry in Eq. (13.59p . and the term 

£G.nse = -\FrF^,, (3.66) 

that describes the self-interaction and the dynamics of the gauge field. 

3.4 Quantum Chromodynamics 

We apply the tools developed in the last sections to describe the theory of Quantum 
Chromodynamics (QCD), and provide a theory to describes the mechanism of the strong 
interactions. The QCD theory is nonabelian since gluons carry color charge causing them to 
interact with each other in a more complicated structure: in fact, the underlying gauge group 
for QCD is the nonabelian group G = SU{3). As we will see, a yet-to-be-solved problem 
arises in QCD, the so-called strong Charge-Parity (CP) problem. In fact, as in the theory 
of weak interactions, the CP symmetry is expected to be violated in strong interactions, 
whereas experiments show that the CP symmetry is preserved by strong interactions to a 
high degree of precision. The lack of a CP-violating term then requires a fine tuning of 
QCD parameters in its mathematical description. 

3.4.1 The QCD Lagrangian 

We denote the QCD gauge field (the gluon tensor field) living in the adjoint represen- 
tation of the group as GJ^^, where a G {1, 8}. The Lagrangian obtained using the SU{3) 
gauge group with Nf quark (the QCD Lagrangian) is [143j 

Cqcd = qin''D^-M)q-\ G^'^^G^, (3.67) 
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where q is the collection of quark spinor fields, M is the quark mass matrix |138j and G^j^ 
is the gluon tensor field, which can be written in terms of the QCD field strength and 

fabc 
'5(7(3) 



of the structure constant /cr?co\ of the SU (3) group as 



G^, = d^Al-d,A';, + g ffjj'^,^ A"^ At. (3.68) 

The covariant derivative appearing in Eq. (j3.67p is defined in terms of the Gell-Mann 
matrices ^5^(3) by 

D^ = d^-igA''esui3y (3-69) 

In the limit of vanishing quark masses M — )• 0, the QCD Lagrangian is invariant under the 
global transformation UiiNf) 05 Uii{Nf). However, such symmetry implies that hadrons 
come in doublets |144j . which is not seen experimentally. In fact, it turns out that the axial 
part of the chiral symmetry is spontaneously broken into the subgroup SU{Nf)A U{1)a- 
Because of the symmetry breaking, nearly massless Goldstone bosons arise: these bosons 
are experimentally identified with the pseudo-scalar octet formed by the u, d, and s quarks 
and explain the SUA{Nf) breaking. However, there exists no candidate in the particle 
spectrum that would play the role of the Goldstone boson associated with the spontaneous 
breakdown of U{1)a- In fact, such a particle should be a light pseudo-scalar boson of 
mass iria < VSm.,^: the latter requirement excludes the rj' meson which would otherwise 
be a suitable candidate having the proper quantum numbers. The absence of this light 
pseudo-scalar particle is known as the U{1)a problem |189j . 

3.4.2 Solving the U{1)a problem 
As proposed by t'Hooft [17211173] . the U{1)a problem might be solved by adding a term 
to the QCD Lagrangian in Eq. (j3.67p that explicitly breaks the U{1)a symmetry, 

^QCD — ^ ^QCD + (3.70) 

with 

^o = ^^G'^'"'Gt.u- (3.71) 
Here, we have defined the dual of the field strength tensor 

^fiu ~ 2 ^M'^'^p ^ ' (3.72) 

with e^i,o-p the Levi-Civita completely antisymmetric tensor having eoooo = 1, 5 is a coupling 
constant, and 6 is an observable |9H I33j that might take values form zero to 27r. The origin 
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of the variable 9 in the QCD Lagrangian can be explained by the following argument. 
The classical gluon field equations admit an anti-instanton solution, which satisfies the 
antiduality condition 

Gfj,,/ = G^j/, (3.73) 

and which has an integer index 

^ / d^xG^l^G^'^^ (3.74) 



n 



The angular variable 9 parametrizes the linear combination of different |n)-vacua in the the- 
ory corresponding to different values of the integer n, with the particular linear combination 
known as the 0-vacuum, 

\9)= e-^"^'|n). (3.75) 



Here, every value of n characterizes the winding number of the U{1)a symmetry in Eq. (j3.74p . 
Although the term in Eq. (|3.7ip is a surface term, it is possible to show that if the U{1)a 
problem is solved and none of the quark is massless, a nonzero value of 9 implies that the 
CP symmetry is broken |157j . Thus, QCD physics depends on the value of 9, or better, 
when electroweak interactions are included, on the combination 

^ = ^ + ewcak, (3.76) 

with the extra term arising in the QCD Lagrangian of the form 

A = O^cak G'^"" G^^. (3.77) 

The parameter ^weak is expressed in terms of the quark mass matrix M as 

^weak = arg(DetM). (3.78) 

Eventually, the QCD Lagrangian including the 'tHooft and electroweak interaction terms 
reads 

^QCD+e = -Cqcd + C§ = £qcd + ^l^f {3.79) 
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3.4.3 The strong CP problem 

As we discussed in Sec. I3.4.2t it is expected that the CP symmetry is violated by strong 
interactions due to the CP-breaking term Cg. However, experimentahy no violation of the 
CP symmetry is observed, and the strong interactions preserve this symmetry to a high 
degree of accuracy. This problem in merging theoretical motivations and experimental 
observations is known in the literature as the strong CP problem, which can be restated by 
asking why strong interactions do not violate the CP symmetry when CP violation is not 
forbidden in the theory. 

The most stringent constrain on the violation of the CP symmetry by the strong 
interaction comes from the measurement of the electric dipole moment of the neutron, 
defined as 

d7v~^^^, (3.80) 

where e is the absolute value of the electron charge, niN is the neutron mass and is 
given in terms of the masses of the up and down quarks m^, nid as 

rriud = ■ • (3.81) 

ruu + rud 



The latest experimental bound on the neutron electric dipole moment is [15] 

I^Tvl < 2.9 X W^'^^ecm at 90% C.L., (3.82) 

while we take the value of cIn- from the review by Kim and Carosi [104] . 

dAT = 4.5 X lO^i^^ecm. (3.83) 

Eqs. (|3.82|) and (|3.83p imply the constraint 

|0| < 0.7 X 10"". (3.84) 

From the definition of 6 in Eq. ()3.76p , we see that this quantity is the sum of the two terms 
9 and ^wcaki whose physical origin is completely unrelated; it is thus a mystery why two 
quantities which are naturally of order one have to cancel out in a sum with such a high 
degree of accuracy as in Eq. (j3.84p . 
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3.4.4 Possible solutions to the strong CP problem 

As of today, three main solutions to the strong CP problem have been proposed. These 
are, respectively: calculable 6, massless up quark, the axion particle. As the axion is the 
most favorable solution and the subject of this thesis, we will briefly review here only the 
first two cases and concentrate on the axion throughout the subsequent chapters. 

3.4.4.1 Calculable 6 

The underlying idea of the proposed calculable 9 is to impose CP invariance in the QCD 
Lagrangian, setting 9 = 0. Taking into account that the only source of CP violation comes 
from weak interactions, as explicit in phenomena like the neutral K- meson oscillations and 
the B-meson decay, the term ^wcak is still nonzero and can be calculable in an underlying 
theory, with a small value constrained by observations as in Eq. (|3.84p . An example of 
such a model is the Barr-Nelson CP violating model, in which heavy singlet quarks are 
introduced and whose vacuum expectation value lies much above the electroweak scale. 
These quarks mix with the ordinary light quarks after symmetry breaking, so that at low 
energy they can be integrated out and the light quarks acquire the required amount of CP 
symmetry through the Cabibbo-Kobayashi-Maskawa mixing. In this theory, the condition 
arg(Det M) = is imposed at tree level. A nonzero value of ^wcak arises through higher-loop 
corrections, ^weak ^ ctloopi with the fine structure constant aioop bound by the stringent 
bound in Eq. i^M\i . 

3.4.4.2 Massless up quark 

If the lightest quark (the up quark) is massless, it is possible to eliminate ^ by a rotation 
of the quark fields q, 

q ^ Exp (^^ ^7^^ q, (3.85) 

so that 9 is no longer an observable. In fact, the QCD Lagrangian in Eq. ()3.79p transforms 
under the rotation above as 

^QCD+e ^QCD, (3.86) 

making the term Cg disappear. The massless up quark possibility has been ruled out by 
measurements, as quoted in Eq. (|4.10p below that states m^/m^ = = 0.568 it 0.042. 
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THE AXION 

4.1 The axion as a solution to the strong 
CP problem 

The third and the most compeUing solution to the strong CP problem discussed in 
Sees. 13.431 and [3.4.4l introduce an additional chiral symmetry that promotes the 9 parameter 
to a dynamical field: the dynamics of this new field naturally relaxes its expectation values 
towards arbitrarily small values, thus eliminating the CP problem. This mechanism is 
achieved by introducing a new global, chiral symmetry J7(1)pq which is spontaneously 
broken at an unknown energy scale fa, known as the PQ energy scale. The axion is the 
Goldstone boson resulted from such symmetry breaking. 

The Lagrangian term resulting after the spontaneous symmetry breaking of the [/(1)pq 
symmetry reads 

= -\ [d^a) {d,a) + ^^G'^'^^GI,,, (4.1) 

where a = a{x) is the axion field and = Xf is a new parameter called the PQ color 
anomaly, defined as the sum of the PQ charges Xf over the fermions in the theory. In 
Eq. ()4.ip . the first term represents the kinetic term, while the second term describes the 
interaction of the axion with the gluon field. 
When considering the Lagrangian 

^QCD+e+a = ^QCD + + Ca, (4.2) 

an effective potential V^fr(a) for the axion appears, whose minimum is reached when 

{a{x) + ^e)=0. (4.3) 

The expectation value of the axion field in Ca at its minimum cancels out the 9 term in Cg, 
thus eliminating the CP-violating term. Since we can shift the axion field by a constant 
amount without changing the physics, we define the (axion) misalignment angle 
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Summing up, with the PQ solution to the strong CP problem, the free parameter 9 has 
been traded for a dynamical field a{x) that evolves to its CP-conserving minimum 6 = 
through the spontaneous breaking of a U{1) symmetry. Essentially, 6 can be seen as the 
phase of a new complex scalar field 

^a = \^a\e'', (4.5) 

that takes a nonzero expectation value at the energy scale fa/N. The axion potential thus 
has the usual Mexican hat form 

^pqW = ^ fl^ap-(^)'') +m2(r)/f (1-COS0), (4.6) 



where A is a coupling constant. Here, the second term in Eq. ()4.6p comes from nonpertur- 
bative QCD effects associated with instantons [77], that break the C/(1)pq symmetry down 
to a Z{N) discrete subgroup |159j . As we will see later in Sec. 14.2. 2^ the presence of this 
second term justifies the fact that the axion mass ma{T) depends on the temperature of 
the plasma considered |181j . After the spontaneous breaking of the PQ symmetry, we are 
thus left with the residual potential (but see Refs. [62l [6ll [63] for different potentials) 

V{e) = ml{T)f!{l-cose). (4.7) 

This potential is approximately quadratic for small values of 9, but quickly shows deviations 
from a pure quadratic potential for ~ 1. 

4.2 On the mass of the axion 
4.2.1 Axion mass at zero temperature 

We can expand the term Cg + Ca around the minimum in Eq. (j4.3p to obtain the value 
of the axion mass as 



« " 327r2 /, da \ 



(4.^ 

{a{x))=-fa e/N 



The axion mass has been computed in Ref. |190| including the axion mixing with the other 
neutral pions through the axion-gluon-gluon coupling, and is written in terms of the quark 
ratios Zd = niu/md, Zg = mu/'ms and of the vr*^ mass and decay constant m^o and /^o as 

"''^ = ^-^jf- ^ [i + zd){i + zd'^y ^^-^^ 
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Using the experimental values |116| 

^rf = 0.568 ± 0.042, = 0.029 ± 0.003, (4.10) 

and the values [9] 

m^o = (134.9766 ± 0.0006) MeV, f^o = (130 ± 5) MeV, (4.11) 

V2 

we obtain the numerical value 

ma = (5.91 ± 0.53)/ieV , V^ . (4.12) 

Ja,12/JV 

Here and in the following, a number y indexing some quantity with units of energy indicates 
that such quantity has been divided by lO'^ GeV: for example, in the equation above, it is 

/a,12 = /a/10l2 Q^y. 

The axion mass in Eq. (|4.12p is exact in the context of the simplest hadronic axion model 
|103t I157j (see Sec. 14.4. 2. 1|) . where Standard Model fermions do not possess a PQ charge, 
and axions only interact with an exotic heavy quark. Above the QCD phase transition 
nonperturbative effect arise, giving the axion a temperature-dependent mass as we now 
discuss. 

4.2.2 Finite temperature effects 

In a quark-gluon plasma of finite temperature T, nonperturbative effects modify the 
properties of the components of the plasma itself. Nonperturbative effects in QCD have been 
studied via lattice simulations, see ex. Ref. [38], or phenomenology, see ex. Refs. [50 } 1158] . 
For axions, the mass acquires a temperature dependence at sufficiently high T due to the 
instanton effects: In the high temperature limit, T ^ Aqcdi the axion mass is given by an 
integral over noninteracting instantons of all sizes [77] as (see also Refs. [178^ [60]) 

33 - 2iV A 



with 



mliT) « c (^^) (^) ' Det (M) X(T), (4.13) 

I(T)= dxx^+^ {\n- e--^(^\ (4.14) 

Jo V ^QCDxJ 



12 a , / vr^ ^ 
In iH x' 



(4.15) 

(l + 7(7rx)-3/2)" > ■ o , - 
and the numerical values c = 0.130078, x = 1-33876, a = 0.01289764, 7 = 0.15858 and 
a = (153 — 19A'^y)/(33 — 2Nf). The parameter Nf gives the number of quarks that are 



40 

relativistic at temperature T, that is, whose mass is < T. It is worth stressing here that, 
using different computation techniques, various authors have obtained different results for 
ma{T). For later convenience, we parametrize the general expression for the temperature- 
dependent mass with a general broken-law function, 

m fa In r>AQCD, 

7na{T) = TTla < \ ' ' (4.16) 

[ l> ^ ^ Aqcd, 

with the parameter c = 7/2 + Nf/6 from Eq. ()4.13p . Before showing our computation for 
the parameters a and d, we will briefly review the literature on the subject for the nontrivial 
case T > Aqcd- 

Earlier work on the temperature-dependent axion mass are in Refs. |1451 [2t [53], where 
the authors use the method by Gross, Pisarski, and Yaffe in Ref. [77] in the dilute gas 
approximation. Preskill et al. |145| obtain, for the case A'^^ = 3, the expression 



1/2 / , s 4 



m.(T) = 2 X 10-2 ^ ( ^ ) 

Ja \ TTl J 



9 In 



ttT \ 
A-QCD/ 



3 



(4.17) 



which we put in the form of the first line in Eq. (j4.16p with Aqcd = 200 MeV and Det (M) = 
1000 MeV^ around T ~ 4Aqcd [US] as 

ma{T) = aprcskill "T-a ^ ^^" ^ , (4.18) 

with Opreskill = 0.035. 

Sikivie [161] quotes the general result from Refs. [145^ O [53] as 

».(r)^4xlO-»ev(12^)(^|X)\ (4..) 

We notice that this expression can be written as 

ma{T) = asikivic nia ( ^^^ ] , (4.20) 



T 

with asikivie = 0.417/ A^ 2 and A0.2 = AQCD/200MeV. 

Turner |178| uses the tools developed in Ref. ^7] to compute ma{T) in the case T S> 
^QCD and gives the general expression 



ma{T) = rua a A'^q 2 



^QCD \ ' 



1 - In- 



A, 



QCD 
T 



-I d 



(4.21) 



T J 

where rua = rna{T = 0) is given in Eq. (I4.12p . Table. 14.11 shows the dependence of the 
parameters a, b, c, and d on the number Nf of relativistic quarks at the temperature Ti at 
which the axion field starts oscillating, see Eq. (j5.15p 
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Table 4.1. Values of the parameters a, b, c, and d appearing in Eq. (j4.2ip . as a function 
of the number of relativistic degrees of freedom Nf. 



Nf 


a 


b 


c 


d 


1 


0.277 


3/2 


3.67 


0.84 


2 


0.0349 


1 


3.83 


1.02 


3 


0.0256 


1/2 


4.0 


1.22 


4 


0.0421 





4.17 


1.46 


5 


0.118 


-1/2 


4.33 


1.74 


6 


0.974 


-1 


4.5 


2.07 



The exact value of Nf is not well determined, because the value of Ti is of the order of 
the current mass of the charm quark, nic ~ 1 GeV, so it is possible for it to be Nj = 3 or 
Nf = 4. Because of this uncertainty, Ref. |178j considers the function 

/A N 3.7±0.1 

m.(r) = 7.7 X 10-^±°-5 j , (4.22) 

where the uncertainties interpolate between the various values of the parameters in Table HTT] 
for different plausible values of Nf. 

DeGrand et al. [49], using the formula in Ref. [77], quote an expression for the axion 
mass valid for T » AqcD) as 



In- 



Aqcd 



(4.23) 



This is basically the same expression found by Preskill in Eq. (j4.18p . 

Bae, Huh, and Kim |14j improved the results by Turner by considering updated values 
of current quark masses and the effects of the QCD phase transition. They parametrized 
the axion mass with 

'"•'^) = °'°-' /j(T/GeV)- - (4.24) 



For the parameters ainst and z^, the authors find a Aqcd dependence like in Table 14^2] 

Wantz and Shellard [1871 1188] . within the interacting instanton liquid model [153] . 
obtained a numerical expression for the temperature dependence of ma{T) valid at all 
temperatures; in particular, they were able to follow the ma{T) function around the QCD 
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Table 4.2. Values of the parameters ainst and v in Eq. ()4.24p for different values of Aqcd- 
From Ref. [H] 



Aqcd 




u 


320 


0.9967 


6.967 


380 


3.964 


6.878 


440 


12.74 


6.789 



scale AqcD) where analytic methods break. The analytic fit in Refs. |1871 118S] reports, for 
T > Aqcd, 

with aa = 1.68 x 10^7 and i/ = 6.68. 

Here, we consider the case Nf = 3, because, as we will show later, the axion field 
starts oscillating at a temperature Ti that lies below the GeV; with only three quarks 
being relativistic at Ti we have Nf = 3 and thus, from Eq. ()4.13p and from Table 14.11 the 
axion mass decreases with T as to the leading order. For this reason, we set c = 4 
in Eq. (j4.16p . To compute the parameters a and d, we proceed as follows. Based on the 
latest available data for the masses of the quarks u, d, s, which are the relativistic quarks 
at temperature Ti, we write the determinant of the quark mass matrix as 



Det (M) = m„ md = — - = (1.64 ifc 0.29) x 10'^ MeVM ^" J , (4.26) 

ZrfZs V3MeV/ 

where and Zg are given in Eq. (j4.10p and m„ lies in the range (1.7 - 3.3) MeV |138j : here 

we take the value = 3 MeV. To evaluate ma(T) we use Eq. (|4.13p with Nf = 3. We are 

particularly interested in evaluating the function I{T), whose plot is shown in Fig. 14.11 

Since the integral in the function X(T) in Eq. (|4.13p depends on the temperature, we fit 

the numerical solution of such function and obtain 

/ T \ ^-^^ 

I{T) « 5.84 X 10"*^ 1 + In . (4.27) 

V Aqcd / 

Using this fit for X(T) with the values of the parameters as above and Nf = 3, we find for 
the temperature-dependent axion mass in Eq. ()4.13p the expression 

ma(r) = (0.017 ± 0.003) A^/^' (^)' (l- In ^)''', forT>AQCD, (4.28) 

while ma{T) = nia in Eq. (j4.12p for T < Aqcj^. The uncertainty in the prefactor a = 
0.017± 0.003 is due to the uncertainties in z^ and Zg only, since we neglected the uncertainty 
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Figure 4.1. The function I{T) in Eq. ()4.14p as a function of temperature T for Aqcd = 200 
MeV. 

in the mass of the up quark. Comparing our Eq. (j4.28p above with the result by Turner in 
Table [4T| we see that we confirm the exponent d = 1.2, but we differ from Turner's result 
because of the prefactor a = 0.017, which is about 50% lower than a = 0.0256 in Table HTl 
We attribute this discrepancy to the different value of the up quark mass used in Ref . [178j . 
Notice that the value of our prefactor a is compatible with the value a = 0.018 used in 



For T of the order of Aqcd, we neglect the logarithmic term; in the following, we will 
always assume this is the case, adopting the following expression for the axion mass 



Thanks to the Lagrangian term Ca, axions couple at tree level to gluons and photons. 
However, depending on the specific particle model, axions may couple differently to the 
Standard Model fermions and other gauge bosons, and the coupling to photons might be 
suppressed. 



Refs. (2il[60l[88l[T85]. 




(4.29) 



with a = 0.017 ± 0.003. 



4.3 Coupling of standard model particles 
with the axion 
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4.3.1 Coupling of axions to gluons 

In all axions models, axions couple to gluons through the second term in the Lagrangian 
in Eq. gl]), 

£a, = i^J^G'-G;„, (4.30) 

where we included the expression for the fine structure constant of strong interactions 
Us = /4iT. Thanks to this coupling to gluons, the axion mixes with pious and acquires 
the zero-temperature axion mass in Eq. (j4.9p . 



4.3.2 Coupling of axions to photons 

The axion-photon coupling is described by the term 



077 



gayaE- B. 



where (7^^ is the axion-photon coupling, given by 



9a-y 



a 



271 fa 



E_ 2{4 + z + Zs) 
N 



a 



27r fa 



E 

N 



1.92 ±0.08 



(4.31) 



(4.32) 



3(1 + z - Zs 

where a ~ 1/137 is the fine structure constant. The ratio E/N can be written in terms of 
the PQ charges Xf, the electric charges Qj and the number of colors in a multiplet Df as 

Depending on the model, ga-y can be enhanced when the number of quarks Df is large, or 
suppressed when E/N ~ 2. In the following, we indicate the Feynman diagram associated 
to the Lagrangian term in Eq. (14.32P with the graph shown in Fig. 14.21 

Models in which there exists an axion-photon coupling at tree level predict that the 
Primakoff effect for axions might have an important role in a number of astrophysical 
processes like the cooling of stars, and enable us to attempt the detection of axions through 
a "shining through a wall" experiment. 




Figure 4.2. The axion-photon vertex in an axion theory in which ga-y 7^ 0. 
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4.3.3 Coupling of axions to fermions 

In the Lagrangian term describing the interaction of axions with the fermionic field 
^j, the axion field always appears in the derivative term d^a{x), so that the invariance 
a — >• a + const, proper to Nambu-Goldstone bosons is imposed: 

C 

^aff = i7f^fl''To^fdf,a. (4.34) 

In the expression above, is the fermionic field of mass mj and Cj is a model-dependent 
parameter that gives the PQ charge of the fermion /. Writing the interaction Lagrangian 
as 

^a,ff = -i^^^fl5^fa, (4.35) 
it appears that the combination 

Saf = (4.36) 

la 

plays the role of a Yukawa interaction, with Oaf = g^^/A-K being analogous to the fine 
structure constant in this model. Fig. 14.31 shows the Feynman diagram associated with 
Caff for the fermion /. 

4.4 Axion models 
4.4.1 The "visible" axion model 

The first viable axion model is known as the PQWW model from the initials of the 
authors Peccei and Quinn [141j . Weinberg |190| . and Wilczek |194j . or the "visible" axion 
because this theory would lead to sizable effects in various physical phenomena due to the 
strong axion interaction with light and light mesons. In the PQWW model, the PQ energy 
scale fa is related to the weak interaction scale A^cak ~ 250GeV. A pair of Higgs fields <I>i 



Figure 4.3. The Feynman diagram showing the axion-fermion vertex for a fermion field /, 
in an axion theory in which gaf 7^ 0. 
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and ^2 is introduced, the first giving mass to the up quark and the second to the down 
quark, with each Higgs field having a vacuum expectation value (^'1,2) = ^1,2/^^7 with 
^y Af + A2 = A^cak- From these premise it is possible to bound the PQ scale as fa < 42 
GeV and so the axion mass ttIq > 200 keV. 

The PQWW axion model has long been abandoned because its very premise fa ~ A^cak 
has been ruled out by various astrophysical considerations (see Sec. I4.5p and because this 
model predicts a lifetime of various mesons like or J/^' that would be too short compared 
to observations. 

4.4.2 Models for the "invisible" axion 

After experiments showed that the PQ energy scale cannot be related to the weak 
interaction scale, authors have either considered fa as a free parameter or they have 
considered more complicated models in which fa relates to the Grand Unification Theory 
(GUT) scale Aqut ~ 10^^ GeV or to the Planck scale Apianck ~ 10^^ GeV. In the first case, 
the axion is studied in a model-free theory and its generic properties are constrained by 
astrophysical considerations, direct and indirect searches, and accelerator experiments. In 
the second case, the axion is embedded in a more fundamental theory such as a GUT or a 
string theory. 

All these models share a PQ scale much higher than the A^cak energy scale, so that 
the coupling of axions to Standard Model particles is considerably weaker, thus making the 
experimental detection of the axion even more challenging. 

4.4.2.1 The KSVZ model 

The first invisible axion model has been proposed by Kim |103| and by Shifman, Vain- 
shtein, and Zakharov |157j . and it is thus known as the KSVZ model. Since, in this model, 
the axion does not couple with SM fermions at tree level, being the Cf = in Eq. (j4.34p for 
SM fermions, in the literature this model is also referred to as the hadronic axion model. 
The KSVZ axion couples at tree level to gluons and to an exotic heavy quark Q; interaction 
with the other SM particles occur at one loop through the intermediation of gluons and of 
the new Q quark. Depending on the charge of Q, the ratio E/N in Eq. (I4.33P may take a 
value in the range [0, 6]. In particular, the KSVZ model allows the value E/N = 2 at which 
the axion-photon coupling gaj in Eq. (j4.32p is suppressed. 



47 

4.4.2.2 The DFSZ model 

In the DFSZ model, introduced by Dine, Fischler, Srednicki [5l], and by Zhitnitskii 
|202j . the axion couples at tree level to SM photons and charged leptons, besides nucleons. 
One of the main purposes for the introduction of this model, and an advantage for it, is 
that it can be easily embedded in GUT models. Any GUT model predicts E/N = 8/3, and 
so a precise value of the axion-photon coupling, 

10"^ 

= (8.67 ±0.93) X -— . (4.37) 

la 

The model also predicts the value of the axion-electron coupling Ce appearing in Eq. (|4.34p 
with / = e, as 

Ce = (4.38) 

with cot the ratio of two Higgs vacuum expectation values of this model. The coupling 
of axions to nucleons are also calculable in the model, and they are related by generalized 
Goldberger-Treiman relations to nucleon axial- vector current matrix elements. Finally, the 
axion coupling to the up and down quarks are 

C. = ^, and Q = ^^^. (4.39) 

4.4.3 Lifetime of the axion 

Similarly to the neutral pion, the coupling of the axion to two photons arises through the 
electromagnetic anomaly of the PQ symmetry, and the decay of an axion into two photons 
is allowed with lifetime 

647r , , 94 (eV/iTia)^ 
= ^r^'^'^'' ^ I(i;/A'-1.92±0 08)/0.75P - ("») 

In the last expression, we used the values of fa and in Eqs. ()4.12p and (j4.32p . In analogy 
with the theory of pion decay, the lifetime of the axion depends on the fifth power of the 
axion mass. When considering light axions of mass rua ^ 1 eV, the lifetime ra-5.27 ~ 10^^ s 
implies that these particles are cosmologically stable. 

4.5 Astrophysical bounds on axions 

If axions existed, they would be produced in the hot plasma that constitutes stars and 
other astrophysical objects. The presence of axions in this dense environment would open 
up additional channels for the occurrence of well-studied astrophysics processes, and thus 
it would alter star evolution. 
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4.5.1 Constraints from the cooling time of white dwarfs 

When helium-burning stars reach their latest stages of helium consumption, they ascend 
in the Hertzsprung-Russell (HR) diagram through the red giant branch and evolve to 
the asymptotic giant branch (AGB). It is then possible that a AGB star, consisting of 
a degenerate carbon-oxygen core and an outer helium-burning shell, to evolve into a white 
dwarf star by first cooling down due to neutrino emission and later by surface photon 
emission. 

If axions existed, an additional channel for the cooling of AGB into a white dwarf star 
exists given by the axion bremsstrahlung process 

e + Ze-)- e + Ze + a. (4.41) 

Computing the theoretical luminosity function by taking into account the processes de- 
scribed above and comparing it to the observed cooling rate derived from the measured 
decrease of the rotational period P/P, it is possible to constrain any new physics contribu- 
tion to the cooling process and in particular the contribution from axion bremsstrahlung. 

The constraint on the axion-electron coupling thus obtained is 

Qae < 1.3 X 10-l^ (4.42) 

which set the most stringent bound on the axion-electron coupling constant. We derive a 
constraint on the PQ scale fa by inserting the expression for Qae in Eq. (j4.36p and (j4.38p 
for Cg in the bound in Eq. (j4.42p . to obtain 

/a > 1.3 X 10^ GeV cos^ /3h, (4.43) 

where cot /3h is the ratio of the two Higgs vacuum expectation values. 

Recently, a method for constraining the axion-photon coupling Qa'y using the amount of 
linear polarization in the radiation emerging from magnetic white dwarfs has been proposed 
|71j . This method sets an upper limit on ga-y that depends on the axion mass ma- For 
nia ^ 10~^ eV, the authors in Ref. [71] derive the bound ga-y < 10^^'^ GeV^^: this bound is 
much more stringent than the one obtained using measurements on the lifetime of horizontal 
branch stars in globular clusters, see Sec. 14.5.41 
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4.5.2 Constraints from SN1987A 

A supernova (SN) event of type II consists of a core collapse of a massive star which 
subsequently leads to a proto-neutron star. Axions may be produced in SN events via 
axion-nucleon bremsstrahlung 

N + N^N + N + a, (4.44) 

and this additional process can affect the cooling time of the SN and shorten the duration 
of the burst. The pattern of the cooling time as a function of the axion-neutron coupling 
QaN is described in Ref. |147j . For very low values of the coupling, the axion emission does 
not affect the burst duration. As gaN increases, the burst duration shortens because the 
emission of bremsstrahlung axions increases, reaching its minimum value when the mean 
free path of axions in the medium is of the order of the size of the SN. For even larger 
values the axions are trapped in the medium and the axion emission decreases, until the 
burst duration becomes unaffected by the presence of these new particles. 

In 1987 it was possible to test models of supernovae explosions, due to the observation 
of the emission from a distant supernova named 1987A. The flux of antineutrinos v'e com- 
ing from SN1987A was detected at both Kamiokande II and Irvine-Michigan-Brookhaven 
experiments, allowing us to compare the data with theoretical expectations. Due to this 
analysis, it is possible to exclude axion models with an axion-nucleon coupling in the range 

3 X 10^10 < 5aiv < 3 X 10^^. (4.45) 

This range corresponds to the high yield emission of axions excluded by the statistical 
analysis on the SN1987A data. 

4.5.3 Axions from the Sun 

In the Sun, axions might be produced through the Primakoff process: an incoming 
photon interacts with the electromagnetic field of a nearby electron or nucleus, converting 
into an axion, see the Feynman diagram in Fig. 14.21 Since this new channel can in principle 
shorten the lifetime of the Sun, it is possible to constrain the coupling of axions to photons, 
electrons, and nucleus by imposing that the presence of axions does not spoil the standard 
solar model (SSM). 

To see this in mode details, we consider solar limits on the axion-photon coupling 
ga'y- in the approximation in which the electron or nucleus is nonrelativistic, the rate 
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of conversion of photons of energy E into axions of the same energy in a nondegenerate 
plasma of temperature T is 



(4.46) 



Here, the screening factor kg is given by the Debye-Huckel formula, 



A^? = 4vra^ j^y, + Z| y, j , (4.47) 

where is the nondegenerate baryon density and Y^, Yj are the fractions of electrons and 
nuclear species j per baryon. The energy loss rate in axions per unit volume is then 

/(fi k E 
(2^ g-s/T _ I ^7-^1' (4-48) 

where = 2 and k is the three- momentum of the incoming photon of energy E. Eventually, 
using the solar temperature distribution obtained from the SSM, the flux of solar axions on 
Earth is evaluated as 

La = 1.7 X 10-3 f-^^-^)^ Lo, (4.49) 



.10-10 GeV-i 

with Lq indicating the Sun luminosity. Here, we use the conservative limit on the solar 
axion luminosity La ^ 0.2 Lq [155j which yields 

ga^ < 1.1 X 10-9 GeV-^ (4.50) 

Using the fact that axions affect the sound speed-profile of the Sun, the authors in 
Rcf. [155] obtained a similar bound, 

ga^ < 1.0 X 10-9 GeY-\ (4.51) 

Raffelt and Gondolo ^5] derived a restrictive bound on ga-y using the sensitive dependence 
of the production rate of neutrinos from ^B. Using the measurements of the solar neutrino 
flux by the Sudbury Neutrino Observatory, the authors in Ref. [75] obtain 

< 7x IQ-io GeV-^ (4.52) 



To conclude, we remark here that the best constraint on ga-y from astrophysical measure- 
ments comes from considering the production of axions in globular clusters. 
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4.5.4 Axions and globular clusters 



A globular cluster (GC) is a gravitationally-bound ensemble of stars which formed 
at about the same epoch. Since these stars share the same age, but not other physical 
parameters like the mass or the surface temperature, a GC is particularly suitable for testing 
models of stellar evolution. A typical method used for comparing stars within the same GC 
is to plot the color luminosity vs. the total brightness of each star in a color- magnitude 
diagram. The typical quantities used are the color B — V and the brightness V. 

When axions are included in the model, helium-burning stars may consume helium 
faster than their expected rate because of the extra axion-production channel. The lifetime 
of stars which are in the horizontal branch is reduced by a factor 



due to the production of axions. When a statistically significant ensemble of helium-burning 
stars is considered, it is found that the lifetime of these stars agrees with theoretical 
expectations within 10%, which in turns leads to the bound 



This bound on the axion-photon coupling is a much more stringent bound than those 
obtained from considerations on the solar activity. 



As first proposed by Sikivie |160] . it is possible to directly search for axions in labo- 
ratories by using the conversion of an axion into a photon in the presence of an external 
electromagnetic field. This process, which is the inverse process of the Primakoff effect, 
has been extensively used in a number of experiments. Direct searches of axions fall into 
three primary categories: search of axions of galactic origin (axion haloscope), produced in 
the Sun (axion helioscope), and pure laboratory experiments, in which virtual axions are 
produced by shining lasers in strong magnetic fields. 

4.6.1 Axion haloscope 

Axion haloscope searches rely on the possibility that relic axions from the Big Bang 
would be gravitationally bound to the Milky way, having a nonrelativistic velocity v with 
dispersion Av ~ 10~^ c, corresponding to the rotational velocity of a virialized object in the 




da-y 



< 1 X 10^^" GeV 



(4.53) 



4.6 Direct axion searches 
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vicinity of the solar system. This would correspond to an axion mean energy and energy 
dispersion 

E ^ nia (^l + , and AE = ruaV Av. (4.54) 

In order to detect these particles, a sensitive technique known as the microwave cavity 
has been developed. In such a cavity, a strong electromagnetic field is produced, with a 
frequency related to the size of the cavity. For a given frequency, there exists a narrow range 
of the axion mass for which the axion would interact with the electromagnetic field and 
convert into a light pulse which would be eventually detected by a receiver. In order to 
search for different values of the axion mass, the size of the microwave cavity is adjustable. 
The search of axions through a microwave cavity assumes that the totality of the CDM is 
in the form of axions. As we will explore in depth in this thesis, this is possible if the mass 
of the axion is of the order of the // eV, corresponding to the energy of a microwave with a 
wavelength of approximately 10 cm. 

The probability that an axion of energy E converts into a detectable photon while 
traveling in the homogeneous and transverse magnetic field i? of a microwave cavity is 

^"-^ = l'Z^+tV^) + " ' "^^^ ' ^'-''^ 

where L is the length of the path, V is the inverse absorption length for photons in the 
medium and q is the momentum transfer between the axion and photon. In the vacuum, 
the momentum transfer q reads 

2 2 

^ (4.56) 



2E 

where the result is expressed in terms of the effective photon mass in a plasma with electrons 
density ng, 

(4.57) 

Notice that the result in Eq. (j4.56p is similar to that obtained in the neutrino oscillations 
theory. When the absorption F is neglected, Eq. (|4.55p reduces to 

F„,= (*l^)%m^(S^). (458) 

The first experiments of this type were performed at Brookhaven Laboratories |195j and 
at the University of Florida [82] : the axion mass is excluded from the former and the latter 
experiments in the ranges 5.4 5.9 /ieV and 4.5 16.3 /ieV, respectively. As of today, the 
most sensitive axion haloscope is the Axion Dark Matter eXperiment (AD MX) at Laurence 
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Livermore National Laboratory (LLNL). In its first direct scan of the parameter space, 
ADMX directly excluded axions with mass in the range |12j 

1.9 ^eV <ma< 3.3 fieV. (4.59) 

A later upgrade of the experiment [l3] exploited a Superconducting QUantum Interference 
Device (SQUID), which replaced the microwave receiver and allowed us to extend the scan 
to higher values of the axion mass. The improved ADMX excluded the region [13] 

3.3 ^eV <ma< 3.53 /ieV, (4.60) 

so that the overall region excluded by the ADMX experiment so far is 

1.9 /ieV < rria < 3.53 /ieV. (4.61) 

4.6.2 Axion helioscope 

The Sun produces axions through nuclear interactions within its core; the expected axion 
flux at Earth due to the solar activity is expressed in Eq. (j4.49p . The search for solar axions 
is mainly conducted via two types of experiments, which respectively exploit the intense 
Coulomb field in a crystal to convert axions into photons (Bragg scattering experiments), 
or use intense magnetic fields that point at the Sun in which solar axions can convert into 
photons (magnet helioscope experiments). 

4.6.2.1 Search of solar axions through Bragg scattering 

This type of experiment takes advantage of the fact that the mean energy of solar axions 
is around 4keV, so that the axion wavelength is of the same order as the lattice spacing in 
a typical crystal. The expected Bragg scattering of the resulting photon would increase the 
signal of the axion interaction with the lattice to around 10^ with respect to the interaction 
of an impurity in the same crystal. Moreover, the signal would be distinctive because of 
the relative movement of the Sun within 24 hours (daily modulation of the signal). Various 
experiments undergoing such a search use a different crystal such as Sodium Iodide at 
DAMA [30], Germanium at SOLAX [65] and COSME [132], or Germanium and Sihcon at 
CDMS [6]. All of these experiments constrain the axion-to-photon coupling to about the 
same level, ga-y ^ 2 x 10~^ GeV^^, which however is more than one order of magnitude 
above the bound from haloscope searches. 
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4.6.2.2 Solar cixions through cixion telescopes 

The main component of an axion telescope consists of a magnet in which the strong 
magnetic field can be pointed at the Sun. An axion in the strong magnetic field may 
convert to a low-energy X-ray, through a mechanism similar to the conversion of CDM 
axions in a microwave cavity. In fact, the probability of conversion is given by Eq. (j4.55p 
also in this situation, provided that the same conditions in which the latter equation has 
been derived are met. 

Searches of solar axions via this method were conducted by Lazarus and collaborators 
|115j and by the Tokyo Axion Helioscope collaboration |133| : the latter group constrained 
5a7 < 6 X 10"^*^ GeV~^ for axion masses rua < 0.03 eV. 

A parallel search has been conducted by CERN Axion Solar Telescope (CAST) exper- 
iment, which to date is the most sensitive axion helioscope in use. In fact, CAST in its 
two phases has been able to constrain ga-y more stringently than astrophysical bounds, 
setting < 8.8 x 10~" GeV~^ at 95% C.L. and for < 0.02 eV in Phase I [203], and 
gay < 2.2 X 10~^° GeV~^ at 95% C.L. and for the mass range 0.02 eV < ma < 0.4 eV 
during Phase II 

4.6.3 Production of axions by laser 

Axion searches conducted in laboratories make use of an intense laser beam that might 
partially convert into axions or other pseudo-scalar particles in the presence of a strong 
magnetic field, via the Primakoff effect. We briefiy discuss two different types of experiments 
that exploit this technique. 

4.6.3.1 Polarization experiments 

When axions are produced by the interaction of a strong laser beam with an external 
magnetic field, the polarization of the laser beam is affected because of such a nonzero 
coupling with the pseudo-scalar particles. Moreover, when a photon-axion-photons re- 
conversion takes place, the emerging photon has a component of the E and B fields 
retarded with respected to a freely-propagating photon. This signal had been claimed 
by the Polarizzazione del Vuoto LASer (PVLAS) collaboration [35], but the event was later 
labeled as artificial. 
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4.6.3.2 Regeneration experiments 

When a polarized laser beam propagating in a constant transverse magnetic field is 
blocked by a thick absorber (a wall), it is possible to detect photons of the same wavelength 
as the laser beam on the other side of the blocking medium itself. This process is called 
photon regeneration, and the experiment is colloquially referred to as the shining light 
through a wall experiment. This unusual phenomenon is possible if some photons from the 
laser beam convert into axions in the magnetic field, before being absorbed by the wall. In 
this case, the produced axions are able to travel through the wall with little absorption. 
Applying a second magnetic field on the other side of the wall makes it possible for some 
axions to reconvert into photons through the inverse Primakoff process. 

This type of experiments has been conducted by the Brookhaven-Fermilab-Rutherford- 
Trieste collaboration [34] and at CERN with the Optical Search for QED vacuum birefrin- 
gence, Axions and photon Regeneration (OSQAR) experiment. In addition, the shining light 
through a wall experiment is able to put limits on more general models in which photons 
couple to axion-like particles (ALPs), theoretical particles that share similar properties with 
the axion but for which the mass-energy scale relation in Eq. (j4.12p does not apply. We will 
talk in detail about ALPs in Sec. [3 The bound on the ALP-photon coupling constant Ga-y 
from this experiment is 

Ga-y < 1 X 10^1° GeV"^ (4.62) 



CHAPTER 5 



REVISING THE AXION AS THE COLD 
DARK MATTER 

5.1 Introduction 

The recent measurements by the WMAP mission |109| have estabhshed the relative 
abundance of dark and baryonic matter in our universe with great precision. About 84% 
of the nonrelativistic content in the universe is in the form of CDM, whose composition is 
yet unknown. One of the most promising hypothetical particles proposed for solving the 
enigma of the dark matter nature is the axion |190| I194j . This particle was first considered 
in 1977 by R. Peccei and H. Quinn |142j in their proposal to solve the strong-CP problem of 
the QCD theory. Although the original PQ axion is by now excluded, other axion models 
are still viable [Ml ^57\ [5311202]. 

The hypothesis that the axion can be the dark matter particle has been studied in 
various papers (see e.g. [HSl [21 [53l [Ml [IMl [l23l [Ml [88] and the reviews in [Ml [M]). 
Here we examine the possibility that the invisible axion may account for the totality of 
the observed CDM, in light of the WMAP5 mission, baryon acoustic oscillations (BAO), 
and supernovae (SN) data. We also upgrade the treatment of anharmonicities in the axion 
potential, which we find important in certain cases. We consider invisible axion models, in 
which the breaking scale of the PQ symmetry fa is well above the electroweak scale. The 
axion parameter space is described by three parameters, the PQ energy scale fa, the Hubble 
parameter at the end of inflation Hj, and the axion initial misalignment angle 9i. 

5.2 Production of axions in the early universe 
5.2.1 Thermal production of axions 

In the early universe, a population of thermal axions originates together with standard 
model particles and possibly other exotic components. Scattering and annihilation processes 
keep this population of axions in thermal equilibrium with the remaining of the hot plasma 
at temperature T, with the axion number density n^^ following the Boltzmann equation 
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^> + 3if(t) n*^ = r,,„ ( - n*M . (5.1) 



dt 

Here, the number density of axions at thermal equilibrium is, see Eq. (j2.75p . 

^eq,th ^ C(3) .5 2) 

In Eq. (jS.ip . the rate at which axions annihilate and are created in the plasma is 

^ann = ^ ni{aiv), (5.3) 
i 

where Ui is the number density of the particle species i. The index i runs over processes 
of the type a + i ^ other particles, with cross section fjj, and with v the relative velocity 
between i and the axion. Finally, the angular brackets in Eq. (j5.3p indicate an average 
over the momentum distributions of the particles involved. Using the conservation of the 
number density at equilibrium, 



J eq,th 
u iLa 

Jt 

we can rewrite Eq. (j5.ip as 



+ 3H{t) nl'i^'^ = 0, (5.4) 



di 



a'{t) (^'i'*^ - n^^)] = -ra„na3(t) {nl^'''' - nf) . (5.5) 



The solution of Eq. ()5.5p implies that a thermal population of axions reaches its equilibrium 
value n^'^^ exponentially fast whenever 

Tann > H{t). (5.6) 

The time tchcm at which axions chemically decouple from the plasma is then defined with 

Eann — -f^(^chcm)- 

Depending on the axion model, axions couple differently to standard model fermions, 
photons, and exotic particles, whereas the coupling to gluons is model-independent and is 
given in Eq. (j4.ip . Thermal axions were first studied by Turner |179| . who considered the 
Primakoff and photoproduction processes. Given a heavy quark Q of mass mq, whenever 
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the temperature of the plasma T < niQ the Primakoff process dominates and axions are 
created via q + ^^q + a{q is a Ught quark) , with cross section 

CTPrimakoff = "Tf ^^''^^ 
J a 

At higher temperatures T > niQ, photoproduction of axions via Q + ^ ^ Q + a dominates 
with cross section 

c^7-prod = a (^YY^ ■ (5-8) 

Ref. jl79j concluded that a thermal population of axions exists today whenever fa ^ 
10^ GeV; moreover, for fa^^>^ 10^ GeV, the thermal population is greater than the 
nonthermal. 

These considerations were revised in Ref. |125j . where the authors extended the work 
in Ref. |179j by including other processes than Primakoff and photoproduction, using the 
model-independent axion-gluon vertex derived from Eq. ()4.ip . In detail, the authors in 
Ref. [125] include the following reactions: 

a + g^ q + q 

a+q^g+q 
a+g^g+g 

These processes have a cross section of the order of Uagg = a^^/ f"^- With these additional 
reactions included, axions thermalize if < 10^^ GeV |125j . which is a broader condition 
than the one found in Ref. |179] because axions are tighter bound to the primordial plasma. 

Assuming that axions thermalize at temperature Tth, their number density today is 
obtained by assuming the conservation of the comoving number density, 

<^To) = <\T,,) ' = 7.5 cm-3 (5.9) 

where g^s{T) denotes the relativistic number of degrees of freedom at temperature T. This 
thermal population of axions is very diluted today, once compared for example to the CMB 
photon number density n^(To) 410 cm~'^. 
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5.2.2 Axions from the misalignment mechanism 
5.2.2.1 Equation of motion on the FRW metric 

The generic Friedmann-Robertson- Walker (FRW) metric is written as 

ds"^ = df - a'^{t)dx- dx, (5.10) 

where a{t) is the scale factor, which relates to the Hubble expansion rates by H(t) = 
d{t)/a{t), and x are comoving spatial coordinates. On this metric, the angular variable 
6 = 9{x) has equation of motion 

e{x) + ^V'{6) = 6 + 3H{t) e - -^Vl 0+^Ve = 0, (5.11) 

Ja ^ V') Ja 

where a dot indicates a derivative with respect to t and Vq = dV/dO. Using the axion 
potential in Eq. (|4.7|) 

y(0) = m2(r)/2(l-cos0), (5.12) 

that originates from nonperturbative effects after the spontaneous breaking of the ?7(1)pq 
symmetry, we obtain 

e + 3H{t) 9 - -^rr'^l ^ + sin 6 = 0. (5.13) 

a [t) 

From now on, we only focus on the mode with zero momentum, for which 9 = and 
which constitutes the CDM component. A general treatment on axionic modes with nonzero 
momenta can be found in Ref. |161j . For zero modes, it is 

9 + 3H{t) 9 + ml{T) sine = 0. (5.14) 

In the early universe, when the plasma temperature is T > Aqcdj the axion mass and hence 
the potential is negligible, see Eq. ()4.29p . and a solution to Eq. ()5.14p is 6* = 9i =constant, 
with 9i being the initial value of the misalignment field. As the temperature drops below 
Aqcd, or, more precisely, when the axion field starts to oscillate at a temperature Ti = 
O(GeV) given by [178] 

3H{Ti)^ma{Ti), (5.15) 

the mass term in Eq. (I5.14P becomes important and the solution is no longer trivial. 

The axion field starts oscillating in the nonperturbative potential V^(^) with initial 
condition 9 = 9i. We can solve the equation above analytically by assuming that the 
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axion mass nia is constant and that 9i be much smaller than one, so that sin^ ~ at all 
times. For a generic cosmology in which the scale factor depends on time as 



a{t) ~ t^, 



and thus the Hubble parameter is 



H{t) 



a{t) _ (3 
a(t) ~ I' 



we obtain 



(5.16) 



(5.17) 



(5.18) 



with Ta = mat and where Ji,{x) is the Bessel function of the first kind of order v in the 
variable x and T[x) is the Euler Gamma function. 

The general solution to Eq. (I5.14P that does not assume Oi and takes into account a 
generic temperature-dependent mass ma{T) can be found only numerically. In Fig. 15.11 we 
compare the numerical solution of Eq. (j5.14p (solid red line) with the approximate solution 
expressed in Eq. (j5.18p (blue dashed line) for 9i = 2. We see that anharmonicities sensibly 
modify the behavior of the field dynamics. 




T = m„ t 



Figure 5.1. Comparison between the numerical solution of Eq. (j5.14p (red solid line) and 
the approximate solution given in Eq. ()5.18p (blue dashed line), with initial condition 9i = 2. 
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5.2.2.2 Computing the oscillation temperature 

When the potential term in Eq. (|5.14|) can no longer be neglected, coherent oscillations 
in the axion field begin. This happens when the temperature of the universe drops to the 
value Tf^'^, given implicitly by Eq. (I5.15P : we have added a label "std" that refers to the 
fact that this temperature is computed in the standard cosmological scenario, in which the 
universe is radiation-dominated right after the end of inflation. For T > A, we obtain an 
equation involving Tf^'^ by inserting Eq. (|4.16p in Eq. (j5.15p and using Eq. (j2.87p for the 
Hubble expansion rate in a radiation-dominated universe. 

In this way, we find that the parameter 

XI ^ -j^, (5.20) 

must satisfy the expression 



rpStd 



""^°"^"^'^-° ) =(607±27)MeV/„-;,f, T > Aqcd, 



(5.22) 



, f ^ =(3300± 1700)/-i^2 J ■ . V^std, forxi>l. (5.21) 
(1 + lnxi)^-^ \j sm9ig^{Tl^'') 

Eq. (|5.2ip defines implicitly the value of Tl^'^. We have checked numerically that the 
logarithmic term contributes a 15% contribution to the final solution, so we neglect it 
because of the huge uncertainty in the parameter a in Eq. (|4.28p . Neglecting the logarithmic 
term in Eq. ()5.2ip is equal to solving Eq. ()5.15p with the axion mass given by Eq. (|4.28p 
instead of Eq. (jdUj) . We obtain [TH5] 

1/6 

4.98 y/g-^Tf^) 

\ 1/2 

— "^''^^P' = (67 ± 3) MeV /"^/^ T < Ac 

In the last expression, we have used the numerical value a = 0.017 ± 0.003 and the values 
for g^{T) in Eq. I^M\i . 

5.2.2.3 Axion energy density from vacuum realignment 

A generic scalar field $(x) of mass m<j, in a quadratic potential has energy density 

p^ = ^ (ml<^^ + ^^y (5.23) 

In the case of the axion, the scalar field a{x) = faO/N does not move in a quadratic potential 
but rather in the residual potential given by Eq. (|5.12p . This sinusoidal potential behaves 



»-QCD, 
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like a quadratic potential if 0j ^ 1, but when the initial value of the misalignment angle is 
of order one, anharmonicities cannot be neglected, as shown in Fig. 15.11 

To begin with, we will consider the case in which the initial value of the misalignment 
angle is small, which allows the expansion 



V{e)^-mi{T) [i^^j for0,«l. (5.24) 

This potential turns on at the temperature Tf''^ defined by Eq. (|5.15|) : for higher values of 
T, the angular variable 9 is frozen to its initial value Oi, so that we can assume 9 = 9i when 
T = Tf^'^. Then, the axion number density is initially 

ua (rf<i) = IxmaiTn (1)' {9^). (5.25) 

Here, the factor x models the temperature dependence of the axion mass around [178], 
and depends on the number of quark flavors Nj that are relativistic at temperature Tf^'^ 
through the index c appearing in Eq. (j4.29p as 

21 + Nf 

X = 0.44 + 0.25 c = 0.44 + (5.26) 

24 

Here, we consider Nf = 3 because the axion field starts to oscillate at temperatures below 
the GeV, so that we take x = 1-44. The number density at the present time is found by 
imposing the conservation of the comoving axion number density Eq. (j2.83p . in the form 
HDH]: 

'"'^^^^^-0, (5.27) 



s{T) 

where the entropy density is given by Eq. (j2.79p . 

s{T) = ^gMT)T\ (5.28) 

and g*s{T) denotes the entropy number of degrees of freedom at temperature T, see 
Eq. (I2.80p . We approximate g*s(T) = g*{T) when T > 4 MeV, while the present value 
for g*s{To) = 3.91 differs from g*{To) = 3.36 |108| . For the range of interest here, for the 
relativistic degrees of freedom, we use the values in Eq. (|2.86p . Using Eqs. ()5.25p . (j5.27p 
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and ()5.28p . the present axion energy density from the misahgnment mechanism p^^^(Tq) 
maUaiTo) is 



Prin) = ]-xmama{Tt'')^^ (^] {ef). (5.29) 



2 

Dividing the last equation by the critical density pc = ?>HqMpJSit and using Eq. (|4.29p for 
ma{Tl^^), the cosmologically relevant ratio r^™'' = Pa/ Pc is 



2pe %(rn ^ 1^ T<Aqcd. 

For the temperature Tf^'^, we substitute the expression we obtained in Eq. (j5.22p in Eq. (j5.30p 
with the numerical values x = 1-44 and a = 0.017 it 0.003, obtaining 

^mis^2 ^ f (0.23 ± 0.03) {ef) fl%, fa < fa, ,^ 

\(6.1±0.8)X 10-3(^2^^3/2^ 

The parameter fa is defined as the energy scale at which the two expressions in Eq. (]5.3ip 
match: 

fa = 5.4 X 10^^ GeV. (5.32) 

The two cases in Eq. ()5.3ip reflect the dependence of the axion mass ma{T) on the temper- 
ature T in Eq. (029]). 

5.2.2.4 The role of anharmonicities 

In the general case in which 0j ~ 1, the expression in Eq. (j5.3ip underestimates the 
actual axion energy density from the misalignment mechanism. In fact, the average value 
of of over a Hubble horizon when the axion potential is treated as a pure harmonic potential 
like in Eq. ()5.24p gives 

(^.'> = ^ f^de,6f = ^. (5.33) 

Instead, the potential in Eq. (I5.12p is not harmonic, and the effects when the value of the 
angle 6i lies close to vr become rather important. Turner |178j found that the effect of these 
anharmonicities raise the value of the integral in Eq. (j5.33p by a factor 1.2-1.6. 

In order to take into account these effects, we parametrize anharmonicities with a func- 
tion f{6i) that behaves like f{Oi) — )• 1 for 0j — )■ and diverges when 9i — )• vr. This function 
has been considered by various authors. Turner |178j integrated Eq. (I5.14p numerically 
and described how the anharmonicity factor can be computed, but did not give an explicit 
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formula. Lyth |122j followed the idea by Turner and performed an explicit calculation, 
obtaining the behavior 

1.175 



fi9^ 



In 1 



TT 



for > 0.97r. 



(5.34) 



Lyth also commented that his result differs from Turner's by a factor of two. The exponent 
1.175 crucially depends on the dependence of the mass on temperature, see Eq. (j4.29p : in 
fact, Lyth assumed that ma{T) oc T~^''^ . 

Strobl and Weiler |169| and Bae et al. |14j performed a more precise numerical analysis 
following the methods in Ref. |178j . and confirmed the result in Ref. |122| for the behavior 
of f{6i) around 6i = vr. 

Here, we look for an analytic function f{Oi) that extends Lyths formula in Eq. (j5.34p 
to values of < 0.97r and symmetrically to negative values of Oi. Such function must show 
the limiting behaviors 



lim f{ei) = 1, and lim /(^i) = +oo. 



Moreover, the integral of the function must give 



1.4 

3 ' 



(5.35) 



(5.36) 



where the extra factor 1.4 accounts for Turner's result in Ref. |178j . We parametrize the 
anharmonicity function as 



In- 



(5.37) 



[1 - {ed^?r 

where the factor 7 is fixed by the exponent c appearing in the temperature dependence of 
the mass ma{T) oc T^^ via the formula 

' c + 2 33 + iV/ ^ ' 

and Lo is fixed numerically by imposing the condition in Eq. (j5.36p . When Nf = 3 and thus 
c = 4, we obtain 



f{Oi 



In- 



7/6 



(5.39) 



This expression for f{6i) is analytic on [— 7r,7r] and has the correct asymptotic behavior in 
Eq. dSSil). 
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5.2.3 Axions from string decays 
5.2.3.1 The domain wall problem 

Since 9 has period 2TrN, the potential in Eq. (j5.12p has N distinct minima: when the 
axion field oscillates around one of these minima, the expectation value of the misalignment 
angle tends to zero and the CP violation in the QCD sector lies to experimentally acceptable 



values. However, as first noted by Sikivie }159j , when > 1 there is a domain wall problem, 
that is, when the PQ symmetry breaks in the early universe the axion field takes different 
minima in causally disconnected regions, in contrast with standard cosmology |201j . In order 
to avoid this problem, various models have been implemented in which = 1 [1141 169 1 l^lj- 
If > 1, we can avoid the problem by assuming that the PQ symmetry breaks before or 
during inflation: a single domain with the same value of the PQ vacuum would be inflated so 
that the region enclosing our visible universe would lie within such a bubble. An alternative 
solution to the domain wall problem in the nontrivial case A^ > 1 consists on introducing a 
small breaking of the PQ symmetry |159j . so that the different vacua are slightly tilted and 
tunneling of the axion field singles out the energetically favored lowest true vacuum. If this 
is the case, the axion walls predominantly decay into gravitational radiation [36] . 

5.2.3.2 Parameter space of cixionic strings 

The spontaneous breaking of the Peccei-Quinn symmetry leads to the formation of 
topological defects such as axionic strings |183j . In fact, if the PQ symmetry breaks after 
inflation, cold axions produced by the decay of topological strings in the early universe 
contribute a large fraction of the axion energy density. Instead, if the PQ symmetry breaks 
before or during inflation, topological defects are washed out by inflation so axions from 
axionic strings are not present. 

The present axion energy density produced by axionic strings p^'^{To) is proportional 
to the present axion energy density produced by the misalignment mechanism /3™"'(To) 
[Ml Ell [80]. The ratio Q between pT{To) and pT^Tq) can be put in the form 



Q.^ = ii?l. (5.40) 

Here, following the notation in Ref. [8^ (however, we use C for their parameter x to avoid 
confusion with our x in Eq. (I5.25P ): is the number of degenerate QCD vacua, f is the 
factor by which the axion comoving number density increases due to string decays, averaged 
over all possible processes that convert strings to axions, ^ is a constant factor depending 
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on the string network model, and ^ accounts for the uncertainties in the low-energy cutoff 
of the radiated axion field. In the following, we add a suffix "std" when these parameters 
refer to the value in the standard cosmology. In the following, we discuss each parameter in 
Eq. (|5.40p separately in the standard cosmology, leaving the extension of previous theoretical 
results in generic nonstandard cosmologies to Sec. 16.31 

The values of the parameters f^'*^, C^^'^j ^-nd i^'^*'^ have been discussed extensively 
both theoretically and via numerical simulations of string networks [461 Wl\ [86l UU \^5\ [36} 
Ell EOl ESa [231 [221 [EQ [Ml [ig [H] . The number of QCD vacua is usually assumed Nd = l, 
because a higher value A^^ > 1 might lead to a domain wall problem |159j : here, we assume 
Nd = 1 as well. In the standard cosmology, C ~ !> although the theoretical uncertainty on 
( is around 50% (SOj; here, we take C = 1 without including uncertainties. 

In the literature, there is still disagreement about the numerical values of f^^'^ and 
Here, we discuss results from different authors separately. The value of f depends on the 
details of the axionic string relaxation toward lower energy configurations and on the energy 
spectrum of the radiated axions |156] . In the standard cosmology it is (see Ref . [80] ) 



Here, ti is the time at which the axion field starts to oscillate and 6 is the string core size 
[36j . In Eq. ()5.4ip . the first line corresponds to the string emission model in Refs. [4611471 I^Hl 
[13 [231 [221 [ESI [ES] , while the second line corresponds to the model in Refs. [36l [86l IHTl [80] . 
In the first line in Eq. (j5.4ip . we have obtained the value f^*'^ ~ 70 by considering the 
illustrative case 6 = (10^^ GeV)~^ and Tf^'^ = 1 GeV: this is approximately the same value 
found in Refs. [Ml [IZl [Ml [IS [231 [22l [EH [ES]. We remind here that in the standard 
cosmology, the time t relates to the temperature T of the universe through the Friedmann 
equation l/2t = H{T). On the value of numerical simulations for an evolving string 
network in Ref. [25l[8] obtain f ~ 13, while simulations in Refs. [8ll[80l[E6] yield ^"^'^ ~ 1. 

In this chapter, when showing numerical results, we will use the numerical values f^^'^ = 
0.8 and ^'^^'^ = 1, corresponding to the string model studied by Harari, Hagmann, Chang, 
and Sikivie. 




\n{ti/5) 70, for a slow-oscillating string, 
0.8, for a fast-oscillating string. 



(5.41) 
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5.3 Axion isocurvature fluctuations 

If the PQ symmetry breaks before inflation ends and if the maximum temperature after 
inflation is less than /„, the axion field is present during inflation and is subject to quantum 
mechanical fluctuations. Since at such high temperatures the axion is massless, the variance 
of quantum fluctuations 6a{x) in the axion scalar field a{x) is given by Eq. (j2.54p . which 
we write in the form 

= (^) • (5.42) 
It follows that fluctuations in the misalignment angle field 9{x) = a{x)/fa have variance 

.|^(^)^. (5.43) 

These fluctuations are over the initial misalignment angle 9i , which has a single value within 
one Hubble volume since it was causally connected at the onset of production. 

Under the condition in Eq. (j5.54p . axion isocurvature perturbations are present during 
inflation and are constrained by WMAP-7. As originally shown in Ref. |124j . this type of 
measurement can also set an upper limit on the inflationary scale Hj. Here, we revise this 
computation using the latest measurements in Ref. [llOj . Defining the power spectrum of 
axion perturbations A^(fc) = {\6pa/Pa\'^), one finds 

The axion entropy-to-curvature perturbation ratio is then 



or, introducing the axion adiabaticity ao(^o) like in Refs. [109^ 1110] . 

Al{ko) ao{ko) 



(5.45) 



A^{ko) 1 - ao{k 







(5.46) 



The adiabaticity is constrained by the WMAP-7 + BAO + SN data to be [TTO] 

ao < 0.064 at 95% CL. (5.47) 



Using the central value for A'^{kQ) in Eq. ()2.52p . this bound is rephrased as 



(5.48) 
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5.4 Parameter space of the cosmological axion 

Depending wether the PQ symmetry breaks during inflation or after it and wether 
such symmetry is even restored after inflation ends, the cosmological history of the axion 
field greatly differs. In fact, if the axion energy scale is lower than the Gibbons-Hawking 
temperature Tqh = Hi/2t: [70] . 

fa < ^, (5.49) 

the PQ scale breaks after inflation and the axion field does not participate in inflation. 
Furthermore, if the axion energy scale lies below the maximum temperature that the 
universe reaches after inflation ends due to reheating, 

fa < TuAX, (5.50) 

the PQ symmetry is restored after inflation, to be broken again when the universe subse- 
quently cools down in a postinflationary scenario. We remind that Tmax is related to the 
reheating temperature Trh by [42] 

Tmax ~ (T^h Hi Mpy) . (5.51) 

When either Eq. (j5.49p or Eq. (j5.50p is satisfied, the axion field exhibits specific features: 
different values of the misalignment angle Oiix) are present within one Hubble volume, 
giving rise to fluctuations that are adiabatic as observed in the CMB spectrum, and there 
are no nonadiabatic fluctuations because Eq. (|2.54|) does not apply. Moreover, topological 
defects like axionic strings are present because inflation could not wash them out or because 
the defects generated after the PQ symmetry were restored. We refer to this portion of the 
parameter space as Scenario I. 

On the contrary, two necessary conditions for the axion field to show nonadiabaticities 
have to be met: (a) the PQ symmetry breaks before inflation ends and (b) it is not restored 
afterwards. Condition (a) requires 

fa > ^, (5.52) 



while condition (b) requires 



fa > Tmax, (5.53) 



When both these conditions are met, the PQ symmetry breaks before inflation ends and it 
is never restored afterwards: we call this portion of the parameter space the Scenario II. 
Under these conditions, the salient features of the axion field consist in that nonadiabatic 
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fluctuations are developed because the axion participates in inflation, the value of the initial 
misalignment field has a single value Oi over one Hubble volume, corresponding to the 
inflated domain we live in, and no defects are present. Putting together the two requirements 
for Scenario II, we obtain that nonadiabatic axion fluctuations arise when 



/a > max{rG//,rMAx}- (5.54) 

Here, we assume Tmax < Hi/2t^, which leads to Trh ^ 10^^ GeV. This requirement is 
consistent with that obtained from theories in which the inflaton decays into fermions, 
Trh < 10^2 QgY [X20]. Because of our choice, Eq. (ISliOll suffices to identify Scenario I, 
while Eq. (j5.52p describes the parameter space of Scenario II. For the reasons explained 
above, axionic strings in Sec. 15.2.31 contribute to the axion CDM energy density only in 
Scenario I, and the bound computed in Sec. 15.31 applies only in Scenario II. 



5.5 The axion as the CDM particle 

The major interest for axions in astrophysics consists in the fact that these particles 
may account alone for all of the observed CDM. If this is true, axions must be in highly 
nonthermal equilibrium and form a Bose-Einstein condensate [161^ I162j . We explored the 
two leading mechanism that produce an axion population with zero momentum, namely 
the misalignment production in Sec. 15.2.2.31 (see also Refs. [145^ [21 [53| 11671 1178| ) and the 
axionic string decay production in Sec. 15.2.3] (see Refs. [l6l [HI [361 ESI [81] ) • The population 
of thermal axions discussed in Sec. l5.2.T] does not constitute a fraction of the CDM, because 
thermal axions are still relativistic today, contributing to the HDM population. 

5.5.1 Axion CDM in Scenario I 

We first explore the expression for the energy density of CDM axions in Scenario I, 
fa < Hi 1271. As discussed in Sec. 15. 4^ in this region of the parameter space, there are no 
axion isocurvature fluctuations, but the energy density gets contributions from both the 
misalignment and string decay mechanisms. We have evaluated the energy density from 
misalignment mechanism p™"^ in units of the critical density in Eq. ()5.3ip which, when 
anharmonicities are included with the function f{Oi) in the expectation value of Of, reads 

^ f (0.23 ± 0.03) fie,)) fl%, fa < L 

\{6.i±o.8)xio~'{eff{e,))fS„ fa>L 
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with fa = 5.4 X 10^^ GeV. We have seen in Sec. 15.2.2.41 that in Scenario I, the average of 
of over its possible values on a Hubble horizon gives 

(02/(^,)) = 1.4^ =4.61. (5.56) 
We thus use this value in Eq. ()5.55p to obtain 

\(0.028±0.004)/„'g, /„>/„. 

The relative contribution from axionic strings to the axion energy density has been consid- 
ered in Sec. 15.2.31 and it is summarized in Eq. (|5.40p written here as 

nfh^ = nf'h^ (5.58) 

These parameters depend on the model used to describe the axionic string oscillation and 
the radiation spectrum of the radiated axions. Using the parameters found for the scenario 
of a slow-oscillating axionic string first described by Davis, Shellard, and Battye |481 H6| 
l47t [23l [22] (from here on DSB Scenario) or, alternatively, the parameters computed in a 
fast-oscillating string scenario by Harari, Hagmann, Chang, and Sikivie [861 \8T\ [36| [80] 
(HHCS Scenario) we have 

^,,,,2 idlOn'^'^'h'^ DSB Scenario; 

\0.8nT''h^ HHCS Scenario. 



The total energy density = Q^^^h'^ is thus 



f^tot;^2 ^ f 911 ^rh' DSB Scenario; 

\1.8nT^h^ HHCS Scenario. 



a 



Setting the total axion energy density in Eq. (j5.60p equal to the measured CDM energy 
density, 

ncDMh'^ = 0.1126 ± 0.0036 at 68% CL, (5.61) 

we compute the value of the PQ scale f^^'^ for which the axion condensate accounts for the 
observed CDM in the universe in the standard radiation-dominated cosmology, 

ystd _ / (4-3 ± 0.7) X 10^ GeV, for the DSB Scenario; 
~ \ {8.8 ± 1.5) X 10^° GeV, for the HHCS Scenario. 

In computing /q**^, we have only considered the case f^^'^ < fa, because the PQ scale fa^'^ is 
generically bound on Scenario I by /^*^ < and the breaking scale / lies above such bound. 
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Thus, we do not need to consider the possibihty that f^'^ > fa here. Notice that the 
two energy scales in Eq. (j5.62p differ by two orders of magnitude, because of the different 
contribution from axionic strings in the two scenarios. In fact, in order for the total axion 
energy density to equate the CDM energy density, the energy scale f^^'^ and the prefactor 
from axionic strings are inversely proportional; since cold axions are copiously produced 
in the DSB Scenario, the energy scale f^^'^ is much lower in the DSB Scenario than in the 
HHCS Scenario. We also remark that the result 

is one order of magnitude lower than the usually quoted value for the PQ energy scale, 

one order of magnitude comes from equating $7*°* with 
the precise measurement in Eq. (j5.6ip for rJcDM) that is 0,a ~ 0.1, whereas seminal works 
in axion cosmology computed f^^'^ from the equation $7^ = 1 that assures that axions do 
not overdose the universe. Also, in the computation of ^la we included the contribution 
from axionic strings decay and the careful derivation of the anharmonicities effects. 

5.5.2 Axion CDM in Scenario II 

We now discuss the parameter space of CDM axions when fa > Hi/2tt, namely Scenario 
II. We have seen in Sec. 15.41 that fluctuations in the axion field arise with variance given in 
Eq. (15.430 : using the definition of variance, we thus write 

{ef) = {e,f + al = ef + al (5.64) 

When the anharmonicity function f{Oi) is included, we approximate {9f f{Oi)) ~ {Of) f{9i): 
this is a good approximation in the region of the parameter space that concern Scenario 
II, because in this region the term that contains an average over a Hubble volume is much 
smaller than the isocurvature fluctuations term, 9f ^ Cg, and the term is multiplied by 
the anharmonicity function with no further averaging involved. In this approximation, it is 
then 

{eff{6.^) = {6f + al) m), (5.65) 

and the total axion energy density is 

^tot^2 ^ f (0.23 ± 0.03) {ef + aj) /(0.) fi;,% fa < L 

\ (6.1 ± 0.8) X 10-3 {ef + al) f{e,) fl% fa>fa. 

In Eq. (j5.66p . the total axion energy density equates the density from the misalignment 

mechanism because no contributions from axionic strings are present in Scenario II. A 
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peculiar fact about Scenario II is that, in this region of the parameter space, the PQ scale 
can evade the bound fa ^ 10^^ GeV when the initial value of the misalignment field is 
smah, Oi < 1 [TT9] . 

5.6 Constraining the parameter space of the axion 

The axion parameter space is labeled by the axion energy scale fa, the initial misalign- 
ment angle 9i, and the Hubble parameter during inflation Hj. Results valid within the 
HHCS axionic string model are shown in Fig. 15.21 

The region labeled "Tensor modes" shows the constraint on Hj in Eq. (j2.65p coming 
from the WMAP7 plus BAG and SN observations [TTO] . The newly launched PLANCK 
satellite will improve the actual measurement on the tensor-to-scalar ratio r by at least one 
order of magnitude, see Ref. |185| . If PLANCK will not detect gravitational waves and sets 
a new limit r < 0.02, Hf will be bound by Hj ^ 2 x 10^^ GeV. This forecast measurement 
is shown by a vertical dashed line labeled "PLANCK" . 

The region labeled "White dwarfs cooling time" is excluded because there one would 
have an excessively small cooling time in white dwarfs. We have used the bound in Eq. (I4.43P 
with cos^ = 1 , giving the bound 

fa > 1.3 X 10^ GeV, (5.67) 

corresponding to an upper bound for the axion mass nia < 15 meV. We remark here that 
this bound on the PQ scale is slightly below the constraint often used > 4 x 10^ GeV. 
The line 

fa = Hi/2tt (5.68) 

divides the region in which the PQ symmetry breaks after the inflationary epoch has ended, 
or fa < Hi/2-K (Scenario I), from the region where the axion field is present during inflation, 
fa > Hi /In (Scenario II). 

The region marked as ADMXI has been excluded by a direct search of axions CDM in 
a Sikivie-type microwave cavity detector [160| by the ADMX Phase I experiment [12\ I58j. 
The window shown corresponds to 1.9 /ieV < nia < 3.3 /ueV, see Eq. ()4.59p . valid for the 
KSVZ axion model. A narrower DFSV axionic window, not shown in Fig. 15.21 has also 
been ruled out. The dashed line labeled "ADMXII" shows the forecast axionic region to be 
probed in the ADMX Phase II, which would search for axions with mass up to 10/ieV, or 
fa = 6.2 X IQii GeV. The proposed CARRACK II experiment is a cavity search that will 
look for axions with mass in the range (2 50)^eV jl99j . 
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Hj [GeV] 



Figure 5.2. Region of axion parameter space where the axion is 100% of the CDM. The 
axion mass scale on the right corresponds to Eq. (j4.12p with U{1)pq color anomaly = 1. 
When the PQ symmetry breaks after inflation {fa < Hj/2tt), the axion is the CDM particle 
if fa = (8.8ibl.5) X 10^° GeV, or lUa = (67ibl7) fJ-eV, which is the narrow horizontal window 
shown on the right. If the axion is present during inflation (/„ > Hj/2Tr), axion isocm'vature 
perturbations constrain the parameter space to the region on the top left, which is marked 
by the values of 9i necessary to obtain 100% of the CDM density. Other bounds indicated 
in the figure come from astrophysical observations of white dwarfs cooling times and the 
nonobservation of tensor modes in the Cosmic Microwave Background fluctuations. Dashed 
lines and arrows indicate the future reach of the PLANCK satellite and the ADMX and 
CARRACK microwave cavity searches. 



For Scenario I, since {0ff{6i)) = 1Att^/3 and does not depend on Hj or 9i, the value of 
$7™"' depends on fa only and it is given by Eq. (j5.57p . In the HHCS model and for fa < fa, 
we have obtained 

17*°t/i2 = (1.90 ±0.25) /J(2' for the HHCS string model. (5.69) 

Demanding that $7*°* accounts for the totality of the CDM, and indicating with /™*^^ the 
value of the PQ scale for which the axion is 100% CDM in the HHCS string model, we have 
obtained 
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/HHCS ^ 8 ^ 5) ^ ;^q10 QgY^ (5 70) 

The band corresponding to f^^^^ is drawn in Fig. 15.21 as the horizontal window on the 
lower right. Assuming a U{1)pq color anomaly = 1, the = ilcDM band corresponds 
to an axion mass 

ma = (67 ± 17) /ieV. (5.71) 

In Scenario II, axion isocurvature fluctuations are present and lead to the bound in 
Eq. (|5.48p . Together with the expression for ^2™*^ in Eq. (15.66^ and the condition that axions 
constitute 100% of the observed CDM, the adiabaticity bound excludes the shaded region in 
the center of Fig l5.2[ The PLANCK satellite is expected to improve the current bounds on 
the axion isocurvature fluctuations by at least one order of magnitude. The dashed line on 
the left of Fig. 15.21 shows the new bound on the allowed region when Eq. (|5.47p is replaced 
by ao < 7 X 10^^. 

The leftmost boundary of this region contains two kinks and can be approximated with 
'7.9 X 10^ GeV , if/ < 9.6 x 10*^ GeV, 

10» GcV^ 

H, 
IQiu GeV^ 

The upper kink occurs at = fa and is due to the change in the dependence of the axion 
mass on the temperature, Eq. (j4.29p . The lower and smoother kink around fa ~ lO^^GeV 
arises from the fact that the anharmonicity function f{9i) differs from one at values of fa 
smaller than 10^^ GeV, see Eq. (I5.39p . Notice that the simple proportionality fa oc Hi at 
small Hi < 10'' GeV is independent of the detailed form assumed for the function f{9i) 
near 6i = vr, and derives in a straightforward way from Eq. (j5.48p . 

In the remaining region on the left of Fig. 15.21 the axion can be 100% of the CDM, 
provided the value of the initial misalignment angle 9i is chosen appropriately. The Oi 
contours in Fig. 15.21 indicate the appropriate values of 9i for a given Hi and fa- Notice that 
the 6i contours are horizontal, i.e. independent of Hi, since in that region the contribution 
from fjg in Eq. (I5.43P is negligible compared to Of. This allows us to show the full relation 
between fa and 9i imposed by the constraint = VtcDM- We find 



fa 



2.1 X 10^3 ( Hi \ ^2/5 ^ ^ <Hi < 2.7 X 10^ GeV, (5.72) 



1.1 X 10^9 GeV [tkSt^] , Hi > 2.7 x 10^ GeV. 



fa, 



12 



: (0.23S.03)^^/(..) )'^' = (0-54 ±0.07) {eUmr'^\ fa < fa or 9, > 10-3, 
[j^j^^^^^^^^f = (7.0 ±0.8) {9ffi9.))-'/\ fa > fa or 9. < 10^3. 



(5.73) 

This function is illustrated in Fig. 15. 3[ 
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Figure 5.3. The misalignment angle 0i necessary for the axion to be 100% of the CDM in 
Scenario 11 {fa > if//27r), as a function of fa- Above the curve, ila > ^cdm- For fa > 10^^ 
GeV, one has 0i ~ IQ-^ifa/W^'^ GeY)'^/^; in particular, for fa > lO^^ GeV, the initial 
misalignment angle 9i has to assume values 6i < 10"^. 



For fa > lO^"^ GeV (m^ < IQ-^o eV or Oi < W^), one has fiOi) ~ 1 and Eq. (ISlTSD 
simplifies to 

/ f \ -3/4 

:^ (0.77 ± 0.07) x 10'=^ [ loUQ^y ) ' for > 10^^ GeV, (5.74) 

or 

9, (1.2 ± 0.1) X 10"3 ( ^Q-To"ev )'^' ' "'-^ ^ ^^-^^^ 

In particular, for > 10^^ GeV, the initial misalignment angle 9i has to assume values 9i < 
10"^ This was also noted in Refs. |60l I14j . These small values of 9i may be uncomfortable 
in a cosmological scenario. 

In the other hmit of 9i^ tt, the form of the function f{9i) assumed in Eq. (15.391) gives 

4 

TT 6 

-K — 9i^ Exp 

So, as 9i approaches vr from below, the corresponding fa approaches 0. This gives rise to 
the linear dependence of fa on Hi in the lower left corner of Fig. 15. 2[ 



(3.0 ± 0.4) X 10'^ GeV 

fa 



for /a < 2 X 10^° GeV 



(5.76) 
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5.6.1 Comparison with previous work 

Other authors find different expressions for the axion energy density. Turner |178] 
carefully analyzed the problem of the temperature-dependent axion mass, as reviewed in 
Sec. 14.2.21 cind derived the parameters in Table 14.11 We differ from Ref. |178| because we 
use a different value of the parameter a, which we computed from fitting a specific function 
of the temperature as we discussed above. For Nj = 3, we confirmed the values of the 
parameters 6, c, and d in Table WA\ 

Hertzberg et al. [88] use different values for X) which range from x = 1 (which they 
call a "moderate" value) to x = 1/20 (which they call a "conservative" value). However, in 
Ref. |88j . anharmonicities in the axion potential are not accounted for when Oi is close to vr. 
As shown in Sec. 15.2.2.41 the anharmonicity function f{6i) is essential to obtain the correct 
behavior of the isocurvature fluctuation bound at relatively small values of fa-, because the 
initial misalignment angles that give the correct axion density are not very small. We choose 
X = 1.44, but most importantly we differ from in that we account for the behavior of the 
anharmonicity function f{9i). With this prescription, it turns out that the figures shown 
in [88] are modified when the function f{9i) ceases to be of order one. This happens for 
fa < 10^^ GeV, as shown in Fig. \5M 

Sikivie [161] studies the axion energy density in the case T|**^ > A, and finds = 
0.15ef ifa/W^^GeYy/^. His expression differs from ours in Eq. (|5.3ip because of the 
different numerical factors used in Ref. |161j . namely x = li ^ = 5/12, and Tf'*^ = 
lGeV(10i2GeV//a)i/'^. 

5.7 Discussion and conclusions 

We have seen that depending on the ratio fa/Hj, two scenarios are possible for the 
axions to be 100% of the CDM, which we called Scenario I {fa < Hi/2tt) and Scenario II 
{fa > Hj/Itt). In Scenario I, one needs a cosmology in which Hf > 4.57 x 10^^ GeV, and 
for the HHCS axionic string model, one finds a PQ scale f^^^^ = (8.8 ± 1.5) x 10^° GeV, 
corresponding to an axion mass = (67 it 17) /ieV. In the DSB string model instead, one 
has /P^^ = (4.3 lb 0.7) x 10^ GeV, a value that is barely acceptable from the constraint on 
the PQ scale from astrophysical considerations, Eq. (15.67p . We remind that it is the copious 
production of axions from axionic strings in the DSB model that leads to this extremely low 
value of the PQ scale. The HHCS and the DSB string scenarios can be considered as two 
extreme models, and a realistic measurement of the PQ scale might lead to an intermediate 
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value between f^^^^ and /P^^. The measured value of the PQ scale would then yield 
information on the fraction of the total CDM axions coming from axionic strings. 

In Scenario I, the Hubble scale Hj is bounded from below to Hj > lO^^GeV (leftmost 
edge of the 0^ = ^cdm horizontal window in Fig. 15. 2p . This window can disappear 
completely if the limit from tensor modes moves to the left beyond Hi > lO^^GeV. If 
the limit on r would become more stringent than r ~ 10~^, Scenario I would have to be 
abandoned in favor of Scenario II. 

We remark that for very large values of fa in Scenario II, the initial misalignment angle 
9i has to be chosen very small for the axion to be 100% of the CDM, see Eq. (j5.73p and 
Fig. 15. 3i This may undermine the axion field as a dynamical solution to the strong CP 
problem, in that it would have to be fixed to a small value as an initial condition. 

Grand Unification Theory (GUT) models that contain axions predict that fa should be 
of the order of the GUT scale, ~ lO'^^GeV |171j . In a variety of string theory models, the 
PQ energy scale results in the range 10^^ GeV < fa < 10^^ GeV. From Fig. 15.21 we see 
that this range of fa values cannot be reconciled with axions as 100% CDM in Scenario I, 
while they can be in Scenario II provided < 10^ GeV. 



CHAPTER 6 



AXION CDM IN NONSTANDARD 
COSMOLOGIES 

Results presented in Chapter [5] greatly change when the possibility that the thermal 
history of the universe might have differed from the standard picture is taken into account. 
The aim of the present chapter is to discuss the effect of such modifications for two particular 
modified cosmological scenarios. 

6.1 Motivations in considering modified cosmologies 

The standard cosmological model has been tested up to a temperature T ~ 1 MeV, 
or down to times as short as ~ Is, when Big Bang nucleosynthesis (BBN) occurred. 
The success of the BBN theory is due to its great precision in predicting the primordial 
abundance of light elements D, ^He, and ^Li. For the success of BBN, the universe must be 
radiation-dominated at temperatures T > 4 MeV [Ml EHl ESI [90l [31] . However, due to lack 
of data prior to BBN, the history of the universe in the pre-BBN epoch T > 4 MeV is only 
indirectly inferred. 

In the standard cosmology, radiation has dominated the energy density of the universe 
before BBN since the very early time at which inflation ended. How radiation was produced 
at the end of inflation from a state of negligible temperature is still a topic of active research: 
models include the decay of the inflaton field |56t [1] and parametric resonance \55\ 11761 1107] . 

In alternative cosmological models, inflation may have ended at times close to BBN 
[21[ I112j . or there might have been a period after inflation in which the dominant energy 
density was not in radiation but in some other exotic form, like the energy density of a 
scalar field [1771 1154| . or still more there could have been an injection of entropy into the 
radiation field [5311168]. 
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6.1.1 Cosmological probes 

A good probe of the history of the universe before ~ Is is a rehc particle that has 
survived from that period. Indeed, in the standard cosmology, dark matter relics like 
Weakly Interacting Massive Particles (WIMPs) (see Refs. [Ml IMl [2H] for reviews) or like 
the axion are produced when the universe was between ~ 10ns and ~ lOps or ^ Ifis old, 
corresponding, respectively, to the age of the universe at the WIMP freeze-out for WIMP 
masses between 10 GeV and 1 TeV or at the QCD phase transition A ~ 200 MeV. These 
dark matter relics are therefore excellent candidates to test the cosmological history of the 
universe at very early epochs [19] . 

Here we show that the CDM axions may work as useful probes of the pre-BBN epoch. 
We analyze two nonstandard pre-BBN cosmologies: the low-temperature reheating (LTR) 
cosmology [Ml [Ml (Ha [111 [111 [TMl (ESI [IMl [11^ and the 

kination cosmology [591 M\ [Ml [92l [931 IIZS [Ml [US [Ml [Ml [Ml [Ml [iQl El [ZH • In the 
LTR cosmology, the expansion of the universe after inflation is driven by a massive scalar 
field <I> that decays and reheats the universe. This stage lasts down to a temperature ?rh, 
after which standard radiation-dominated cosmology applies. In the kination cosmology, 
the energy content of the universe is dominated by the kinetic energy of a scalar field 
which evolves without entropy release down to a temperature T^in, after which the standard 
radiation-dominated cosmological scenario begins. 

We find that with respect to the standard cosmology, if the Peccei-Quinn symmetry in 
the axion theory breaks during inflation, the allowed axion parameter space is enlarged in 
the LTR cosmology and restricted in the kination cosmology. Instead, if the Peccei-Quinn 
symmetry breaks after inflation, the mass of CDM axions in the LTR cosmology is orders of 
magnitude smaller than its standard-cosmology value, and it is orders of magnitude greater 
in the kination cosmology. 

Past work on axions in nonstandard cosmologies has examined the parameter space of 
hot dark matter axions in the LTR and kination cosmologies [76], and the cosmological 
bound Oq/i^ < 1 for cold (i.e. nonthermal) axions in the LTR cosmology assuming that the 
Peccei-Quinn symmetry breaks after inflation [53^ 1168^ 1113^ [T98l 1100^ I72j. 

Our work considers axions as 100% of the CDM, allows the PQ symmetry to break after 
or during inflation, and includes anharmonicities in the axion potential. 
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6.2 The misalignment mechanism in nonstandard 

cosmologies 

6.2.1 Hubble rate in nonstandard cosmologies 

We review the production of cold axions by the vacuum reahgnment mechanism, follow- 
ing the conventions in Refs. |185l I186j . The details of a cosmological axion model depend 
on the details of the cosmology before BBN only through the dependence of the Hubble 
expansion rate H[T) on the temperature T. Different pre-BBN cosmologies differ in the 
choice of H(T). This can be explained by noticing that Eq. (|5.14|) . which describes the 
evolution of the axion field in a FRW universe, depends on the specific cosmology through 
H{t) only, the initial condition 6i being unspecified and the axion mass being given by 
Eq. Km . 

We assume a scenario in which the expansion of the universe right after inflation is 
described by a nonstandard cosmology, in which the Hubble rate is H(T), until the universe 
transitions to the standard radiation-dominated cosmology. Such a transition happens 
at a characteristic critical temperature Tc which is bound by the details on primordial 
nucleosynthesis to be > 4 MeV [Ml EH ESI [90l EI] . When discussing speciflcally the LTR 
and kination models, we will refer to Tc as Trh or T^in, respectively. This said, the Hubble 
rate in this modified cosmology reads 



H{T) = F,ad(Tc) <^ 



ffB(T/Tc)\ forT<r„ 

[|^(r/r,r, forr>r,. 



where H^^(i{Tc) is the Hubble rate during the radiation-dominated epoch, Eq. (j2.87p . at 
temperature T = and v specifies the modified pre-BBN cosmology. For the LTR 
cosmology, it is = 4, and for kination cosmology it is u = 3. 

We parametrize the time-dependence of the scale factor a{t) with a generic relation, 

a{t) oc t^, (6.2) 

where /3 is a constant that depends on the details of the modified cosmology. This is a 
notation that has been already used consistently in the present thesis. For example, (3 = 2/3 
for the matter-dominated and the LTR cosmologies, (3 = 1/2 for the radiation-dominated 
universe, and (3 = 1/3 for the kination cosmology. The Hubble rate at a given time t is then 

m = ^ = (6.3) 

a{t) 

We now see how this modification affects the value of the temperature Ti at which the axion 
field starts to oscillate. 
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6.2.2 Oscillation temperature in nonstandard cosmologies 

As discussed in Sec. \5.'2.'2.'2\ coherent oscillations in the axion field begin at the temper- 
ature Ti^oiiStd defined by 

3^(^NonStd) = ma{T^°''^^'^). (6.4) 

When inserting Eq. (j6.ip in Eq (j6.4p . we obtain a relation for T^NonStd ^^lis modified 
cosmology. Since we have already discussed the outcome of Eq. (j6.4p for a standard 
radiation-dominated cosmology in Sec. 15.2.2.21 here we only consider the case T > Tc. 
Using the expression for the axion mass in Eq. (j4.29p we obtain 



jiNonStd rp ^ 



nig g*(Tc) 



l/v 



l/(w+4) 

, for T^^-^Std < 
for °nStd > ^ 



(6.5) 



6.3 Axions from string decays in nonstandard 

cosmologies 

We have discussed the effect of cosmic strings on the present abundance of CDM axions 
in Sec. 15.2.31 The main equation giving the relative abundance of string axions Q to the 
abundance of axions from the misalignment mechanism [361 [8T| [80] is expressed in Eq. ()5.40p , 

In this section, we discuss more on the parameters appearing on the right-hand side of 
Eq. (j6.6p . and extend previous theoretical results presented in Sec. 15.2.31 to obtain the 
values of the parameters in a generic nonstandard cosmology. 

Parameter Nd- In the standard cosmology, it is usually assumed that A^^ = 1, because 
for Nd > 1 a domain wall problem may arise [159j . In modified cosmologies, we take the 
same value N^, = 1 as in the standard cosmology, because N^, describes a property of the 
axion field. 

Parameter f. The value of f, which depends on the energy spectrum of the radiated 
axions |156j . was expressed in the standard radiation-dominated cosmology by Eq. (j5.4ip . 
To extend Eq. ()5.4ip to nonstandard cosmologies, we repeat its standard-cosmology deriva- 
tion in Refs. |80l [86], but we change the relation between time and Hubble parameter to 
that appropriate for a nonstandard cosmology, using Eqs. (|6.ip and (|6.3p . 
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Harari and Sikivie [86] derive the axion number density from string decays nf'^{t) from 
the equations 



dnf'{t) 1 dpr 



dt u}{t) dt 

and 



3^(t)nf(t), (6.7) 



^ - 1 - (-) 

Here, pr is the energy density of the radiated axions, oj{t) is the average energy of axions 
radiated in string decay processes at time t [80], and ps is the energy density in strings, 
given by 

.. = if^^»^ (6.9) 
From Eqs. ([SZD, dESl), and ([SS]) we obtain 

where ti is the age of the universe when its temperature is 7'^°"^*'^, and ipg <^ ti is the time 
at which the PQ phase transition occurs. The formula to obtain the parameter f fohows 
from Eq. (2.13) in Ref. [80], 

nfit) = ^^f^, (6.11) 

which, evaluated at time ti, gives 

The function uj{t) depends on the model for the energy spectrum of the emitted axions. 
For slow-oscillating strings, Davis [461 H7] argues that the energy spectrum of the axions 
radiated at time t is peaked around 2n/t, and finds uj{t) = 27r/t. Using this expression of 
u{t) in Eq. (f02D gives 

,^lj0lMti/S), for /5/ 1/3, 

\|ln(tiApQ)ln(ti/5), for/? = 1/3. 

In particular, for the standard cosmology /3 = 1/2 and f^*'^ = ln{ti/S), as in the first line of 
Eq. 

For fast-oscillating strings, Harari and Sikivie [86] argue that the energy spectrum of the 
radiated axions is broad, with a low-energy cutoff at energy -rr/ti and a high-energy cutoff 
at energy tt/S. They find uj{t) = {271 /t) \ii{t/5). Eq. (I6.12P with this expression of u}{t) leads 
to f^^'^ = 1 for the standard cosmology. Numerical simulations \36\ \8T\ [80] favor a slightly 
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smaller value of f^**^, namely 0.8 as quoted in Eq. ()5.4ip . Therefore, we decided to multiply 
Eq. (|6.12|) by 0.8. Hence, for the fast-oscillating strings. 



0.8 i^, for 1/3 

0.8| ln(ti/tpQ), for /3 = 1/3 



3/3-1' — - / (g_i4) 



In the particular case of the LTR cosmology (/? = 2/3), Eq (j6.14p gives f^"^^ = 0.27, while 
in the kination cosmology (/? = 1/3), Eq. (j6.14|) gives f^^"^ = 0.53 ln(ti/tpQ). 

In presenting our results, we use the value of f for fast-oscillating strings, Eq. (j6.14p . 
We discuss the alternative choice of Eq. (j6.13p in Sec. 16.61 

Parameter S^. The value of S^^^'^ in the standard cosmology has been discussed in 
the literature and different authors quote different results |8H \8U[ \196\ \25\ [8] . Numerical 
simulations for an evolving string network in Ref. [251 [H] yield ^^^'^ ~ 13, while simulations 
in Refs. ISU W\ lES] give f ~ 1. 

The value of ^ changes in a modified cosmological scenario. The authors in Ref. |197] 
outline a method to estimate ^ in a modified cosmology in which the universe is matter- 
dominated from the value ^^^'^ for a radiation-dominated universe. Here, we generalize the 
results of Ref. |197j to a generic cosmology with arbitrary /5, following their method. We 
define the characteristic length L for an axionic string of energy density p and tension fx 
per unit length through p = p/L?- The parameter ^ appears in the time dependence of the 
string length L as L = t/^/^. The method of Ref. |197j consists in computing ^ as 

e = (7o^t)', (6.15) 
where 70 is the fixed-point value of Eq. (14) in Ref. |197j . 

Here 7 = {H{t)L)-^, and c> is a constant determined by the value of ^ for a radiation- 
dominated universe. 

For a generic /3, we have 

w = -r 

and the fixed-point of Eq. ()6.16p follows from setting its right-hand side to zero. 



2-3/3+ V(4c + 9)/32- 12/3 + 4 
70 = . (6.18) 
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Then from Eq. (j6.15p and H = /3/t, we find 

e = ^ (2 - 3/3 + ^{ic + 9)/32 _ 12/3 + 4)^ (6.19) 

The constant c is fixed from the requirement that ^ = ^^^'^ in a radiation-dominated universe, 
/3 = 1/2. This gives 

1 + 2^^ 



— 4^std • (6-20) 

This is to be substituted in Eq. (j6.19p to find the value of ^ in the nonstandard cosmology. 
For S,^^'^ = 1, we find c = 3/4 and 



e = - (2 - 3^ + 2^/3/32 - 3/3 + ij (forftd^i). (6.21) 

Then, for the LTR cosmology, in which /3 = 2/3, we find = 16/27 = 0.5926; for 
the kination cosmology, in which /3 = 1/3, we find S,^^"" = 4(7 + 4\/3)/27 = 2.0634. 
The choice S,^^'^ = 13 is discussed in Section \6M 

Parameter In the standard cosmology, (" ~ 1 |80J. To find (ma generic cosmology, 
we use the fact that on dimensional grounds, ( is of order [8^, so that to a change 

there corresponds a change AC/2. However, the theoretical uncertainty on C is around 
50% [80], higher than the difference A^ due to the change in the cosmology used. We thus 
consider = 1 in all cosmological scenarios. 

6.4 Axion CDM in the LTR cosmology 

In the LTR cosmology [5^ flM fTTTl fTM m fTM M flM fTim fTM [721 [7^ ^ 
\57\ , the universe after inflation is dominated by a massive decaying scalar field <I> down 
to the reheating temperature Trh- The reheating temperature is defined as the temperature 
at which the decay width r$ of the scalar field <I> is equal to the Hubble expansion rate 
H{T) pl[T0Tl[T00l [i2l[72l[73l[67l[68l[66l[57]. 
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1/4 



At T > Trh, the universe follows the LTR cosmology; at T = Trhi it transitions to the usual 
radiation-dominated era. At T < Trh, the Hubble expansion rate is given by Eq. ()2.87p for 
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a radiation-dominated universe; at T > Trhi the Hubble rate H{T) depends on the scale 
factor a^^^{t) as in a matter-dominated epoch [1771 UMl ESI [Ml [Ml [IH [Ml [T00lB2]. 

/„LTR/'j' )\^/^ 

H{T) = g,ad(TRH) f ^LTR(r) ) ' for r> Trh. (6.23) 

Here -ffrad(rRH) is in Eq. (j2.87p for the radiation-dominated epoch. The relation between 
the scale factor aF^^{T) and the temperature T during the LTR epoch |113j is 

{Tf'^ T^/^ a^™(r) = constant, for T > Trh, (6.24) 

which shows that the evolution of the universe during the LTR stage is nonadiabatic because 
the entropy density in Eq. (j2.79p does not scale with (a'^™)^. Using Eq. (j6.24p and the 
usual relation during the radiation-dominated epoch, 

5,(r)i/3 Ta'*'^(r) = constant, for T < Trh, (6.25) 

we find the Hubble expansion rate 



H{T) = i/rad(TRH) <^ 



OT^fe/' for T< Trh, 
for T> Trh, 



(6.26) 



which is Eq. (16. ip in the case of the LTR cosmology v = 4, with replaced by Trh. 

In the standard cosmology, the axion field starts to oscillate at a temperature Tl^'^ given 
by Eq. (j5.22p with the standard radiation-dominated expansion rate Hts_^(T) on the left- 
hand side. In the LTR cosmology, H{T) differs from the standard expression at T > Trh, 
and the axion field may start oscillating at a different temperature TY^^ . More precisely, if 
the standard temperature Tf''^ is less than Trh, then the axion field starts to oscillate when 
the universe is radiation-dominated. Moreover, since H[T) is the same in both cosmologies 
at T < Trh, the oscillations start at the temperature tY^^ = Tf''^ given in Eq. (I5.22p . In 
this case, the results of Sec. 15.2.21 apply. 

On the other hand, if Tf*'^ would be larger than Trh, then the LTR temperature T^"™ 
will be smaller than Tf''^. In this case, the axion field starts to oscillate when the universe is 
dominated by the decay of the massive scalar field <I>. The temperature T^^R follows from 
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Eq. ()6.4p with H(T) given by the first line of Eq. ()6.ip . Since the dependence of H{T) on 
T steepens from to as T becomes greater than Trh, it follows that 

pLTR ^ pstd_ ^g_27) 

Introducing the notation 

g*{TRH) 



5rh(T) = ^1^. (6.28) 



and with Trh.McV = ^rh/McV, we find 

1/8 



am^MpiTl^A^QC^ / J^, ') , for T^tR > Aqcd, 
maM„r2Hy=I=:j , for TLTR < ^Q^^. 

Numerically, 

LTR^ f(160±2) MeV5^r(rr^)<Mev/Mf , for T^™ > Aqcd, 

\(350±8)MeV5;;f (Tr^)r^gMev/a:if> for T^™ < Aqcd- ^ 

To summarize, if Trh < Aqcd, 



jnLTR 



' (67 ± 3) GeV <7*"'^'(rf , for Tf ^ < Trh, 

(350 ±8) MeY 9-ll'{Tl^^^)T^'f-ll\ for Trh < T^™ < Aqcd, (6-31) 
I (160 ± 2) MeV5;^l/^'(rL™) T^A /"V^, for Aqcd < T^™; 



if Trh > Aqcd, 

jnLTR 



f (67 ± 3) GeV5*"'/'(Tf d)/-,^^ for Tf ^ < Aqcd, 

(607 ± 27) MeV5*"'/''(Tf d)/-^f , for Aqcd < Tf" < Trh, (6.32) 

t(160±2) MeYg-f'{Tl^^^)T'^'f-li\ for Trh < T^. 



Also, the present axion energy density is modified from the standard case if Tl^'^ > Trh. 
We examine the misalignment mechanism and the production in string decays separately. 

String decays in the HHCS give the contribution to the present axion energy density 
expressed in Eq. (|6.6|) . 

j^LTR,str ^ gLTR j^LTR,mis ^ q -^g f^LTR,mis^ (g 33) 

where we have used the values = 1, f^'^^ = 0.27, S,^"^^ = 16/27, and C = 1 in Eq. ^^6^, 
consistent with the discussion in Sec. 
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For the misalignment mechanism, the axion number density at the present time can be 
found from considering the expression for n™'^(Ti) in Eq. (j5.25p . with Ti = TY^^, then use 
the conservation of axion number in a comoving volume, na{T) oc a~^{T). This gives 



[n.(T-K) forT->TKH. 

Here na(Ti) is the function given in Eq. (j5.25p . One clearly has 

nL™(ro) = <*'i(To) forrfd<TRH. (6.35) 

In the case Tf^'^ < Trh coherent oscillations begin in a radiation-dominated universe, 
and the number density of cold axions is the same as in the standard cosmology. Instead, 
for T|*'^ > TrH) one obtains a different axion density due to both the release of entropy and 
to the different function for the scale factor a^'^^{T) in this modified cosmology. To better 
understand the origin of the difference, it is convenient to introduce the ratio between the 
present density n^"'"^(To) in the LTR cosmology, and the present density ?i^*'^(To) that would 
ensue if the cosmology were standard at temperatures T > Trh- We write, for Tf^'^ > Trh, 

nL™(ro) _ AtLTR yUTR 

where 



.34) 



.36) 



m^'^ na(rf d) V a^*^(rf d) J ^^-^^^ 

is the standard-cosmology ratio of the comoving number of axions N^"^^ at the temperature 
T^'^^ to the comoving number of axions N^^^ at the temperature Tf^'^, and 

yLTR /„LTR(^LTR)X3 

- ' ^ ^ ' ^ (6.38) 



is the ratio of the LTR-cosmology volume V^"^^ to the standard-cosmology volume V^^'^ at 
the temperature T^"^^. 

The ratio N^'^^/N^^'^ accounts for the fact that coherent oscillations in the axion field 
start at a different temperature in the LTR cosmology compared to the standard cosmology. 
The ratio yLTR^ystd accounts for the fact that at temperature T^"^^ the scale factors, and 
so the volumes, in the LTR and in the standard cosmologies differ due to entropy production 
from the decay of the scalar field in the LTR case. 
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Using the relations between temperature and scale factor during the radiation and LTR 
epochs, Eqs. ()2.8ip and ()6.24p respectively, we find 

iv"^^_ 9*s{Tn ( Tr y 

and 

(6.40) 



V-''' <7*5(Trh) 5*2(TL™) VtL™; • 
Prom T^LTR < 2^std (g^g gq. (fOTD l. we find that N^^^ > N""^"^. But for T^^^ > Trh, 
yLTR ^ ystd^ The latter factor dominates, and ?i]^'^^(To) is less than n^*'^(ro). 

The present axion energy density from the misalignment mechanism follows as, in units 
of the critical density. 



,LTR,mis _ ^anLTR(ro) _ / 0^ '"^^ for < Trh, 

" Per: -\nr'-^^, forrr>TKH. 



Here Qa*'^'™^'^ is given in Eq. (j5.55p . 

The first line of Eq. (|6.4ip . valid for Tf^'^ < Trh, is numerically equal to Eq. ()5.55p . The 
second hne of Eq. (lOTT) . valid for Tf^ > Trh, is 

^LTR,.. , ^ V^)''' ^*(^o)9^riTF^) .Hv^gg, (^^m)^ for/. < /.(Trh), 

¥9*(ro)^g||^(e?m)>, for/.>A(TRH). 

(6.42) 



a 



Numerically, 

j^LTR.mis;^: 



LTR,mis,2 _ 1 (2-0 ± 0.2) X lO^^ {9f f^Q,)) g^^^T^^'') fj^ TIumcV^ for < /.(Trh), 
\(7.5 ± 0.8) X 10-8 {9f f{9i)) /2,2 TRH,MeV, for/. > /.(Trh), 

(6.43) 



where 

/a(TRH) = (5 ± 1) X lOi^ GeV9^^/^TL™)T2H,MeV, (6-44) 

is the PQ scale at which the two lines in Eq. (j6.43p match. 

The present axion energy density in the LTR cosmology, for Scenario I, is given by the 
sum of the misalignment mechanism and the string decay contributions 

OLTR _ oLTR,mis , oLTR.str _ J ' + Q^^'^)j fo^^ ^f*'^ < Trh, , . 

^ " " ^ ' "l^^r-'^^^d + Q^™), for Tf^> Trh. ^'-''^ 
Here, Q"*'^ and Q^TR ^j^j^^g ^.j^g ^.^^^^^ p^^' (Tq) / p]^'"" {Tq) in Eq. (ICTD in the 

standard and LTR cosmologies, respectively. In the HHCS model, we obtained 

Q"*'^ = 0.8, and Q^™ = 0.16. (6.46) 
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Numerically, the first line of Eq. (j6.45p corresponds to the total (strings plus misalignment) 
energy density in axions found in the standard cosmology, Eq. (|5.57p . while the second line 
is valid for Tf^'^ > Trh and gives 

^LTR^2 ^ f (2.3 ± 0.2) X 10-6 {9f /(aO) ^RH^'mLTR) /fg r|H,MeV> ^or < /.(Trh), 
\ (8.7 ± 0.9) X 10-8 {ef f{e,)) fl,2 TRH.MeV, for fa > UTrr), 

(6.47) 

This equation, where axionic strings contribute to the value of ^2^"'"^, is valid only in Scenario 
I {fa < H 1/271); otherwise, Eq. (j6.43|) suffices and yields the total axion energy density in 
the LTR cosmology. 

6.4.1 Results for the LTR cosmology 

We now derive the regions of axion parameter space where the axion is 100% of the 
CDM in the LTR cosmology. We then compare them to the standard-cosmology regions. 

The CDM axion parameter space in the standard cosmology depends on the PQ energy 
scale fa (or alternatively the axion mass nia), the initial misalignment angle 6i, and the 
Hubble parameter during inflation Hj. In the LTR cosmology, an additional parameter is 
included, the reheating temperature Trh- 

If the PQ symmetry breaks after the end of inflation (Scenario I, fa < Hj/2'k), there 
is only one PQ scale fa for which the totality of the CDM is made of axions. There 
correspondingly is also a single value for the axion mass m^. In the standard cosmology, we 
referred to this scale as fa^'^, whose value is given in Eq. (j5.63p for the HHCS axionic string 
model. Here, we repeat the same computation in the LTR cosmology: we use the equate 
the value for Ocdm h"^ in Eq. (j2.90p with the expression for Oq"""^ /i^ in Eq. (|6.47p setting 
{Of fiOi)) = 1.4 7rV3, obtaining the PQ scale 

= (4.7 ± 0.6) X 10^^ GeV5i'S(TL™) T-^^^^, (6.48) 



or 



ma 



LTR 



(13 ± 2) neV <7]^f (rr^) T^gMeV (6-49) 



In Fig. 16.11 we plot fa"^^ as a function of Trh- The jumps and kinks in the fa"^^ line 
are due to the different values of 5*(Trh) and g^iTj""^^) in Eq. (I2.86p . There is also a 
(visually small) discontinuity between the fa^^ and f^^'^ lines at Trh = ^i*'^ = T^^^ due 
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Trh [MeV] 



Figure 6.1. The Peccei-Quinn scale ^ ^ function of the reheating temperature Trh 

for the axion to be 100% of the CDM in Scenario I {fa < Hi/2tt). Also shown are the PQ 
scale Z^**^ in the standard cosmology, and various constraints (shaded regions). 



to different contributions from string decays. In fact, from the first line of Eq. (|5.55p and 
from Eq. (j6.45p we have 

LTR _ / (0.23 ± 0.03) {ef f{e,)) fI% (l + Q-'% for Tf d < Trh, 

- 1(0.23 ± 0.03) (^.^m))i^S^^(l + Q^™), forTf'i>TKH. ^ ^ ^ 

Equating the two lines in Eq. (|U3U|) at Trh = TT*"^ = TY^^, where A^^"^^ = N''^'^ and 

/ 1 , r)std \ 6/7 



Numerically, /a'^^(rRH = rf'i) = 1.45 /f'^, which is slightly higher than /f^. 

In Fig. 16.11 we also shade out the following bounds: the bound from white dwarfs 
cooling times in Eq. (]5.67p : the indirect bound on the PQ scale from the nondetection 
of primordial gravitational waves arising from fa < Hi/2tt and Eq. (|2.65p (region labeled 
"Tensor Modes"); the bound on Trh > 4 MeV from Big Bang Nucleosynthesis; and the 
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bound from the ADMX experiment \12\ [T3] excluding a KSVZ axion with a mass ma 
between 1.9 /xeV and 3.53 /xeV, see Eq. (j4.6ip . The dashed hne marks the requirement 
that the axion starts to oscillate in the LTR cosmology, Trh < Tf*^, with Tf^'^ given by 
Eq. (j5.22p . The ADMX bound can be rephrased as an exclusion bound for the reheating 
temperature Trh- Using the expression for the axion mass in Scenario I, Eq. ()6.49p . the 
ADMX result corresponds to an exclusion of the region 85 MeV < Trh < 140 MeV, valid 
for KSVZ axions. 

Depending on Trhj fa'^^ may differ from /^^^^ in Eq. (j5.63p by orders of magnitude. 
The maximum value of is achieved for Trh = 4 MeV and is, with (?*(Trh) = 10.75, 

/a™ — (2.0 lb 0.3) X lO^'* GeV. This value is three orders of magnitude larger than /^^^^ in 
Eq. ()5.63p . As discussed in Sec. 16.61 these large values of the Peccei-Quinn scale correspond 
to axion masses that are beyond the reach of current DM axion search experiments. 

In Scenario II {fa > Hj/2-k), the parameter space is bounded by the nondetection of 
axion isocurvature fluctuations in the CMB spectrum, see Eq. (|5.48p . For Trh > Tf^'^, 
the isocurvature bound has the same expression as in the standard cosmology, namely 
Eq. (15:72]) . For Trh < Tf^, we eliminate Oi in Eq. by equating I^cdm /i^ with the 

expression for f^^™ derived previously in Eq. (I6.43p . with no contribution from axionic 
strings. The LTR isocurvature bound is then, for Trh < Tl^'^ and (7*(T|'*^) = 61.75, 



Hi. 



12 



'2.0 X 10-2 [/(0,)]-i/2 ^1/4 ^ 8 2 X 1014 GeV < fa < faiTnn), 

5.0 X 10-2 T-^%^y, for /.(Trh) < fa- 

(6.52) 

Here ^(Trh) = (5 ± 1) x lO^^ GeV g"]/' (Tf™) T^^^MeV is given by Eq. ^M- The LTR 
isocurvature bound can be approximated by 

■ 1.3 X 10-4 fa,i2, for fa < 8.2 X 10^4 T-^/^^v GeV, 

2.0 X 10-2 T-^ fl% for 8.2 x lO^^ T'^/^ GeV < < ^(Trh), (6.53) 



Hi. 



12 



RH,MeVJa,12' j-w^ u.^ /n -^RH.MeV 
5.0 X 10-2 r-j^/]^^v. for fa > /a(TRH) , 



As in the case of the standard cosmology, there are two changes in the power-law dependence 
of Hi 12 on fa^i2 in the LTR cosmology, the first one being at fa = 8.2 x lO^^ Tj^^'mcV GeV 
and the second one at fa = /a(rRH)- Notice that at large fa, the isocurvature bound is 
independent of fa- 

happens for 



When Trh = Tf^'^, the LTR and the standard isocurvature bounds coincide. This 



/, = 4.6 X 1021 Tp-2 ,^^v GeV, for Tf^ < Aqcd, ,^ 

fa = 5.3 X 1028 T^l Qgy^ fo^ j^std > J^^^^_ 
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Fig. l6.2l shows the regions of the parameter space (/a, Hj) where the axion is 100% of the 
CDM in the LTR cosmology. The axion mass scale on the right is Eq. (j4.12p with = 1. 
The region labeled "Tensor Modes" is excluded by the nonobservation of tensor modes in 
the CMB fluctuations, Eq. ()2.65p . The region labeled "White Dwarfs Cooling Time" is 
excluded from astrophysical observations of white dwarfs cooling times for KSVZ axions, 
Eq. (I5.67P |147] . The line fa = Hi/2tt divides the region where the PQ symmetry breaks 
after inflation (Scenario I, fa < Hj/2tt) from the region where it breaks during inflation 
(Scenario II, fa > Hi/2Tr). In the lower right region (Scenario I), the axion is the CDM 
particle if fa equals the value fa"^^ given by Eq. (j6.48p . In the upper left region (Scenario II), 
we plot the isocurvature bounds to the allowed parameter space for the standard cosmology 
(thick line) and for Trh = 4 MeV (dotted hue), 15 MeV (dot-dashed line), and 150 MeV 
(dashed line). For a given Trhj the isocurvature bound with the LTR cosmology lies below 
the standard line, because the entropy dilution term ~ (Trh/^i)^ lowers the axion energy 
density. Thus, the LTR cosmology allows more CDM axion parameter space than the 
standard scenario. 




Figure 6.2. In the LTR cosmology, axions are 100% of the CDM in the white region on the 
left (limited by a different line for each Trh) and in the narrow bands marked by horizontal 
lines in the lower right triangle (one line for each Trh)- 
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In the allowed region of parameter space for Scenario II, the axion can be 100% of the 
CDM provided the value of 9i is chosen appropriately. In the standard cosmology, 6i is a 
function of fa only [185^ [8^. see Eq. ()5.73p . Instead, in the LTR cosmology, 6i depends on 
both fa and Trr. Taking 5*(rfTR) = 61.75 we find 

2/3 



fa, 



12 



\' ' -5/2 

g2 f(e.)r^ , for fa < /a(rRH) Or Oi > ITTpH MeV 

i J \ ' RH,McV / 

for fa > /a(TRH) or Oi < 17 T^^/' 



(6.55) 



The relation in Eq. ()6.55p is plotted in Fig. 16.31 for the standard cosmology (thick line) and 
for Trh = 4 MeV (dotted line), 15 MeV (dot-dashed line), and 150 MeV (dashed line). 
As Trh decreases, one departs from the standard cosmology. The value of Oi at fixed fa 
increases when Trh decreases. 



rria [eV] 



IQ-i IQ-^ IQ-^ IQ-^ IQ-^ 10"^ IQ-^ 10"^ IQ-'^IG- 




fa [GeV] 



Figure 6.3. The initial misalignment angle 9i as a function of the Peccei-Quinn scale fa for 
the axion to be 100% of the CDM in Scenario II {fa > H 1/271): standard cosmology (black 
solid line), LTR cosmology with Trh = 4MeV (red dotted line), 15MeV (green dot-dashed 
line) or ISOMeV (blue dashed line). 
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6.5 Axion CDM in the kination cosmology 

We now discuss axion CDM in the kination cosmology pi Ml [Ml ESI EHl ^HM ^iM 
[1461 [UHl \IM [Uni [MI Uni HH [H]. The Hubble parameter for this pre-BBN cosmology is 
(see Ref. I39]) 



H{T) = H,^^{T^^) { V ^^x'-^J V^kiny ^ (6.56) 



9. (Tkin) V^ki: 



forT > Tkin- 



Here Tkin is the temperature at which the universe transitions from kination domination 
to radiation domination. Entropy is conserved during the kination cosmology, so the scale 
factor during kination a'*™(T) follows the same temperature dependence as a^^'^{T), 

alsi'^) T a^^{T) = constant. (6.57) 

For this reason, here we do not make a distinction between the scale factor in the standard 
and in the kination cosmologies. We write a'*™(T) = a^^'^{T) = a{T). 

If Tf''^ < Tkin, coherent oscillations in the axion field start in the radiation-dominated 
universe, at the temperature Tf^'^ given in Eq. (j5.22p . On the contrary, if Tf''^ > Tkin, 
coherent oscillations in the axion field start when the universe is in its kination stage. In 
this case, the temperature T}^™ at which axion oscillations begin is given by the following 
expression: 

ykin ^ f 2.48 GeYg-^^'{T^-) T^^%^^ /"^f , for T^- < Aqcd, 
' \331 MeY g-^/'\T^^) T^/,^^,^ Cil' ^ for Tf- > Aqcd- 

Here Tkin,McV = ^kin/MeV. 

The axion energy density in the kination cosmology has contributions from string decays 
and from the misalignment mechanism. 

String decays give a contribution to the present axion energy density 

j^kin,str ^ gkin j^kin,mis ^ q 43 f^"" (6.59) 

where to compute Q^^^ in Eq. (|5.4U|) we used A'^^ = 1, = 2.06, C = 4.9, and f^^^ is given 
in Eq. (KIM as 

In the last expression, we used the relation t oc 1/H{T), the fact that at the time of the 
PQ transition tpQ, the temperature of the universe is T = /„, and the fact that at the time 
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ikin 



ti, the corresponding temperature is T^. Using the expression for the kination Hubble 
parameter in Eq. (|6.56p and neglecting the term ln(A / ^'"'l^j'i^L ) ~ 1, we obtain 



<?kin{riK'") 

= 1.6 Hfa/T^'^). (6.61) 

The temperature T}^™ is greater than Aqcd for any value of Trh and any value of fa for 
which there are contributions from string decays (Scenario I, fa < Hj/2'ir). Thus, using the 
expression for T^^™ in the second line of Eq. (j6.58p . we obtain 

fl^i- = 57 + y In fa,12 - InTkin.McV- (6.62) 

In the region of the parameters of interest for kination, f^^"^ ~ 35. Thus, axions from 
strings dominate the total axion population, the energy density JIq™''^*'^ being one order of 
magnitude larger than ila™'™'^. We notice that this is opposite to what we obtained in 
the standard and LTR cosmologies, where the radiation of axions from axionic strings is a 
subdominant production mechanism for cold axions. 

The contribution from the misalignment mechanism results from the conservation of the 
axion number in a comoving volume, na{T) oc a~^(T). This gives 

n^(ro) 



^-(^f ') ( w) ' for Tf d < Tkin, 



(6.63) 



Here na(Ti) is the function given in Eq. (|5.25p . One clearly has 

n^-(ro) = (To) for Tf^ < Tkm- (6.64) 

For Tf <^ > Tk in, one obtains a different axion density. As for the LTR cosmology, we 
introduce the ratio between the present density nJ^™(To) in the kination cosmology, and the 
present density n^**^(To) = /3^*'^(To)/ma, see Eq. ()5.29p . that would ensue if the cosmology 
were standard at temperatures T > Tkin- We write, for Tf^'^ > Tkm, 

n«td(To) ATstdystd' 

where iykin/^std ykin/ystd defined as follows. The ratio iykin/^std ^^le standard- 
cosmology ratio of the comoving number of axions N^^^ at the temperature T^™ to the 
comoving number of axions A^'^*'^ at the temperature Tl^'^. Using Eq. (j6.57p . we write it as 

ATk- _ na{T^n ^a(^k"^)^^' _ uaiT^n gMTf") /^fd^^=' 



^std „^(j.std) \^„(2.std)y n,(Tf d) 5,5(rki") VTk- 
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The ratio ykmyystd -g ^j^^ ratio of the kination-cosmology volume V"^™ to the standard- 

'kin 
1 ' 



cosmology volume V^^"^ at the temperature Tf^^ 



17^= („std(2^kin)J =1' (6-67) 

The last equality follows because no significant entropy is released during the kination stage 
[T50] . so aki'i(r) = d'^^iT). 

The present axion energy density from the misalignment mechanism, in units of the 
critical density, is therefore 

CQStd,mis fnrT'^'^^T, ■ 

Qkin.mis _ J ' iUI ^ -tkin, /„ 

Inserting numerical values, the first line of Eq. (j6.68p is given by Eq. ()5.31|) . while the 
second line reads 

j^kin,mis;^2 ^ HgO ^2 ^^2 f^^^^ T^^^^y. (6.69) 

Due to the peculiar dependence of the Hubble rate with temperature in kination, H(T) ~ 
T^, there is no distinction in Eq. (j6.69p between Oa'"'™"^ for Tj^'" > Aqcd and for T^^'" < 
^QCD- 

Finally, the present axion energy density in the kination cosmology is given by the sum 
of the misalignment mechanism and the string decay contributions 

{QStd,mis ^std\ f T^std ^ J- . 

Here, Q"*'^ and Q^'"" are the values of the ratio p''^'^ (Tq) / p]^''%To) in Eq. ([S^iO]) in the standard 
and kination cosmologies, respectively. 

6.5.1 Results for kination 

We now derive the regions of the axion parameter space where the axion is 100% of 
the CDM in the kination cosmology. We then compare them to the standard-cosmology 
regions. 

The axion parameter space in kination cosmology depends on /„, Hj, 9i, and the 
additional parameter Tkm- 

If the PQ symmetry breaks after the end of inflation (Scenario I, fa < Hj/2tt), there 
is only one PQ scale fa for which the totality of the CDM is made of axions. There 



97 

correspondingly is also a single value of the axion mass ma- In the kination cosmology, 
using the observed value of r^cDM^^ in Eq- (I2.90P . and the expressions for Vl)^^ derived in 
this section, we find 

ff^ = (7.9 ± 0.2) X 10^ Ge\gl'\T^r.) ^ le . ^^^^'^ ^ > (6-71) 

57 + In f^J^ - In Tkin^MeV 

and 

m^- = 739 ± 22meV3,-^/^Tkin) t^'^T^Im.w (6-72) 

In Eq. ()6.7ip . we used the explicit expression for f'^™ derived in Sec. 15.2.31 

In Fig. 16.41 we plot as a function of Tkin- The function does not present jumps, 

because both 5'*(r^™) and g*{T\^h\) do not change in the domain of /a™- 

The discontinuity between the /^'" and f^'^ lines at T^in = Tf*^ = is due to different 

contributions from string decays. In fact, from Eqs. (j5.3ip and (j6.70p we have 
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Figure 6.4. The Peccei-Quinn scale as a function of the kination temperature T^in 
for the axion to be 100% of the CDM in Scenario I (/a < Hi/2tt). Also shown are the PQ 
scale f^'^ in the standard cosmology, and various constraints (shaded regions). 
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n 



ikin 



'a 



i.32g;^/^2(rf d) {ef fie,)) {f^^y/^ (i + q^'^) 



for T\ 



<rki 



(6.73) 



Equating the two lines in Eq. at Tkin = Tf^ = T^'"", where N^''' = N''^'^, we obtain 



In Fig. 16.41 we also shade out the following bounds: the bound from white dwarfs cooling 
times in Eq. (j4.43p : the indirect bound on fa from the nondetection of primordial gravita- 
tional waves arising from fa < Hi/27r and Eq. ()2.65p (region labeled "Tensor Modes"); the 
bound on Tkin from Big Bang Nucleosynthesis; and the bound from the ADMX experiment 
excluding a KSVZ axion with a mass ttt-q between 1.9 /zeV and 3.3 /leV. The dashed 
line marks the requirement that the axion starts to oscillate in the kination cosmology, 
Tkin < Tl^"^, with Tfi given by Eq. (15:221) . 

The PQ scale /^'° is orders of magnitude lower than the PQ scale fa^'^ in the standard 
cosmology. The low values of in comparison with fa^'^ is due to two different reasons. 
The first reason is that, since coherent oscillations of the axion field start later in the 
kination cosmology than in the standard cosmology, the initial comoving number of axions 
^ykm -g j^igj^gj- than iV*^**^. The second reason is that the contribution from axionic strings 
to fi^™ in the kination cosmology is much higher than the same contribution to O^*'^ in the 
standard cosmology. Then, at a given PQ scale fa, the energy density J^^" > Ql^"^. A lower 
PQ scale is thus required in order to have the same CDM energy density J^cdm- 

The PQ scale can be so small as to violate the limit from the white dwarfs cooling 
time in Eq. ()4.43p . This imposes the requirement Tkin > 217 MeV if axions are 100% of the 
CDM. This requirement is more stringent than the BBN constraint Tkin > 4 MeV. 

In Scenario II (/„ > Hf/ln), the parameter space is bounded by the nondetection of 
axion isocurvature fiuctuations in the CMB spectrum, Eq. (I5.48p . For Tkin > Tf^'^, the 
isocurvature bound has the same expression, Eq. (15.72p . as in the standard cosmology. For 
Tkin < T|**^, we eliminate 6i in Eq. (|5.48p by equating Ocdm with the expression for ilj^™ 
derived in this section. The resulting kination isocurvature bound for Tkin < Tf ^ is 




(6.74) 



We find /^''^(Tkin = Tfd) = 2.04 x 10^ GeV. 



< 7.48 X 10 




(6.75) 



This bound can be approximated by 
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^ ^ fl.31 X 10-Va,i2, for/, < 3.26 xlO^GeVTkin.MeV, 

'''' \7.48 X lO-VaS^kYn'MeV ^ > 3-26 X 10^ GeVTkin.MeV- 

Contrary to the cases of standard and LTR cosmologies, in the kination cosmofogy, there 
is only one change in the power-law dependence of Hj 12 on /a^i2, namely at = 3.26 x 
10^ GeVTkin,McV- This change is due to the effects of anharmonicities. 

When Tkin = Tf'**^, the kination and the standard isocurvature bounds coincide. This 
happens for 

f/„ = 4.6 X 1021 r-2^^^GeV, for Tfd < Aqcd, 
\/a = 5.3 X 1028 TiTjj^^^GeV, for Tf ^ > Aqcd- 
Fig. 16.51 shows the regions of the parameter space where the axion is 100% of 

the CDM in the kination cosmology. The axion mass scale on the right is Eq. (j4.12p with 
= 1. The region labeled "Tensor Modes" is excluded by the nonobservation of tensor 
modes in the CMB fluctuations, Eq. (j2.65p . The region labeled "White Dwarfs Gooling 
Time" is excluded from the bound in Eq. (|4.43|) . The line fa = Hj/2tt divides the region 




Figure 6.5. In the kination cosmology, axions are 100% of the CDM in the white region 
on the left (limited by a different line for each Tkin) and in the narrow bands marked by 
horizontal lines in the lower right triangle (one line for each Tkin). 



100 

where the PQ symmetry breaks after inflation (Scenario I) from the region where it breaks 
during inflation (Scenario II). In Scenario I, the axion is the CDM particle if fa equals 
the value /^'" in Eq. (1677X1) . In Scenario II, we plot the isocurvature bounds to the allowed 
parameter space for the standard cosmology (thick line) and for T^in = 4 MeV (dotted line) , 
300 MeV (dot-dashed line), and 700 MeV (dashed line). In the allowed region of parameter 
space for Scenario II, the axion is 100% CDM if 

„ _ ^^CDM^i^ ^^^(^km) 2kin,MeV '7Q\ 

TiEoefm • ^ ' 

We plot this relation between fa and 9i in Fig. 16.61 for the standard cosmology (thick line) 
and for Tki^ = 4 MeV (dotted line), 300 MeV (dot-dashed line) and 700 MeV (dashed line). 
As Tkin decreases, one departs from the standard cosmology. The value of 9i at flxed fa, or 
of fa at flxed 9i, decreases when Tkin decreases. 




Figure 6.6. The initial misalignment angle 9i as a function of the Peccei-Quinn scale fa 
for the axion to be 100% of the CDM in Scenario II (fa > Hi/2tt): standard cosmology 
(black solid line), kination cosmology with T^in = 4MeV (red dotted line), 300MeV (green 
dot-dashed line), or 700MeV (blue dashed line). 



101 

6.6 Discussion 
6.6.1 Comparison to previous work 

Axions in kination cosmology were studied only in Ref. [7^ and only as hot dark matter 
(i.e. thermally produced in the hot primordial soup). To the extent of our knowledge, CDM 
axions in kination cosmology were not examined before. 

Axions in the LTR cosmology were studied before [531 [Ml IIZZI (IMl ESI [Ml [TOTl 
11001 172] . but only to determine the cosmological bound on the PQ scale in what we call 
Scenario I, namely fa < Hi/2tt, in which the Peccei-Quinn symmetry breaks after the end 
of inflation. Our work can also be used to set an upper bound on the PQ scale by imposing 

n^h^ <ncDMh\ (6.79) 

These bounds can be read off the figures and the equations we presented previously, all 
of which represent the equation Q.ah?' < Q.cDuh'^- We remark that therefore, our work 
extends previous papers in that we have examined also the region > i?//27r, where the 
PQ symmetry breaks during inflation, have updated r^cDM to the current observational 
value, have used an improved constraint on Trh from Big Bang Nucleosynthesis, and have 
included anharmonicities in the axion potential. 

Our numerical result for the highest allowed value of the PQ scale in the LTR cosmology 
/LTR = 8 58 X 10^3 GeV, obtained for Trh = 4 MeV, differs from previous authors. 
Steinhardt and Turner |168j showed that the entropy production due to the decay of a 
massive scalar field raises the maximum PQ scale to ~ 10^^ GeV, but they were 

corrected by Kawasaki et al. [101^ llOOj for using the value of Ti in the standard cosmology 
instead of the LTR cosmology. Kawasaki et al. [lOllllOO] used Ti in the LTR cosmology and 
obtained /a"""^ ~ 10^^ GeV; however, they used VLah? = 1 and Trh = 1 MeV instead of the 
current value ^cnuh^ = 0.1131 ± 0.0034 and the current BBN bound Trh > 4 MeV. 
Giudice et al. [72] find for the maximum PQ scale in the LTR cosmology the value 
/a™ ~ 10^^ GeV, which is higher than ours for the same reasons as for Kawasaki et 
al. 

6.6.2 Effects of changing the string decay parameters 

The computation of /f^ in Eq. dSJQj), fa^^ in Eq. ([HliH]) . and /^'" in Eq. (IHTTl) strongly 
relies on the model used to describe the axionic string evolution and the energy spectrum 
of emitted axions. 
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In Sec. I5.2.3|, we discussed the dependence of f and ^ on the model for the axionic 
string oscillation and radiation spectrum. There, we showed how these quantities in a 
modified cosmological scenario are related to their values in the standard cosmology. For 
the latter, we used the values f'^^'^ = 0.8 and = 1, obtained assuming that an axionic 
string radiates axions in a broad energy spectrum and that the axionic string network is 
a global string network [1591 [36| 11961 \8T\ [80] . We have referred to this set of assumptions 
as the HHCS model. With the values of the axionic string parameters set by the HHCS 
model, the contribution from axionic string to the total axion energy density in the standard 
and LTR cosmologies is subdominant compared to the contribution from the misalignment 
mechanism, while it is dominant in the kination cosmology. 

We now discuss the modification to the axion parameter space when we assume that 
axionic strings radiate axions in a narrow energy spectrum [46[ Wl\ HSj [23| [22| I156j , and 
that the axionic strings network is a local strings network |25[ [H]. In this case, f^^'^ = 70 
and ^'^^'^ = 13. We call this set of assumptions the DSB model. With these values, the 
contribution from strings to the axion energy density in the standard cosmology is dominant 
{Nd = 1, C = 4.9), 

j^std,str_ 200O^*<^'™^ (6.80) 

This affects the value of fi^*'^ = ^2^*'^''**'^ + ila''^'™"*, which in the DSB model is about two 
hundred times higher than in the HHCS model. As a consequence, within the DSB model, 
the value of the PQ scale f^^'^ for which Q^^'^ = Ocdm in Scenario I is 

This value is smaller than that in Eq. ()5.70p . obtained within the HHCS model. It is of the 
order of the astrophysical constraint from white dwarfs cooling times in Eq. (j5.67p . 

We conclude that, depending on the model for the axionic string and its emission 
spectrum, f^^'^ in the standard cosmology can range from the value in Eq. (j5.70p to the 
astrophysical bound from white dwarfs cooling times 4 x 10^ GeV. Correspondingly, m^'^ 
can range from 67 it 17 /ueV to 15 it 4 meV. 

In nonstandard cosmologies, when going to the DSB model, we must redo the calculation 
of f and ^ using the formulas in Sec. For the DSB model value f = 13, Eq. (IMjl 

gives c = (\/52 + l)/52 ~ 0.158 and Eq. (f6T8]l gives 

e ~ 10 (2 - 3/3 + 2^2.41^2 _ 3^ + (for f *<i = 13). (6.82) 
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The values of in the DSB model are then ^^"""^ = 2.8 for the LTR cosmology (/3 = 2/3)and 
^kin _ £qj, ^YiQ kination cosmology (/3 = 1/3). 

For the parameter f, we turn to Eq. (|6.12|) . The results from the HHCS model favor 
the fast-oscillating axionic strings model, which predicts f as given in Eq. (j6.14p . The DSB 
model points toward a slow-oscillating axionic string for which f is given by Eq. (j6.13p . 
Within the DSB model, for the illustrative case 6 = (10^^ GeV)~^ and Trh = ^kin = 
4MeV, we obtain f ~ 70 in the standard cosmology, f'^™ ~ 20 in the LTR cosmology, and 
^kin ^ 4300 in the kination cosmology. 

The LTR and kination axion energy densities from axionic strings in the DSB model are 
then, with Nd = 1 and ( = 4.9, 

j^LTR,str^2 ^ ^-^^ j^LTR.mis ^^^^ J^km.«f/i2 = 3.6 x 10^ (6.83) 

The higher axionic string contributions in the DSB model with respect to the HHCS model 
sensibly lower the values of the PQ scales fa"^^ and f^^"^ for which the axion is 100% of the 
CDM. We have, taking Trh = Tkin = 4MeV, 

/LTR = (1,43 ± 0.04) X 10^3 GeV and /^'" = (4.3 ± 0.1) x 10^ GeV. (6.84) 

6.6.3 Distinguishing nonstandard cosmologies observationally 

One might be able to distinguish different nonstandard cosmologies before BBN by using 
properties of the axion CDM particle. 

Future CMB measurements in the tensor modes spectrum from the PLANCK satellite 
may help with this task in the following way. Suppose that the cosmology is standard and 
the axion is 100% CDM. A detection of tensor modes within the PLANCK range would 
yield the value of Hj, and thus that of ma |1851 184] . If the axion mass is found to be 
different from the one derived from the PLANCK measurement, then either the axion is 
not the CDM particle and the cosmology is standard or the axion is the CDM particle and 
the cosmology is nonstandard. 

PLANCK will also improve the current bounds on axion isocurvature fluctuations, or 
detect them. If the axion is the CDM particle, then a combination of the measurements 
of the axion mass and of the axion isocurvature fluctuations yields a point in the Hj — fa 
plane, which might be excluded in the standard cosmological scenario but not in the LTR 
scenario. Combining the axion searches and the PLANCK measurements thus might lead 
to new information on the pre-BBN history of the universe. 
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Outside CMB measurements, one may try to distinguish nonstandard cosmologies by 
measuring both the axion CDM density J^cdm and the axion mass nia- However, one 
immediately runs into the following problem. 

Assume, for example, that the axion is found to be the main CDM component and the 
axion mass is measured at ttIq ~ 10^'^ eV. These facts can be ascribed to two different 
cosmological models. The first model involves the axion field evolving in the standard 
cosmology, with the dominant contribution to the total axion energy density coming from 
axionic strings and only a tiny fraction from the misalignment mechanism, as the DSB 
model would predict. The second model involves a stage of kination before BBN lasting 
until Tkin ~ 900 MeV, with the contribution from axionic strings and from the misalignment 
mechanism of the same order of magnitude, as in the HHCS model. 

These uncertainties in the production of axion from strings decay prevent distinguishing 
nonstandard cosmologies with this method alone. 

One may complement the measurements of r^cDM and iria with a measurement of 
the axion CDM velocity dispersion 5v. The latter allows nonstandard cosmologies to be 
distinguished, at least in principle. The argument proceeds as follows. 

When axions start to oscillate at the temperature Ti, axions from vacuum realignment 
and axionic string decay have a momentum dispersion of the order of the Hubble scale at 
Ti [Ml], 

Spin) ~ H{Ti). (6.85) 

The momentum dispersion scales with the scale factor as Sp{T) oc l/a{T). In the standard 
cosmology, the velocity dispersion at present is then 



rria 

In the kination cosmology 



ma V a{To) J \maj ' ^ ' 



'7 



It is clear that if one has measured nia, a measurement of 5v will give the value of r^in. 
Similarly, in the LTR cosmology, 

ma V "(^o) ) \ma) ' ^ ' 

A difficulty in measuring 5v may arise from virialization of the axion population within 

galactic dark halos, although it has been claimed that 5v would be preserved in the phase- 
space evolution |161j . 
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6.7 Conclusions 

In this section we have examined the parameter regions in which the axion is 100% of 
the CDM density in cosmologies that are nonstandard before Big Bang nucleosynthesis. 
We have recognized two ways in which these regions change in going from the standard 
cosmology to the nonstandard cases. If the Peccei-Quinn symmetry breaks after the end 
of inflation (Scenario I), the axion CDM regions shift to different values of the axion mass 
rua (or of the corresponding PQ scale fa)- If the PQ symmetry breaks during inflation 
(Scenario II), the axion CDM regions can shrink or expand according to the cosmological 
model. 

We have considered two different nonstandard cosmologies that change the axion CDM 
regions in opposite directions. In the LTR cosmology, the axion CDM regions shift to lower 
axion masses in Scenario I and expand in Scenario II. In the kination cosmology, the axion 
CDM regions shift to higher axion masses in Scenario I and shrink in Scenario II. 

Different axionic string models lead to different quantitative results, but the overall 
modifications from the standard cosmology follow the same trend. 

We have also commented on the possibility to distinguish standard and nonstandard 
cosmologies using observable properties of the axion CDM population. We have tentatively 
concluded that the axion velocity dispersion may be a good indicator of the cosmology 
before Big Bang nucleosynthesis. 



CHAPTER 7 



NATURAL WARM INFLATION 

In this chapter, we turn our attention to the more generic family of axion-hke particles. 
Throughout this chapter, we consider the role of the axion-like field in the inflationary 
paradigm; in particular, we assume that the inflaton field responsible for the primordial 
inflationary period of the universe is an axion-like field. We carry out our computation in 
the context of the warm inflation scenario. 

7.1 Motivations 

In Sec. 12. 3| we have discussed the reasons that led to the formulation of the inflationary 
theory and its embedding in the standard model of cosmology. Briefly, inflation [1021 
11651 rrS} I15H [71 I118| provides a mechanism for generating the inhomogeneities observed 
in the CMBR [1371 [79| [871 11661 [IS] - and explains the observed flatness, homogeneity, and 
the lack of relic monopoles that posed severe problems in the standard Big-Bang cosmology 
[12011108] . Realistic microphysical models of inflation, in which the expansion of the universe 
is governed by the energy density of the inflaton field (j), are complicated by the requirement 
that the inflaton potential U{<j)) be very flat in order to explain the anisotropics observed 
in the CMBR [4]. 

Natural Inflation [611 [31 1152j is a viable model in which the inflaton is identified with 
an axion-like particle: in fact, the shift symmetry (j) ^ (p + const, present in axionic 
theories assures a flat inflaton potential. Although Natural Inflation is well-motivated and 
it is consistent with the WMAP measurements [152j . it is not an easy task to embed this 
model in fundamental theories like string theory [16], the main complication coming from 
the fact that the energy scale / at which the shift symmetry spontaneously breaks must be 
/ > 0.6 Mpi, in order to agree with the constraints on the scalar spectral index iig. 

In this section, we show that the energy scale / for axion-like particles can be as low 
as the Grand Unification Theory (GUT) scale Aqut ~ 10^® GeV, if Natural Inflation is 
considered in the context of warm inflation ([26l[27]; see also Refs. [89l[T35l [1211 [Mil [TTT l ): 
we refer to this as the Natural Warm Inflation (NWI) model. 
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Since axion-like particles arise in generic four-dimensional models [1261 1127^ \^3\ I44j . 
string theory compactifications |171| . and generic Kaluza-Klein theories [37j, and possess 
attractive features for inflation models like a flat potential already embedded in the theory, 
such particles have been extensively discussed in the inflation literature \W5\ \52\ \131\ \128\ 
[96 | \10 \ [95], In particular, other attempts at lowering the Natural Inflation scale / within 
the warm inflation scenario have been discussed in Refs. [131^ I130j . 

This chapter is organized as follows. We first fix our notation for the warm infiation 
scenario and for the axion particle physics in Sec. 17.21 and 17.31 respectively. In Sec. 17. 4^ 
we analyze the dynamic of the infiaton field using the slow-roll conditions and number 
of e-folds for sufficient inflation, and we study the parameter space of our NWI model. 
In Sec. 17.51 we discuss the bounds on the NWI model resulting from the measurements of 
cosmological parameters, in light of the WMAP7-|-BA0-|-SN data jllOj . Finally, discussions 
and conclusions are drawn in Sec. 17.61 

7.2 The warm inflation scenario 

In the warm infiation scenario, the infiaton field appreciably converts into relativistic 
matter (from here on referred to as "radiation") during the inflationary period. This 
mechanism is parametrized by the appearance of a dissipative term T in the dynamics 
of the inflaton fleld. In the following, we reasonably assume that radiation thermalizes on 
a time scale much shorter than 1 /T [261 [27] . 

The energy density in radiation is given by (see Sec. I2.5.3|) 

Pr = ^9.{T)T\ (7.1) 

where g*{T) is the number of relativistic degrees of freedom of radiation at temperature 
T. Here we do not specify g*(T) in the equations, but in the flgures we will always use 
g^:{T) = 228.75, corresponding to the number of relativistic degrees of freedom in the MSSM. 

Warm inflation is achieved when thermal fluctuations dominate over quantum fluctua- 
tions, or p6l[27| 

H{T) < T, (7.2) 

where H{T) is the Hubble expansion rate at temperature T. The effectiveness at which the 
inflaton converts into radiation is measured by the ratio 
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for Q ^ 1 a strongly dissipative regime is achieved, while Q < 1 represents the weak regime 
of warm inflation. Throughout this section, we will present the general equations for warm 
inflation, focusing on the case Q » 1 in the subsequent sections when specified. 

In the following, we model the inflaton field with a scalar field (j) = <j){x) minimally 
coupled to the curvature and moving in a potential U = U{(j)). The evolution of the 
inflaton field in a Friedmann-Robertson- Walker metric is described by 

'(P + {2,H + T)^ + U^ = Q, (7.4) 

where a dot indicates the derivation with respect to the cosmic time t and t/,^ = dU/d(j). 
Imposing the conservation of the total energy of the system, we find that the radiation 
energy density pr satisfies 

Pr + 4.Hpr =T<j)^, (7.5) 

with the term on the RHS of Eq. (j7.5p describing the effectiveness of conversion of the 
inflaton field into radiation. 

The total energy density of the system at any time is ptot = <P'^/'^ + U{(j),T), where 
U{(j),T) is an effective potential that accounts for temperature effects. As discussed in 
Ref. [136] . a requirement for building a consistent model of warm infiation is that finite 
temperature effects on the inflaton potential must be suppressed; this makes it possible to 
separate the effective potential of the inflaton as U {cp, T) = U {(p) + Pr{T) = U + pr- 

The total energy density during the warm inflation period is 

Ptot = ^^'^ + U + Pr, (7.6) 

and the corresponding Friedmann Eq. ()2.19p reads 

Inflation takes place when the potential U is approximately constant and dominates over 
the other forms of energy, assuring the Hubble expansion rate H to be constant. These con- 
ditions are generically imposed by requiring that the inflaton potential and the dissipation 
term satisfy slow-rolls conditions, which in the warm inflation scenario read |175j 

e<l + Q, |r?|<l + Q, |S|<1 + Q. (7.8) 

Here, we have introduced the slow-roll parameters 



109 

Notice that these conditions reduce to the usual slow-roh requirements in the cool inflation, 
where Q ^ 1. 

During the slow-roll regime, the higher derivatives in Eqs. (|7.4p and (|7.5p can be ne- 
glected. So, Eq. (j7.4p in the slow-roll regime reads 

(7.10) 

where here and in the following we use the symbol "~" for an equality that holds only in 
the slow-roll regime. Eqs. (|7.5p and (j7.7p in the slow-roll regime, respectively, read 

PrC^^<P^ (7.11) 



and 



H'^^U, (7.12) 



from which we see that a shallow potential U gives rise to a nearly constant expansion 
rate H. Notice that we have used the fact that U ^ pr, which can be combined with the 
definition of pr in Eq. (j7.ip . the requirement for a warm inflation regime in Eq. (17. 2p . and 
the Friedmann Eq. (j7.12p . to yield the constraint 

U'/' « {-^^^^y ^Pi = 5.57<,.-^/^r) X 10^« GeV. (7.13) 

This is not a stringent bound, since in most theories of grand unification the inflaton 
potential is related to the unification scale, U^^^ ~ 10^^ GeV. 

7.3 Axion-like particles 

Axion-like particles are pseudo-scalars that arise naturally whenever an approximate 
global symmetry is spontaneously broken. Examples include generic four-dimensional mod- 
els |1261ll27"t 143 1 144]. string theory compactifications |171j . and generic Kaluza-Klein theories 
|37j . Contrary to the model of the invisible axion, in which the mass and the energy 
scale fa of the invisible axion are related by Eq. (j4.12p . in a generic axion-like theory, the 
mass of the pseudo-scalar particle and the decay constant / appearing in the model 
are not linked, and are treated as two independent variables. However, it is customary to 
introduce a scale A that relates to the energy scale of the underlying theory as 



A = ^/^. (7.14) 

This definition mimics that appearing in the invisible axion theory, in which V/a is of 
the order of Aqcd- In some specific models, the value of A for pseudo-scalar particles can 
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be predicted and can range up to the GUT scale Aqut ~ 10^^ GeV or even the Planck 
scale Mpj. Here, we do not specify the value of A and we keep and / as two distinct 
parameters. 

Axion-like particles move in the potential, see Eq. (j5.12p . 



U = U{^) = A^ 



1 + cos ( y 



(7.15) 



with integer. It is well known that axion-like particles serve as suitable candidates for the 
inflaton field because the potential U{(j)) is naturally flat whenever (j)/ f <^ ir/N. Moreover, 
as we will show in Sec. l7.5.Tt the self- interaction term derived from U{(j)), 

,2 



A0 = :T7^, (7.16) 



^4 
24/2 

fulfills the requirement from observations [4], < 10"^. 

We assume that axion-like particles couple to radiation with a term similar to the axion- 
photon coupling in Sec. 14.3.21 with coupling constant gp. The decay rate of axion-like 
particles into radiation in this model is, see Sec. 14. 4. 3^ 

This dissipative term does not depend on the axion-like field (j) and we have = B = 0. 



2 m". 



r = 5^^. (7.17) 



7.4 Warm natural inflation 

Prom now on, we identify the axion-like particle (j) with the inflaton, using the flat 
potential in Eq. (|7.15p to describe the dynamics for the inflaton field in the early universe. 
Here, we will take A^ = 1 in Eq. (I7.15p . so that the potential has a unique minimum at 
(j) = vr/; the inflaton field relaxes towards this minimum in its evolution. In the following, 
we indicate the inflaton with (p, bearing in mind that we have assumed that (p is also a 
pseudo-scalar particle for which the formulas in Sec. 17.31 apply. 

Purthermore, from here on, we are interested only in the strongly dissipative regime of 
warm inflation, Q ^ 1. In this limit, the only parameters in the theory are the inflaton 
mass m^, the decay constant /, and the dissipation term P. In Table 7.1, we have collected 
the most important quantities of the theory in terms of these parameters 

We now examine the constraints on the NWI model coming from the slow-roll conditions 
in Eq. (j7.8p . the requirement for sufficient inflation, and the WMAP measurements on the 
power spectrum of density and tensor perturbations. 
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Table 7.1. Expressions for some derived quantities in the NWI theory, vahd during slow- roll 
and Q > 1. 

Quantity Equation(s) used 

=m2/2(l+cos0//), Eq. (mSI); 
U(f){4>) = —rri^f sine/)//, derived from Eq. (|7.15p : 



^-Sfv 3^(1 + cos </>//), Eqs. (ma) and dnSI); 

^2 f 

4) ~ sm^)/f, Eqs. (mil . (17121) and (I7J5]1 : 

2 ^ ^ A/ 24.nAs./n ^ Eqs. mB and dM; 



/_\/ 247r(l+cos 0//) • 
3^m|/Mpi sin^ (/.// 

■y/l+cos <}!)//' 



^ ^Jik-^^^ ^==77' Eqs. (inO|), dm]), (Hm) and dnS]); 



\/l^^7^^^^^=77' Eqs. (IZIJ, (IZIO]), (IZII]), TO 



^ = T6^(t|£w' Eqs. O and (mSI). 

^--S^F^lf^' Eqs. and (inSI). 



7.4.1 Slow-roll conditions 

During warm inflation, the values of the slow-roll parameters e, |r/| must be smaller than 
1 + Q, see Eq. (17. Sp : inflation ends when one of these two conditions is violated. In the 
following, we only consider the case in which e^l + Q~Qis violated in the strongly 
dissipating regime. Writing the parameter Q during slow-roll as 

, ^ , (7.18) 



the slow-roll condition e ^ Q reads 



1 /c/."' 



« ^====, (7.19) 



IGttG \ U J V247rGC/' 
or 

Using the expression for the potential f/(</>) in Eq. ()7.15p . this slow-roll condition gives 

sin2 0// r/ 



(l+cos</.//)3/2 ^ V 3 Mpim^' ^'^■^^^ 
Since the combination on the RHS of Eq. (|7.2ip appears frequently, we define 



a 



^jf^ = ^-^4^12 /i6 (7.22) 

3 Mpi nifj, ^ ^ 
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The slow-roll regime ends when the field (j) reaches a value (pf for which the condition 
e ^ Q is no longer satisfied, 

2 + a2 _ ^0,4 + 8q,2 



cos{(j)f/f) 



(7.23) 



We use Eq. ()7.23p definition of 4>f. 

From Eq. (j7.2ip . we see that it is always a > 1. When a ^ 1, a Taylor expansion of 
Eq. (17:23]) yields 

^ = .-^, (7.24) 

Since observations favor a large value of a > 10 (see below), here we use Eq. ()7.24p to obtain 
(j)f. 



7.4.2 Number of E-folds 

The number of e-folds is defined as 

rt2 



r 

ln(a2/ai) = / Hdt, 
Jti 



(7.25) 



where ai and a2 are the values of the scale factor a{t) appearing in the Friedmann metric 
when inflation begins and ends, respectively. Sufficient inflation requires 



iVp > 60. 



(7.26) 



During the inflationary stage, the value of the inflaton field decreases from the initial value 
to the value at the end of inflation (j)f defined via Eq. (j7.23p . We now derive a relation 
between (pi and 0/, using the definition of above. In the case of a slow-rolling of the 
inflaton and in the strongly dissipative regime, Eq. ()7.25p reads 



HT 



(7.27) 



where we used Eq. (j7.10p with F S> 2>H. Using the axion-like potential in Eq. ()7.15p and 
the Friedmann equation, we obtain 



+ cos X a , X 

Nf. = \ T—— / ax = —j= m tan — 

V 3 m^Mpi J^^if smx ^J2 4 



</-/// 



or 



tan = tan — Exp 
4/ 4/ 



V2iVe 



a 



(7.28) 



(7.29) 



113 



In general, we use Eqs. (|7.24p and (j7.29p to obtain the value of (pi/ f, 



— = 4arctan 



1 — tan 



1 + tan 



V2a 



V2Ne 



a 



(7.30) 



The appearance of trigonometric functions is due to the shape of the potential U{(j)), that 
differs from a pure quadratic one and contains a cosine function itself. 

In Figure [TTH we show the value of (j)i/f given in Eq. (j7.30p as a function of a for different 
values of the number of e- folds Nf. Blue dotted is Ng = 40, red dashed is Ng = 60, and 
green dot-dashed is Nf, = 80. Also outlined is the value (pj/f = vr as a black dashed line. 
Considering the three parameters m,^, /, and T from which a depends, moving towards 
greater values of a corresponds to a decrease in or to an increase in / or F, once the 
other two parameters have been fixed. 

For larger values of a, the initial angle (f)i/ f approaches vr and the axion-like potential 
is not distinguishable from a pure quadratic one. This fact is shown in Fig. 17.21 where we 
compare the value of the angle with Ng = 60 (black dashed line), together with the 
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Figure 7.1. The initial value of the angle f as a function of a = \/ 32tt /3 F/ /?n^ Afpi for 
different values of the number of e- folds Ng. Blue dotted: = 40; Red dashed: = 60; 
Green dot-dashed: A^e = 80. Also shown is the line (pi = irf (Black dashed line). 
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Figure 7.2. The initial value of the angle f as a function of q = y^327r/3Tf /m^Mpi for 
the case of the axion-like potential in Eq. (|7.15p (black dotted line) and for the quadratic 
potential in Eq. (j7.3ip that approximates Eq. (j7.15p around (pi ^ tt f (black dot-dashed 
line). For these lines, we fixed A'^e = 60. Also shown is the horizontal line corresponding to 

4>f/f = TT. 



value of (pi if the inflaton potential were a pure quadratic one (black dot-dashed line). The 
quadratic potential used for the plot in Fig. 17.21 in place of U {(p) in Eq. (j7.15p . is 

For a > A/'e, Eq. (j7.30p is approximated by 

(pi^TT- VSNe/a, (7.32) 
and the dynamics of the field can no longer discern between the two potentials U {(p) and 

f^quad (</>)• 

We now discuss the constraints on the parameter space of the NWI theory, based on the 
inequalities used only. 
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7.4.3 Parameter space of the NWI 

We have derived our results in the strongly dissipating regime F ^ 3H: this constrains 
the parameters as 

a » IGvr ( -r^ ) Vl + cos 0^. (7.33) 

This constraint is not severe for values of / around the GUT scale, which is the region of 
interest here. In fact, for / ~ 10^^ GeV a numerical solution of Eq. ()7.33p is 

a > 4.8 X 10"^ f^Q, (7.34) 

independently of iVg. This bound is much more loose than the requirement of flatness 
of U{(j)), which yielded a > 1. When / ~ Mp\, the constraint in Eq. (j7.34p starts 
depending mildly on the number of e-folds: for example, setting / = Mpi in Eq. (j7.33p 
gives a > 20N^-'^^. 

7.5 Perturbations from inflation 

In this section, we specialize the tools described in Sec. 12.41 to the NWI model. In 
particular, we will consider the scalar spectrum A^(fc), see Eq. (|2.5ip . the scalar spectral 
index ns{k), see Eq. (j2.59p . and the tensor perturbations spectrum A^(fc), see Eq. (j2.60p . 
using results from the warm inflation literature to set constraints on the NWI model. 

7.5.1 Scalar power spectrum 

The spectrum of the adiabatic density perturbations generated by inflation is specified 
by the power spectrum A|j(fc), given in Eq. (j2.5ip . and a scalar spectral index Ug, given in 
Eq. (I2.59p . In both warm and cool inflation models, the scalar power spectrum is generically 
expressed by Eq. (I2.53p . 

A?e(fc)=(|) iM'), (7.35) 

with (f) ~ —U(f)/{3H + r). In warm inflation, the spectrum of fluctuations in the inflaton 
field {5(p) has been computed in Ref. [20] , 

Using Eqs. (|7.35p and (j7.36p for the variance of fluctuations, we obtain the scalar power 
spectrum in the strongly dissipative regime of warm inflation as 



thermal 



A2 (,^- ^ H^/^ry^T 
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which is the same result as in Refs. [1751 1136] . with an extra factor l/(27r)^ that accounts 
for our normaUzation of the power spectrum. For comparison, the power spectrum in the 
usual cool inflation is given in Eq. (j2.56p : density perturbations are larger in warm inflation 
by a factor 



From now on, scalar perturbations are considered only in the warm inflation scenario, so 
we suppress the index "warm" in Eq. ()7.37|) . Using the expressions for H and T in Table 7.1 
we write Eq. (j7.37p as 

3/2 

(7.39) 



In deriving this last expression, we have used the fact that the largest density perturbations 
are produced when (p = (pi j61j . Eq. (j7.39p defines the power spectrum in terms of a, T, and 
A'^e, the latter variable appearing implicitly in the definition of (pi. We equate the expression 
for the power spectrum in Eq. (j7.39p to the measured value from WMAP in Eq. (j2.52p to 
obtain a relation between F and a, as shown in Fig. 17.31 for case A'^e = 60. We have checked 
that different values of A'^e do not modify the curves sensibly. 




Figure 7.3. The dissipation term F as a function of a, given by Eq. (j2.52p with 
A|(/co) = (2.445 ± 0.096) x IQ-^. We used Eq. fTMD for the analytic expression of the 
scalar power spectrum. 
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From Fig. \7.'6\ the dissipation F reaches ct maxiiniiiii valuG Trnax 

at large values of a, 
found analytically by equating Eqs. ()2.52p and (I7.39P in the limit a » 1, 

3.52 ±0.09 lo f9*{T)V^^ 
F„,ax = X 10^=^ GeV . (7.40) 



V 228.75^ 

To study the strength of the quartic self-interaction in the NWI model, we use Eqs. (j7.16p 
and ()7.22p to eliminate F in Eq. (j7.39p and obtain a relation between A,^ and a, 

Using the measured value for A^(A;o) in Eq. (I2.52p . we obtain the dependence of the quartic 
self- interaction on a, as shown in Fig. 17. 4[ We see from Fig. 17.41 that the self-interaction 
is A^ < lO""*^*^ for any a, and the NWI model can easily satisfy the constraint for the quartic 
self-interaction term A,^ < 10~* obtained in Ref. [1]. When a ^ N^, Eq. (j7.4ip with the 
measured value of the scalar power spectrum is approximated by 

A, = l.2xlO-«f^V'V4)-7^V. (742) 

* 1 228.75 / Vlov Veo; ^ ^ 




Figure 7.4. The inflaton quartic self-coupling A^ as a function of the parameter a, 
Eq. (I7.41|) . The parameters A|j(fco) and A'^e are fixed as indicated in the figure. 
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7.5.2 Scalar spectral index 

The scalar spectral index describes the mild dependence of the scalar power spectrum 
on the wavenumber k, as in Eq. ()2.5ip . Here, we do not consider the spectral tilt r in 
Eq. (|2.58p . because this quantity depends on higher orders in the slow-roll parameters. 
With this approximation, the scalar spectral index in the warm inflation scenario results in 
[83] 

ns-l = ^(-^e + lv-7B] . (7.43) 



Q\ A 2' 4 
Using the expressions for e and i] in Eq. (|7.9p . and with B = 0, we find 



or, using U{(p) in Eq. ()7.15p and its derivatives at </> = (pi, together with the values of Q, e 
and ?y in Table 7.1, 

3 [Ym^Mpi 3 + cos(/)i// 1.50 3 + cos (pi/ f 

ns = I \ = 1 . (7.45) 

8V27r Tf ^l + cos(Pi/f a y/l + cos(P^/f 

Figure [731 shows the dependence of the spectral index on a as in Eq. (j7.45p . for different 
values of the number of e- folds Np. 
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Figure 7.5. The scalar spectral index n^, Eq. ()7.45p . as a function of a, for different 
values of the number of e-folds. Blue dotted: A'^e = 40; Red dashed: A'^e = 60; Brown dot- 
dashed: A'^e = 80. Also shown are the C.L. regions from the combined WMAP+BAO+SN 
measurement, Eq. (j7.46p . Yellow: 68%; violet: 95%. 
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We have compared the values of Ug with the combined WMAP +SN+BAO measurement 

Us = 0.960 ± 0.013 at 68% C.L., (7.46) 

so that the shaded yehow region in Fig. l7.5l corresponds to the 68% C.L., while the light blue 
region is the corresponding 95% C.L. region. We see that small values of Nf. ~ 50 and large 
values of a > O(IOO) are favored, although values of A'^g = 60 — 70 can be accommodated 
in the 68% C.L. region for 100 ^ a ^ 300 or in the 95%C.L. region for all values of a. 
For large values of a the scalar spectral index is essentially independent of a, while for 
a < A^e it is the dependence on Ng that vanishes. This behavior can be summed up by 
considering Eq. (j7.45p in these two limits, 



l-n, = "^"^ (7.47) 

ll.5/iV. a>W,. 



This type of dependence of on either or a has been also noticed in models of Natural 
Inflation set in the usual cool inflation scenario, see for example Eq. (12) in Ref. |152j . In 
the case in Ref. |152j . in which there is no dissipation term P, and the inflaton mass is 
adjusted so that A ^ 10^^ GeV, the result analogous to Eq. (|7.47|) is expressed in terms of 
/ instead of a. 

7.5.3 Tensor power spectrum 

Since warm inflation considers thermal fluctuations instead of quantum fluctuations to 
generate scalar perturbations, it is only density fluctuations that modify in this scenario 
while tensor perturbations show the same spectrum as in the usual cool inflation [136j . 
Defining the tensor power spectrum as 



A^(fc)^^:-4^ = A2-(fco) f , (7.48) 



where the tensor spectral index ut is assumed to be independent of k, because current 
measurement cannot constraint its scale dependence. WMAP does not constrain A^(A;o) 
directly, but rather the tensor-to-scalar ratio 

A|(to) 

which qualitatively measures the amplitude of gravitational waves per density fluctuations. 
The WMAP-I-BAO+SN measurement constrains the tensor-to-scalar ratio as |110j 

r < 0.22 at 95% C.L.. (7.50) 
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Since in warm inflation the scalar power spectrum is enhanced by the quantity k in 
Eq. (I7.38P with respect to the value in cool inflation, the tensor-to-scalar ratio and thus 
gravitational waves are reduced in the NWI model by the same amount k. Using Eq. ()7.37p 
for the scalar power spectrum and the expression for the tensor power spectrum, 

2 



the tensor-to-scalar ratio in warm inflation is 

_ 1280F (j)^ 

With the values in Table 7.1, we find 



(7.52) 



(^\"' 1/4,^, (f^\'' (sin^.//)3/^ 



"^^^^HlSoJ ^* ^^\t^MI,) (1 + cos </).//) V8 
- 1.- X 10-3 f (-^-/C, . (7.53) 



V228.75y I ] (l + cos(/.i//)^ 

In the limit a ^ Nf^^ with (pi ^ ir — \/^Ne/a, we obtain 

Using the WMAP measure of the tensor-to-scalar ratio in Eq. (j7.52p . we constrain the 
inflaton mass in the NWI model to 

,n, < 2.7 X 10- GeV (^V'* f^V"'" ( p.55) 
* V0.22/ \m J V228.75y ^ ^ 

where we used the expression for Fmax in Eq. ()7.40p . Future measurements will not 
substantially improve the bound in Eq. (j7.55p . because of the power 1/4 raising r. As 
an example, the forecast PLANCK measurement will constrain the tensor-to-scalar ratio 
by one order of magnitude with respect to WMAP, r < 0.01, thus the bound in Eq. (j7.55p 
will approximately lower by a factor of two. 

Results for the case / = 10^^ GeV are summarized in Fig. 17.61 where we show the 
constrains on the NWI model from the WMAP data in the r-Ug plane. The WMAP data 
are better in agreement with a low number of e-folds A'^p = 40 — 50, although values up to 
A'^e = 80 can be accommodated within the 95% C.L. region. 

We conclude that a value of the axion decay constant of the order of the GUT scale 
/ ~ Aqut can be easily embedded in the NWI model, making it possible to construct 
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Figure 7.6. Prediction from the NWI model and WMAP constraints in the r-ns plane 
for / = 10^6 Bi^g dotted line: = 40; Red dashed line: Ne = 60; Brown 

dot-dashed line: iVg = 80. The yellow and violet regions are the parameter spaces allowed 
by WMAP+BAO+SN at 68% and 95% C.L., respectively. 



microscopic theories of warm inflation with axion-like particles at the GUT scale. For 
this value of the axion decay constant, the expected value of r and thus the amount 
of gravitational waves that are produced with this type of inflation is extremely low. If 
gravitational waves are found with r 10~^^ or above, this model has to be abandoned 
and a Aqut valued axion decay constant / is no longer viable. 

7.6 Discussion and conclusions 

The parameter space of the inflaton field in the NWI model is broader than that in 
the usual Natural Inflation, because of the presence of the extra quantity T describing 
dissipation and of the enlarged slow-roll conditions in Eq. (j7.8|) . 

The NWI model allows us to lower the value of the axion decay constant / from the 
Planck scale resulting in Natural Inflation to the GUT scale; a ratio f /M-p\ ~ 10"^ helps 
overcoming some difficulties encountered in Natural Inflation model-building. 

A decay constant of the order of the Planck scale is still possible, since / is not bound 
from above. In the case / ~ Mpi, Natural Inflation and NWI can be distinguished 
observationally by a measurement of the tensor-to-scalar ratio r, see Fig. 17.61 This difference 
in the value of the ratio r comes from the fact that in the NWI model, the amplitude of 
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gravitational waves is suppressed by the factor k in Eq (j7.38p with respect to the standard 
cool inflation. 

Measurements of the scalar spectral index favor larger values ol a ^ N^, as shown 
in Fig. 17.51 In this region, the value of (f)i in Eq. (|7.30p can be approximated with (pi « 
7r—\/^N(,/a, see Eq. (j7.32|) . and the axion-like potential in Eq. (|7.15p cannot be distinguished 
from a pure quadratic one through the dynamic of the inflaton. In this condition, one 
would need to obtain the value of the self-interaction independently, for example from 
considering density fluctuations in the CMBR; Natural Inflation and NWI models predict a 
precise ratio between the height of the potential and the self-interaction coupling constant, 
U /^(f) = 12 whereas other inflaton models predict different values of this ratio. 

The possibility that the axion-like energy scale / is of the order of the GUT scale when 
warm inflation is considered has been already taken into account in Ref. |131] . However, we 
differ from Ref. |131j in various points: for example, we do not consider the temperature as 
an independent variable, because in warm inflation the radiation bath is thermalized and 
temperature that is expressed is linked to the radiation energy density pr as in Eq. (j7.ip . 
Moreover, here we have included a detailed analysis of the CMBR observables r, A^(A;) 
and Us that was missing in Ref. |131j . 

We have presented a model in which Natural Inflation takes place within the warm 
inflation scenario: in such a model, the decay constant / no longer ties to the Planck scale 
as in the usual Natural Inflation model, but it can be as low as the GUT scale, / ~ 10^^ GeV. 
We have shown the viability of the NWI model and its agreement with current astrophysical 
data. 
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